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ABSTRACT 

This report consists of three parts, A to C. Part A is 
discussed in this document. Parts B and C are results from a 
subcontract with Dr. Lloyd Gay from the University of Arizona and 
are presented in separate documents. Part B is entitled 
Evaluation of Surface Radiation and Energy Balance Stations on 
the Konza Prairie and Part C la entitled Evaluation of 
Atmospheric Effects on Remotely Sensed Surface Temperatures 

Four Surface Energy and Radiation Balance Systems (SERBS) 
were Installed and operated for two weeks in Kansas during July 
of 1986. During the first week a comparative study of various 
equipment was conducted by five research groups. During the 
second week, surface energy and radiation balances were 
investigated on six sites (top of a ridge, bottom land, and east, 
west, south, and north facing slopes) on the Konza Prairie 
located about 3 km south of Manhattan, KS. 

Measurements were made to allow the computation of the 
radiation components: total solar and diffuse radiation; 

reflected solar radiation; net radiation; longwave radiation 
upward and downward. In addition, measurements were made to 
allow the computation of the sensible and latent heat fluxes by 
the Bowen ratio method using differential psychrometers on 
automatic exchange mechanisms. Data were sampled at 30 a 
intervals with battery operated computer controlled data 
acquisition systems. A total of 64 sensors were monitored by 4 
separate systems for a total of 64 system days. 

This report includes a description of the experimental 
sites, data acquisition systems and sensors, data acquisition 
system operating instructions, and software used for data 
acquisition and analysis. In addition, data listings and plots 
of the energy balance components for all days and systems are 
given . 
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INTRODUCTION 


1 . 


This report contains the results of two activities conducted 
during the summer of 1986. First, a comparative study of flux 
densities measured over a stubble surface was conducted for a 
period of one week by five research groups. Second, after the 
comparative study. Dr. Gay and I studied the energy and radiation 
balances on six sites on the Konza Prairie. 

The comparison study was conducted during the period July 12 
to July 19, 1987 on the Ashland Experimental Farm located about 
10 km south of Manhattan, Kansas. The following groups 
participated : 


GROUP NUMBER AND TYPE OF EQUIPMENT USED 


Lloyd Gay 

University of Arizona 

Edward Kanemasu 
Kansas State University 

Leo Fritschen 
University of Washington 

Harrold Weaver 

U. S. Geological Survey 

Burt Tanner 
Campbell Scientific 


4-Bowen ratio 
systems 

2-Bowen ratio 
systems 

4-Bowen ratio 
systems 

2-Bowen ratio 
systems and 
systems 

2-Bowen ratio 
systems and 
systems 


energy 

balance 

energy 

balance 

energy 

balance 

energy 

2-eddy 

balance 

correlation 

energy 

2-eddy 

balance 

correlation 


The data I collected during this test are contained in this 
report. The results of the comparisons are being presented 
elsewhere . 

Following the Ashland study. Dr. Gay and I moved our Surface 
Energy and Radiation Balance Systems (SERBS) to the same sites on 
the Konza (Figure 4.1.) that we studied during June of 1985. 

The Konza Prairie about 3 km south of Manhattan, KS. Data 
collection started on July 19, 1987 and terminated on July 26, 
1987. The data collected during this period are reported hear. 
These studies, in addition to providing energy and radiation 
balance data, were also used to evaluate the systems in hot humid 
conditions with frequent and intense thunderstorms. 

The unique instrumentation, the data acquisition 
systems, the method of analyses, and the data collected during 
the studies referred to above are presented in this report. 


2. THE ASHLAND EXPERIMENTAL FARM SITES 

The experimental equipment of the various groups was arrayed 
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from west to east at 10 m intervals across the north end of the 
lysimeter field on the Ashland Experimental Farm which is about 
10 km south of Manhattan, Kansas. The field was covered with 
stubble which had been knocked down. In addition, some wheat 
sprouting was evident. The University of Washington SERBS were 
located on the western most 40 m of the equipment line. 


3. THE KONZA PRAIRIE SITES 

The locations of the SERBS on the Konza Prairie stations 
operated by the University of Washington are shown in Figure 4.1 
as N, S, E and W. Two additional stations operated by the 
University of Arizona are shown by B and T. All heights are in a 
msl . The sites were chosen to represent the greatest variability 
in the energy balance components. That is N, S, E and W 
exposures at midslopes below the limestone outcrops. The T site 
was on top of the limestone outcrop while the B site was on the 
deep soil in the valley bottom. The systems were located on the 
slopes as follows: 1, north; 7, east; 8, west; and 9, south. 


4 . THE WEATHER 

Two Intense storms occurred during the installation of the 
SERBS at the Ashland Experimental Farm. The site was under water 
for one day after the second storm. The wind shifted to the 
south and blew at 3 to 6 m a~l for the rest of the week. Air 
temperatures increased exceeding 36 °C at 15 cm above the soil 
surface on the last two days of the first week. The relative 
humidity dropped with the increasing temperature and strong 
winds. The skies were very clear. The second week generally was 
dry. There were several cloudy days and a few showers. 


5. DATA ACQUISITION 

A small, inexpensive personal computer (NEC PC-8201 A) was 
used at each site to control the psychrometer Automatic Exchange 
Mechanism (AEM) and to sample, process, and store the data via 
the data acquisition system. Model ADC-1 data acquisition 
systems (Remote Measurements Systems) were used for data 
acquisition. These systems are modified so that the 
instrumentation amplifier can be used with all 16 analog inputs 
(rather than just the first 8), and so that two separate offset 
voltages can be applied to the instrumentation amplifier input 
from an external source. The offset voltages were supplied by an 
auxiliary module attached to each system. These modules also 
supplied constant current sources, regulated output voltages and 
served as the interface between sensors and the ADC-1. Power 
consumption of the system is given in Appendix 9.6.6. 
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Figure 4.1. The location of the aurface energy balance 
stations on the Konza Prairie near Manhattan, KS. The northeast 
corner is located at 4331000mN, 706000mE on the Swede Creek, 

Kans . SE/4 Manhattan IS' Quadrangle map. The sites are about 39° 
05' N latitude and 96° 35' E longitude. 

5.1 Computer 

The computer directed the ADC-1 to sample the data channels 
30-s intervals, with digital information being passed to the 
computer via an RS-232 port. The computer also activated the AEM 
every 6 minutes to interchange the psychrometers . After 
activation, the computer delayed sampling for two minutes to 
allow the psychrometers to attain equilibrium at their new 
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locations. Under computer control, raw data were averaged at 6 
minute intervals and recorded at 30 minute intervals on the half 
hour on a cassette tape recorder (NEC PC-8200). The computer was 
programmed so that the field operator could review the 
instantaneous data (sampled at 30-s intervals) in raw form or in 
engineering units using a single a keyboard command. In 
addition, a third keyboard command displayed calculated values of 
the energy budget components, computed and updated at 6 minute 
intervals, and a fourth display contained the instantaneous, 
present 6-minute, past 6-minute and 12 minute averages of the 
temperatures and the temperature differences. 

The computer, data system, tape recorder and two small 
batteries were housed in a 40 quart Coleman cooler which was 
covered with a space blanket. This was done to keep the computer 
at a reasonably constant temperature; in addition the space 
blanket was used to keep liquid water out of the cooler. 

5.2 Data Acquisition System 

The data acquisition system was connected to the RS-232 port 
of each computer. The ADC-1 contains 16 channels for analog 
inputs, 4 channels for digital inputs and 6 output functions. 

The basic millivolt ranges were ♦ 400 (low gain) and * 20 mV 
(high gain, using the instrumentation amplifier). Schematic 
diagram of the system are shown in Appendix 9.5. 

The basic systems were modified so that two offset voltages, 
nominally 140 and 270 mV, could be applied to the instrumentation 
amplifier signal input. Four possible gain/offset combinations 
resulted, randomly accessible to any of the 16 channels: low 

gain; high gain; high gain with 140 mV offset; and high gain with 
270 mV offset. A schematic diagram of these modifications 
(Figure 9.8) is supplied as an overlay to the basic system 
schematic (Figure 9.7). 

5.3 Auxiliary Module 

Auxiliary modules were constructed for each system and 
attached underneath. These modules supplied the ADC-1 offset 
voltages through the use of one constant current source and a 
series string of precision resistors. Two other constant current 
sources supplied the various temperature sensors. A 5 V 
regulator supplied power to the data system while a 6 V regulator 
supplied power to the computer. The sensors were interfaced to 
the terminal strips supplied on the ADC-1 by means of seven plug 
connectors on the auxiliary module. Primary power was supplied 
by a 12 Vdc deep cycle RV battery. 

A schematic of the sensor interface wiring is given in 
Appendix 9.1. Wiring diagrams, detailed descriptions and 
adjustment procedures for the constant current sources and the 
offset voltages are given in Appendix 9.3.2. 

The potentiometers used for current source adjustment are 
accessible through the side of the auxiliary module. The current 
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source circuit, is described in the National Semiconductor 
Application note "LM-334 3-Terminal Adjustable Current Source”. 
The location and identification of these adjustments are shown in 
Figure 9.4. The offset voltage is adjusted with a ten turn 
potentiometer and ten turn dial located on the aide of the 
auxiliary module. Wiring diagrams, detailed descriptions and 
adjustment procedure for the constant current sources and the 
offset voltages are given in Appendix 9.3. A schematic of the 
sensor interface wiring is given in Appendix 9.1. 

5.4 System calibration 

There are two basic adjustments for the analog input section 
of the ADC-1: the analog to digital converter (A/D> reference 

voltage via trimpot R1 and the instrumentation gain adjustment 
via trimpot R50 (see the ADC-1 Owner's Manual supplied by Remote 
Measurement Systems for component identification and details) . 
Small holes were drilled in the side of the ADC-1 where the RS- 
232 connector is located to make these adjustments accessible 
with the ADC-1 completely assembled. R1 is accessible using a 
long, narrow screwdriver through the hole on the right hand side. 
RSO can be adjusted by inserting a small screwdriver into the 
flexible plastic tubing protruding slightly from the right-hand 
hole on the RS-232 connector. 

5.4.1 Pre-experimental calibration 

The ADC- 1 ' s and the offset voltages were calibrated using a 
precision potentiometric bridge with 1 microvolt resolution, and 
an absolute accuracy of ♦ 0.02X of the reading ♦ 1 digit (Electro 
Scientific Industries model 300PVB. The offset voltages were 
adjusted using potentiometers OS1 and 0S2 (Figure 9.4) to adjust 
the voltage measured by the ESI. The ADC-1 low gain was 
calibrated using potentiometer R1 of the ADC-1 and the ESI as a 
precision voltage source set to 300 mV. All systems were 
adjusted to read 3000 raw A/D units. The standard deviation of 
lO readings was 0 units. Once this was set, high gain was 
selected and calibrated with the input set at 15 mV. All systems 
were adjusted to read 3000 raw A/D units. The standard deviation 
of 10 readings was ♦ 0.7 units. _ 

5.4.2 Previous calibration and temperature coefficient 
determination 

A calibration performed after the 1984 ASCOT experiment gave 
results which agreed to within ♦_ 1 digit of the pre-experimental 
calibrations. Temperature coefficients of the combined ADC- 
1/Auxiliary module package were then determined using a 
control led-temperature chamber. To improve system stability over 
extremes of temperature, components U29, VR, Rl, R33, R34, R50, 
R51 and R52 were upgraded to components with temperature 
coefficients of < 15 ppm o<j-~l (See Appendix E-4, ADC-1 Owner's 
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Manual for part, ident.ificat.lon>. The Instrumentation amplifier 
(U34) supplied a standard system already meet this specification. 
The results of this test are shown in Table 5.1. The temperature 
coefficient (TC) was reduced significantly for the ADC-1 itself 
(Item 1). Item 2 included the TC of the current source being 
measured, while Item 3 included the TC of the constant voltage 
source used to offset the input in this case (270 mV) . The 
smaller or negligible improvement in Items 2 and 3 are due to the 
temperature coefficients of the current sources/ voltage reference 
circuits. 

The constant current sources would reduced the latter two 
TC's. Post ASCOT 84 experimental calibration of the current 
sources indicated the Rset should be 293 ohms and R1 should be 
2404 ohms. Using these values the current varied only 7.3 x 10~9 
A °C“1, or 15 ppm °C“1, over the temperature range from 0 to 55 
°C, and 15 x lO - ^ A °C“1, or 30 ppm °C“1, over the temperature 
range from -25 to 55 °C. 


Table 5.1. Data acquisition system temperature coefficients 
(TC's) in parts per million (ppm> over the range of -25 to +35 
°C. 


Component 

ADC-1 

TC ( DDR 

oc-1) 



range 

standard 

Low TC 

1 . 

Constant 10 mv input 

high 

40 

8a 

2. 

Constant current source 

high 

180 

120 

3. 

Fixed resistance 
(simulates error in 
temperature sensors) 

high 
(270 mV 
offset) 

300b 

300 b 

a 1 

digit on the ADC-1 is the 

limit of 

resolution of 

the test. 


equivalent to 8 ppm. 
b Equivalent to 0.03 °C °C“1. 


6 . SENSORS 

The instrumentation at each site consisted of: upward and 

downward facing pyranometers; pyranometer with shadow band; net 
radiometer; upward and downward facing total hemispherical 
radiometer; three soil heat flow transducer at 5 cm depth; three 
vertical soil temperature sensor, 0 to 5 cm; wind vane; three 
cup anemometer; and two psychrometers mounted on an automatic 
exchange mechanism (AEM) . With these Instruments all of the 
components of the radiation and energy balances were obtained. 
The signals from these sensors were measured and recorded by the 
battery operated data acquisition system. 

The wiring connections between the sensors and the data 
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acquisition system are shown in Appendix 9.1. 




6.1 Solar Radiation 

The pyranometers used to measure total and reflected 
radiation Kipp and Zonen CM2 pyranometers. The pyranometers were 
mounted horizontally, 150 cm above the soil surface, at all 
stations on the end of a horizontal pipe. The horizontal pipe 
was pointed true south. The Kipp sensors were mounted in an 
adapter consisting of a tee and two floor flanges <PVC> . The 
flanges were turned to accept the sensors. This mount was 
threaded on the end of the horizontal pipe. The upper sensor was 
leveled with a bubble level placed on a plastic cylinder located 
around the glass dome of the sensor. A sun shade was used with 
the upper pyranometer . 



6.2 Diffuse Radiation 

A silicon cell pyranometer manufactured by LiCor was used 
with a shadow band to measure diffuse radiation. The shadow band 
was specially fabricated in the shop to allow the sensor to be 
mounted parallel with the slope. The operation of the shadow 
band is given in Appendix 9.4. 

6.3 Net Radiation 



Net radiation was measured with high output miniature net 
radiometers (Hicroaet Systems) . The net radiometers were 
oriented to the south at a height of 150 cm above the soil 
surface. 

6.4 Total Hemispherical Radiation 

The upward total hemispherical radiation was measured with 
specially built radiometers (Micronet Systems) . The radiometer 
consists of two high output thermopiles mounted on either side of 
an aluminum heat sink. The temperature of the aluminum heat sink 
was Measured with a 100 ohm RTD. The thermopiles were protected 
with polyethylene wind shields. The total hemspherical 
radiometers were mounted at 150 cm above the soil surface. 

6.5 Soil Heat Flux Density 

Soil heat flux density at the surface consisted of the sum 
of the change in energy storage of the 0 to 5 cm layer of soil 
and the soil heat flux measured at 5 cm. The soil heat flow was 
measured at 5 cm with three high output heat flow transducer 
(Micromet Systems) . 

The change in energy storage of the 0 to 5 cm layer was 
calculated as the product of the soil heat capacity and the 
change of the mean soil temperature of the 0 to 5 cm layer. The 
mean temperature of the layer of soil above 5 cm was monitored 
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with three 10-cm platinum resistance temperature detector 
connected in series and inserted in the soil at a 45° angle. 

6.6 Air Temperature and Vapor Pressure 

Vertical temperature and vapor pressure gradients were 
measured at each station using a pair of fan-aspirated updraft 
psychrometers that could be interchanged at selected time 
intervals by means of an AEM. The bottom psychrometers were 
located 10 cm above the top of the vegetation while the upper 
psychrometer was 1 m above the lower psychrometer . Each 
psychrometer was aspirated with a small 12 Vdc fan (Micronel 
V581L) . The intake of the psychrometer was pointed to the north 
and doward facing exhausts were installed on the psychrometer to 
reduce the effect of south winds on the ventilation rate. The 
fan drew 50 mA of current and provided 530 1/min air flow. 

The psychrometer temperature sensors consisted of a 500 ohm 
platinum resistance element encased in a stainless steel tube. 
Each tube has a serial number located near the wire end of the 
tube. The four resistance elements were connected in series to a 
constant current source, as described by Fritschen and Simpson 
(1982). With this technique, the same current was flowing 
through all resistance elements. The voltage drop across each 
resistance element was determined after it was offset with the 
270 mV offset voltage. This technique increased the recording 
sensitivity to 0.006 Q C and allowed a 40 Q C temperature range. 
Ceramic wicks were used for the wet bulbs. The wicks were one 
bar low flow ceramic with a constant head water supply of one cm. 

6.7 Wind Speed and Direction 

The wind vanes and anemometers used were manufactured by R . 
M. Young. The E and N stations has shop built anemometers. 

6.8 Automatic Exchange Mechanism 

The AEM utilized was similar in principle to that described 
by Gay and Fritschen (1982). The AEM's were plastic chain driven 
with a small 12 Vdc Brevel reversing motor which drew 450 mA for 
30 s each 6 minutes. Two AEM's were designed to allow for a 
gradient distance to be adjusted from 0 to 100 cm while the other 
could be adjusted from 0 to 200 cm. 

The wiring diagram for the AEM's is given in Appendix 9.2. 

6.9 Batteries 

Three batteries were required at each station. A 12 Vdc 
deep cycle RV battery was used to power the psychrometer fans, to 
operate the AEM and to power the three constant current sources. 
Two voltage regulators were operated from the 12 V battery. A 6 
V regulator supplied the computer while a 5 V regulator supplied 
the power for the data acquisition system. A 6 V gel cell was 
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used to supply the voltage for the offset voltages. A 6 
lantern cell was used to power the tape recorder. These 
additional batteries were need because of ground loop problems 
and were large enough to last longer than the recording period. 
The 12 V battery voltage was monitored and the battery was 
replaced with a fully charged battery when the voltage dropped 
below 10 V. A 12 W solar panel was connected to the main battery 
to keep it fully charged. 

6.10 Sensor Calibration 


All sensors used were calibrated at the University of 
Washington (UW) . These sensors include: platinum resistance 

elements; radiometers; soil heat flow transducers; soil 
temperature probes and anemometers. 

6.10.1 Radiometers 

All radiometers were calibrated against an Eppley PSP 
pyranometer SN 9030D2 on the roof of Bloedel Hall. The 
pyranometers were calibrated by correlation while the net and 
total hemispherical radiometers were calibrated using the shading 
technique. 



6.10.2. Platinum Resistance thermometers 

The 500 ohm platinum resistance temperature sensors were 
calibrated by comparing the resistance of the elements encased in 
the stainless steel tubes located in a constant temperature air 
bath (Delta Design MK2300) against the bath temperature as 
indicated by a 100 ohm platinum resistance element (Laboratory 
Standard) . This was done to adjust the O °C resistance of the 
platinum elements. The universal resistance-temperature relation 
(line 1550 of SAMP. BA, Appendix 9.8.1) was used to compute other 
temperatures. Previous test indicate that the universal relation 
applies quite well (Fritschen and Simpson, 1982) . The results 
from these calibrations are shown in Table 6.1. 

6.10.3 Soil heat flow 

The soil temperature probes were calibrated in the oven 
using the PLS as a standard. 

The soil heat flow transducers were calibrated in a chamber 
at UW. The heat flow calibration chamber consists of two 
aluminum tanks (7.6 x 33.0 x 43.2 cm) spaced by 5.84 cm. The 
facing sides of the tanks consist of aluminum plates (0.635 x 
35.6 x 45.7 cm). The plates are spaced by 1.27 x 5.84 cm pieces 
of plastic creating a void between the tanks of 5.84 x 33.0 x 
43.2 cm. The void has a cross section of 1425.6 cm^ and a volume 
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Table 6.2. Platinum resistance element calibrations. Air bath 
temperature was 21. SO °C for the calibration. 


SERIAL 

NO. 

OUTPUT 

(mV) 

ZERO 

RESISTANCE 

(ohms) 

1 

270 . 80 

499.209 

2 

270 . 92 

499 . 430 

3 

270.68 

498 . 987 

4 

270 . 92 

499 . 430 

5 

271.09 

499.743 

6 

271.08 

499.725 

7 

271.19 

499.928 

8 

271.17 

499.891 

9 

271.12 

499.798 

io 

271.05 

499.669 

11 

271.15 

499.854 

12 

271.14 

499 . 835 

13 

271.09 

499.743 

14 

271.04 

499.651 

15 

271.10 

* 499.762 

16 

271.12 

499.798 

17 

271.13 

499.817 

18 

271.21 

499 . 964 

19 

271.15 

499 . 854 

20 

271.22 

499 . 983 


of 8,325.5 cm3. The plastic spacers have a total cross section 
of 145.3 cm2 and a volume of 848.5 cm3. 

The void between the tanks is filled with glass beads 10 
microns in diameter. The porosity of the beads is 36 percent. 

One junction of a thermocouple loop is cemented onto the surface 
of each aluminum plate so that the temperature difference between 
the inside surfaces of the aluminum plates (or the temperature 
difference across the glass beads) can be measured directly. 

If water of temperature, Tl, is circulated through tank one 
and water of temperature, T2, is circulated through tank two, 
then the heat flowing' from tank one to tank two <T1 > T2) can be 
expressed by: 

G = (k/IXTl - T2) 

Where G is the heat flow in W m~2^ j* j. a the thermal conductivity 
of the water glass bead mixture 0.94 W (m °K>“1 and 1 is the 
distance between the plates (5.04 cm). 

The heat flow is related to the mV signal of the transducers 
placed in the glass bead water mixture. 



7. SOFTWARE DESCRIPTION 

Two categories o£ software will be considered. The first 
consists of programs written for the NEC personal computers that 
were used for field data acquisition, data transfer, and data 
acquisition system testing. The second consists of a series of 
programs used for post-experimental data processing, including 
energy and radiation balance calculations, plots and printed 
summaries used in this report. 

Listings of the various programs are given in Appendix 9.8. 
All software is still in the development stages, and as such is 
not free from errors, nor have all the refinements been 
incorporated to make their operation “user friendly'*. This is 
especially true of the auxiliary, supporting software. However, 
based on the excellent field performance of the primary data 
acquisition and processing program SAMP. BA, it is felt that the 
software is basically sound. 

7.1 Data acquisition (SAMP. DO, INDATx.DO) 

Data acquisition and field procession was controlled by 
program SAMP. BA, a BASIC program written for the NEC portable 
computer. This program, listed in Appendix 9.8.1, is largely 
self documenting. Statements 1000 to 1192 constitute the main 
program, with control of subroutine calls routed through a jump 
table in lines 100 to 300. 

The program is customized for a particular location, system, 
and set of sensors through the use of an input file named 
“INDATx.DO", where x was an identifier set equal to the 
particular system number. Files INDAT1.D0 through INDAT9.D0 for 
the four systems used in the present experiment are listed in 
Tables 7.1 to 7.4. Control parameters set by files INDATx.DO are 
partially identified in lines 1, 3 and 5, and the last two 
columns of the actual file, or by comparison with program lines 
9110-9195 where the data is read by SAMP. BA. A decription of 
these identifiers is found in Table 7.5. 

The psychrometer separation and site elevation are included 
in the INDAT files. A standard atmosphere (101.3 kPa> is assumed 
in the calculation of atmospheric pressure (P) , which is then 
corrected for altitude using a lapse rate of -0.01055 kPa m“l. 
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Table 7.1. Input file for station 1 <W) . 

M N1 N2 N3 N4 N5 NS GO M7 

18, 30,30, 1,30, 3, 0, 2,16 

LG HG HONE REF 01 02 RC NCRTD 

9.997,200.7,7, 0,264.94,159.96, .4981466,9 

DELZ ELEV CSOIL DZ REF HOME 

1.00, 315, .27, 0.05, O, 7 


CN 

RG 

GAIN BIAS 

TYPE 

DESC. 

SER. NO. 

1 , 

2, 

30.742, 

0 , 

4 

G 

16 

2, 

0 , 

22.36, 

0 , 

4 

Q 

QS6004 

3, 

2, 

82.13, 

0 , 

4 

Kdn 

3750 

4, 

2, 

82.92, 

0 , 

4 

Kup 

3701 

5, 

2, 

153.39, 

o. 

4 

D 

1577 

6, 

0 , 

11.85, 

0 , 

4 

Qdn 

THR86004 

7, 

2, 

1.00, 

0 , 

6 

Home 


8, 

o. 

3.488, 

o. 

3 

DIR 

1 GILL 

9, 

1, 

499.96, 

0 , 

2 

TAT 

11 3/4 

10, 

1, 

500.00, 

0 , 

2 

TWB 

11 5/6 

11, 

1, 

500.00, 

0 , 

2 

TAT 

12 3/4 

12, 

1, 

500.01, 

0 , 

2 

TWT 

12 5/6 

13, 

3, 

301.62, 

0 , 

2 

TS 

3000HN 

14, 

0 , 

100, 

0 , 

2 

Tthr 

86004 

15, 

0 , 

0.2027, 

0.7, 

4 

U 

GILL1 

16, 

0 , 

12.28, 

0 , 

4 

Qup 

THR86004 

INDAT7 

.DO 


15:33 

1/16/87 

BASED ON INDAT8 

.DO 

OF 5/17/86 
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Table 7. 

2. 

Input file 

for station 7 

' (S) . 

M N1 N2 N3 N4 N5 NS GO M7 



18, 

6, 

30, 1,30 

3, 

2, 2,16 



LG 

HG HOME 

REF 

01 02 RC NCRTD 

10. 

0,193.3,7, 

o. 

264.67 

,160.54, .49829,9 

DEL2 

ELEV CSOIL 

DZ REF HOME 

l.OO, 

315, .15 

f 

0.05, 

0, 7 

t 

CN 

RG 

GAIN BIAS 

TYPE 

DESC. 

SER. NO. 

1, 

2, 

30.293, 

0, 

4 

G 

76 

2, 

0, 

22.81, 

0, 

4 

0 

Q86021 

3, 

2, 

80.88, 

O, 

4 

Kdn 

773743 

4, 

2, 

80.45, 

0, 

4 

Kup 

773741 

5, 

2, 

88.90, 

0, 

4 

D 

6712 

6, 

o. 

16.63, 

0, 

4 

Qdn 

PYR86002 

7, 

2, 

1.00, 

0, 

6 

Home 


8, 

0, 

3.837, 

0, 

3 

DIR 

73 GILL 

9, 

1, 

499.89, 

0, 

2 

TAB 

71 3/4 

10, 

1, 

500.10, 

0, 

2 

TWB 

71 5/6 

11, 

1, 

500.03, 

0, 

2 

TAT 

72 3/4 

12, 

1, 

499.97, 

0, 

2 

TWT 

72 5/6 

• 13, 

3, 

302.22, 

0, 

2 

TS 

3000HM 

14, 

0, 

1O0, 

0, 

2 

TO 

PYR86002 

15, 

0, 

0.1924 

, 0 

.7, 4 

U 

73 GILL 

16, 

0, 

15.19, 

0, 

4 

QUP 

PYR86002 

INDAT7 

.DO 


1541 11/04/86 

BASED 

ON INDAT8 

.DO 

OF 5/17/86 













Table 7.3. Input file for station 8 (E> . 

M N 1 N2 N3 N4 N5 N8 GO M7 

18, 6,30, 1,30 3, 2, 2,16 

LG HG HOME REF 01 02 RC NCRTD 

10.0,199.9,7, 0,264.55,160.40, .49877,9 

DELZ ELEV CSOIL DZ REF HOME 


l.OC 

* 9 

315, .27 

9 

0.05, 

0, 7 


CN 

RG 

GAIN BIAS 

TYPE 

DESC. 

SER. NO. 

1* 

2, 

39.03, 

o. 

4 

G 

86 

2, 

o. 

22.25, 

0, 

4 

Q 

086022 

3, 

2, 

84.36, 

0, 

4 

Kdn 

838771 

4, 

2, 

85.09, 

0, 

4 

Kup 

838751 

5, 

2, 

92.49, 

0, 

4 

D 

6710 

6, 

0, 

12.66, 

0, 

4 

Qdn 

Thr86003 

7, 

2, 

1.00, 

0, 

6 

Home 


8, 

0, 

1.340, 

0, 

3 

DIR 

83 GILL 

9, 

1, 

499.65, 

0, 

2 

TAB 

81 3/4 

10, 

1, 

499.71, 

0, 

2 

TWB 

81 5/6 

11, 

1, 

500.31, 

0, 

2 

TAT 

82 3/4 

12, 

1, 

500.23, 

0, 

2 

TWT 

82 5/6 

13, 

3, 

302.58, 

0, 

2 

TS 

86 

14, 

0, 

lOO, 

0, 

2 

TQ 

THR 86001 

15, 

0, 

0.2473 

,0. 

3, 4 

U 

83 FRIT 

16, 

0, 

12.26, 

0, 

4 

QUP 

THR86003 

INDAT7 

.DO 



1550 11/04/86 


BASED ON INDAT8 . DO OF 5/17/86 


Table 7.4. Input file for station 9 (N) 

M N1 N2 N3 N4 N5 NS GO M7 
18, 6,30, 1,30 3, 2, 2,16 
LG HG HOME REF 01 02 RC NCRTD 

10.0. 200.0.7, 0,264.67,159.96, .49909,9 

DELZ ELEV CSOIL DZ REF HOME 

1.00, 315, .27, 0.05, 0, 7 


CN 

RG 

GAIN BIAS 

TYPE 

DESC. 

SER. NO. 

1, 

2, 

30.516, 

0, 

4 

G 

96 

2, 

o. 

22.00, 

0, 

4 

Q 

Q86023 

3, 

2, 

82.15, 

o. 

4 

Kdn 

001 

4, 

2, 

82.55, 

0, 

4 

Kup 

60294 

5, 

2, 

167.90, 

0, 

4 

D 

1579 

6, 

0, 

12.89, 

0, 

4 

Qdn 

Thr6001 

7, 

2, 

1.00, 

o. 

6 

Home 


8, 

0, 

1.10, 

o. 

3 

DIR 

93 GILL 

9, 

1 , 

500.10, 

o. 

2 

TAB 

71 3/4 

10, 

1, 

499.93, 

o. 

2 

TWB 

71 5/6 

11, 

1, 

499.88, 

o. 

2 

TAT 

72 3/4 

12, 

1, 

500 . 06 , 

o. 

2 

TWT 

72 5/6 

13, 

3, 

301.38, 

0, 

2 

TS 

96 

14, 

o. 

10O, 

o. 

2 

TO 

THR 86001 

15, 

o. 

0 . 0905 

, o. 

3, 4 

U 

93 FRIT 

16, 

o. 

12.24, 

o. 

4 

QUP 

THR86001 


INDAT7 . DO 1554 11/04/86 

BASED ON INDAT8 . DO OF 5/17/86 

Table 7.5. Description of INDATx.DO control files used in 

program SAMP. BA. 

M Total number of variables in each data record 

N1 Minutes in each averaging period (between Bowen ratio 

interchanges) 

N2 Number of seconds between samples <0 < N2 < 59) 

N3 Maximum number of records allowed in memory storage 

buffer (calculated in program line 9265) 

N4 Minutes between data output to cassette tape (Changed 

in the program into N4/N1, which is the number of data 
records written to the cassette each access. 

N5 Number of times each analog channel is sampled before 
the value is saved. This allows for a longer settling 
time for the A/D converter when sampling low level 
signals using the on board amplifier. 

N8 Number of minutes samples are not taken after the Bowen 

ratio interchange device has operated to allow temperatures 
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to come into equilibrium. 

GO Not used 

M7 Total number o£ analog and digital inputs being sampled 
(M7 = M - 2; 2 variables are used to store date and 
time) . 

LG Gain of low range <mv/AD unit) 

HG Gain of high range (mv/AD unit) . Selected by adding 32 
to the channel number . 

HOME Channel number of AEM Home signal. 

REF Channel thermocouple reference connected to (not currently 
in use) 

01 Offset #1 (nv) . Selected by adding 16 to the channel 
number. 

02 Offset #2 (»v) . Selected by adding 48 to the channel 
number. 

RC<1> Value of constant current through dry and wet bulb 
resistance temperature elements (ma) . 

NCRTD Channel number of the first resistance temperature 
element 


CHAN C(K) 
RANGE Cl (K) 


GAIN GOO 
BIAS B (K ) 
TYPE N (K) 


DESC. X$ 


Array of channel numbers 


0 

= Lo 

gain 


-adding 

0 

to ( 

chan. no. 

1 

= Hi 

gain. 

offset 1 

-adding 

16 

to 

chan, no 

2 

= Hi 

gain. 

no offset 

-adding 

32 

to 

chan, no 

3 

= Hi 

gain. 

offset 2 

-adding 

48 

to 

chan, no 

mv 

gain 

(eng . 

units/mv ) 






bias (eng. units) 


1 * type K thermocouple 

2 = resistance temperature element 

3 = wind direction 

4 = linear calibration 

5 = digital input 

6 = Home signal 


Used in data file for description only 


In addition, the following quantities are calculated in 
connection with the above control parameters. 


16 


G2(K) mv gain for each channel 


B1(K) offset for each channel (if used, otherwise zero) 

C1(K> This is converted to the actual channel number plus the 
offset for use in the A/D routine 

NRTD The number of resistance temperature elements 

NWD not used 

NDIG number of digital channels 
NANLG number of analog channels 

7.2 Data transfer (READT2.BA) 

The procedure used for reading cassette data tapes in the 
field used a combination of the ROM-resident communications pro- 
gram called TELCOM, and the BASIC program READT2.BA. TELCOM was 
used to initially set up the receiving computer, in this case a 
AT compatible microcomputer using DOS operating system. The NEC 
computer is configured as a terminal, and connected to the RS 232 
port on the AT. Cross Talk software was used on the AT to 
capture the transmitted files. READT2.BA was then used to 
transmit the data. 

7.3 Test Programs (ADCTST.BA) 

A program was developed for use in testing the operation 
of the data acquisition system, ADCTST.BA. The ADCTST.BA 
(Appendix 9.8.6) uses the built-in serial port driver. 
Communication with the ADC-1 from a BASIC program via the 
standard serial port driver uses INP and OUT statements. 

7.4 Post Experimental Data Processing and Data Conversion From 
Raw to Engineering Units (SAMPE.BAS) 

The second series of programs were developed for initial 
post-experimental data processing, including energy and radiation 
balance calculations, plots and printed summaries used in this 
report. They are coded in Microsoft BASIC 5.2, and were intended 
to be compiled and run using the Microsoft BASIC compiler to 
reduce execution time. The data conversion (SAMPE.BAS) and 
energy/radiation balance processing (SAMPP.BAS) programs are 
based in part on the field sampling and analysis program SAMP. BA 
(Section 7.1.1). All programs are controlled by an input file 
named PDS.FIL which contains values of certain control para- 
meters, an identification label, and a list of file names to be 
processed (Table 7.6). The meaning of the control parameters 
varies, depending on exact program involved. 










1 

1 

4 

I 




1 

1 

1 


Table 7.6. Sample contents of control file PDS.FIL. 

0,4,0, B: ,D: , .MF,P 

Energy balance, 6 minute data <SAMPP 

1/12/85 > SI 0929 , S50929, S30929G , S40929 , S40S29T , S20929, END 

The ASCII raw data files were edited using a text editor so 
that each contained one day's data for one data system, starting 
and ending at 0000 hours. Since an average was stored every 6 
minutes during data collection, each file contains a maximum of 
241 records. Using this data as input, SAMPE.BAS (Appendix 
9.8.8) converts the raw data (in A/D units) to engineering units 
(e.g. oC, m s-1, etc.). System and time specific data is found 
in lines 9300 - 9400, lines 6300 - 6400 (analogous to lines 9000 
- 9106 and 9110 to 9195 in SAMP. BA) and lines 6200 - 6300. The 
identical files used in the field analysis (Tables 7.1 to 7.5) 
are used with this program. The data is stored in a BASIC random 
file in compressed binary format. 

7.5 Energy balance processing (SAMPP.BAS) 

SAMPP.BAS (Appendix 9.8.9) uses the output from SAMPE.BAS 
to compute the radiation and energy balance. The same basic 
calculational algorithms as those in the field programs are used, 
with some minor additions. The most significant of these was the 
inclusion in the soil heat flux term of the energy storage in the 
layer of soil above the heat flow transducer. Others include the 
setting to zero of small amounts of spurious negative shortwave 
radiation which occasionally occurred at night, the shadow band 
correction for the diffuse radiometer, the processing of the AEM 
"home** signal, and a diagnostic routine in development to detect 
a drying wet bulb. These changes should be considered in any 
reanalyses of the data. 

7.6 Data Plotting and 30 Minute Summary Listings 

The resulting radiation and energy balance data were then 
averaged for 30 minute periods, and the results plotted were 
plotted and/or printed in summary listings. Examples of programs 
used for this purpose are included in Appendix 9.8.10. 
SUMMARYE.BAS prints a table of 30 minute averages and totals for 
a 24 hour period. As input it uses either the 6 minute data 
created by SAMPP.BAS, in which case the 30 minute averages it 
creates are stored on a summary disk file, or it uses the 30 
minute average created by a prior run of SUMMARYE.BAS to print 
the table only. This choice is determined from the value of the 
flag ICFLG (the first parameter in file PDS.FIL; Table 7.7). 
Sample programs which produce line printer plots of the 
radiation and energy balance data are reproduced in Appendices 
9.8.11 and 9.8.12) 

A spreadsheet called SMART was used to create the 30 minute 
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average files and the graphs included in this report. 

8. RESULTS 

8.1 History of Data Records 

During the first week some equipment difficulty was 
experienced. They included: the data system board shorting 

against the chassic on system 1; a wire shorting the exchange 
mechanism on system 7 which burned out the input switch for the 
home signal; one of the fans stuck some time during the week on 
system 9; and a tape recorder battery connected in reverse 
during a battery change. These problems were solved by the end 
of the first week and all stations were operating properly during 
the second week. A rodent cut the wind vane wire on the south 
slope system some time during the second week. 

8.2 Energy and Radiation Balance Data Listings 

Listings of the 30 minute averages, and various totals or 
averages of the energy and radiation balance data plus other 
environmental data are present in Attachment 1 for the first week 
and in Attachment 2 for the second week. The time represents the 
preceding 30 minute period. All times are CDST. 

8.3 Plots of Energy Balance Data 

Plots of the 30 minute average energy balance data are 
present in Attachment 3 for the first week and in Attachment 4 
for the second week. The averages are plotted at the time 
representing the preceding 30 minute period. All times are CDST. 


9. APPENDICES 

9.1 Sensor Wiring Diagram 

The wiring of all sensors to and through the current housing 
box located below the data system to the data system is shown on 
the following page. The wires are color coded as follows: 

B, black 
Br, brown 
Gr, green 
Gy, gray 
0, orange 
R , red 
T, tan 
W, white 
Y, yellow 

Other codes or symbols are defined as: 










1 

« 


« 

I 

I 

I 

I 

I 

I 

I 

I 

I 

1 

I 

1 


Q», net radiation 

THR, total hemispherical radiometer 
THRT, temperature sensor of THR 
K , global solar radiation 
K , reflected solar radiation 
cur, current 


Figure 9.1. Energy balance station sensor wiring diagram (see 
following page) 
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9.2 AEM Wiring Diagram 


The wiring for the automatic exchange mechanism is shown on 
the following page. The wires are color coded as follows: 


B, 

black 

Br, 

brown 

Bu , 

blue 

Gr, 

green 

Gy, 

gray 

0, 

orange 

P, 

pink 

R, 

red 

T, 

tan 

w. 

white 

Y, 

yellow 


Other symbols used are defined as: 

P#, plug 
R#, relay 
S#, Switch 

Switch 1 is used to turn on the power to the AEM. Relay 2 
controls the direction of rotation of the drive motor by 
controlling the polarity of the applied voltage with S7 and SS. 
Relay 2 can be accuated automatically from the data system 
through R1 and S10 or manually with momentary S2. Switches 4 and 
5 limit the length of travel by interrupting the power to the 
drive motor. A line fuse (2 amp) is located in the power line to 
the drive motor. 

Switches 3 and 6, called the home switches, indicate the 
position of the AEM by reversing the polarity of the voltage drop 
across a 30 ohm (some AEM's have 10 ohm) resistor. Current 1 is 
applied to this resistor via pins 5 and 6 while the voltage is 
sensed via pins 1 and 2 of DE9P #7. The convention used is the 
signal is + when the right hand psychrometer (looking from the 
drive motor end) is down, - when it is up and 0 when neither 
psychrometer is in home position. NOTE: BE SURE TO CONNECT THE 

12 VDC PROPERLY. THE RED WIRE IS ♦ AND BLACK WIRE IS -. THE AEM 
WILL NOT FUNCTION AND THE FAN WILL RUN BACKWARDS IF THE VOLTAGE 
IS REVERSED. 


9.3 Auxiliary Module Description and Operation 

9.3.1 Current Source and Offset Voltage Wiring Diagram 

The wiring of the current sources and the offset voltages is 
shown on the following page. The current sources and offset 
voltages are located on a board in the box under the data system. 
The wires are color coded as follows: 
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B, 

black 

Br, 

brown 

Gr, 

green 

Gy, 

gray 

0, 

orange 

R, 

red 

T, 

tan 

w. 

white 

Y, 

yellow 


Three 0.5 mA current, sources were utilized with each data 
system. Current source 1 supplied for the RTD's and the home 
signal. The output of this source is the brown wire attached to 
pin 24 which is connected to socket 1 pin 1. The return line is 
a tan wire attached from socket 7 pin 6 and is attached to pin 
23. Pins 1 and 2 are tie points to monitor current 1. 

Current source 2 supplied the current for the wind vane and 
the vertical soil temperature probe. Its output (pin 22) is a 
yellow wire which is attached to socket 3 pin 3. The return line 
is a tan wire from socket 6 pin 2 and is attached to pin 21. 

Pins 3 and 4 are tie points to monitor current 2. 

Current source 3 supplied the current for the offset voltages. 

The offset voltages are terminated on pins 16, 17, and 18. Pin 

16 (red wire) is the plus side of the small offset voltage, pin 

17 is the negative side. Pin 18 (white wire) is the plus side of 

the large offset voltage and pin 17 is the negative side of this 
offset voltage. The offset voltages can be monitored via pins 9 
and lO which are attached to channel 15. Voltage for current 
sources 1 and 2 was supplied by a 12 Vdc battery attached to pins 

7 and 8 via a black jacketed cable. Voltage for the offset 

voltages was supplied through pins 11 and 12 via a orange 
jacketed cable attached to a 6 Vdc battery. Two voltage 
regulators are also located on this board. The 5 Vdc regulator 
supplied the voltage for the data acquisition system via pins 20 

and IS. The 6 Vdc regulator supplied voltage to the computer via 

pins 5 and 6 which have a gray jacketed wire attached to them. 

Two additional tie points provided voltage to the anemometer 
(socket 3 pins 7 and 8) and to the temperature sensor in the 
total hemispherical radiometer (socket 5 pins 3 and 4) . 

9.3.2 Current Source and Offset Voltage Adjustment 

The current source circuit is described in the National 
Semiconductor Application note LM- 334 3-Terminal Adjustable 
Current Source. Each current source as two adjustments, Rset and 
R1 . Rset is initially set using an ohmmeter to 285 ohms. This 
value is based on minimizing the temperature coefficient of the 
current source. R1 is then adjusted for each (about 2404 ohms) 
so that the current is 0.500 ma . 
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The nominal magnitude of the offset voltages are set by the 
precision resistors in series with potentiometers R3 and R4 . 

These voltages can be adjusted +12.5 mV using these 
potentiometers . 

9.4 Operation of the Shadow Band 
9.4.1 Mounting the Shadow Band 

The procedures for the operation of the shadow band are 
adapted from Publication 8105-22 from Li-Cor, Lincoln, NE. 

First one has to mount the base of the shadow band on a 
horizontal pipe with the clamp supplied. The horizontal pipe 
must be oriented true north and south. A slight error in 
orientation will cause the sun to shine on the sensor during the 
day. Next level the base using the three socket screws in the 
base. After the base is leveled lock it in place with the three 
bolts in the base. 

The rod attaching the shadow band to the base is now 
adjusted for the latitude of the site. Latitude marks are 
scribed on the vertical curved surface attached to the base. The 
upper portion of the rod should be set at the correct latitude. 

The sensor platform is mounted on a gimbal and should be 
leveled in both directions. After the sensor platform has been 
leveled, the shadow band should be adjusted for the solar 
declination (See 9.4.2). The upper surface of the block 
containing the thumb screw should be set on the proper 
declination angle which is scribed on the upper surface of the 
rod. The declinations are marked every 5 degrees and labeled 
every 10 degrees. After these adjustments have been made, the 
shadow band should cast a shadow over the sensing element. If 
this is not the case, check the latitude and declination angles. 
Also make sure that the horizontal pipe and rod are oriented true 
north and south. This can be done by observing the shadow cast 
be the pipe and rod. At solar noon (See 9.4.3) the shadow should 
be directly below the rod or pipe. A plumb bob attached to a 
piece of string will help to determine the location of the 
shadow. 


9.4.2 Shadow Band Adjustment 

9. 4. 2.1. Declination angle versus time of year. 

The declination angle can be approximated as follows: 

= sin C (360°/365 days) <n-n V e^ 

where D = declination angle, E = obliquity of the ecliptic 
(=23.45°, n = day of year and n V e = day °* vernal equinox =81. 
Then 

= 23.45 sin C0.9863(n - 81)3 
or 
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23.45 sin [0.9863(n + 284)3 


Example 1: 

Date: March 16, 1981 <n = 75) 

= 23.45 sin [0.9863(75 + 284)3 

= 23.45 sin [0.9863(359)3 

a 23.45 sin 354° 

= -2.40 


9. 4. 2. 2 Shadow band adjustment versus time of year and latitude 

The shadow band projects a shadow of varying widths upon the 
sensor diffusing eye depending upon the time of year and 
latitude. In addition, the declination angle of the sum changes 
at different rates depending on the time of year. Both phenomena 
need to be considered when determining how ofter to adjust the 
band . 


- Declination angle - 23.45 sin CO. 9863 (n ■*■284)1 

Rate of declination angle change 

d /dn = 23.45 cos CO. 9863 <n +284)3 (0.9863) 2Prad/360° 

= 0.404 cos [0.9863 (n +284)3 degrees/day 

The maximum rate of change is 0.404 degrees/day and occurs 
when n = 81 (vernal equinox) and n = 263 (autumnal equinox) . The 
minimum rate of change is 0 degrees/day and occurs when n = 172 
(summer solstice) and n = 355 (winter solstice) . 

The angle subtended by the shadow band on the sensor 
diffusing eye is 

F = 2 tan-1 C(w cos D-d cos (F -D)/(2r/cos D)3 
for (F - D) < 90° 

where D = declination angle, w = width of shadow band, d = sensor 
diffusing eye diameter, r = band radius and F = latitude. 

For a perfectly aligned shadow band, one should change the 
band as follows: 

1. Determine the day of the year (n). 

2. Determine the latitude (F). 

3. Determine the declination angle (D) . 

4. Determine the rate of declination angle change dD/dn. 

5. Determine the angle subtended by the band (F) . 

6. Determine the number of days shadowed (Ds) as follows: 

Da = (F - Fs) / (dD/dn) 

where Fs = angle subtended by the sun (- 0.5°). However, let Fs 
= 1° as a safety factor. 
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Example 2: 


Date: March 16, 1981 <n = 75) 

2. Q = 41^ 

3. = 23.45 sin CO. 9863(75 + 284)] ■ -2.42° 

4. dD/dn = 0.404 cos CO. 9863 <75 ■•■284)3 = 0.402 degreea/day 

5. F = 2 tan-1 C<0.5 cos (-2.42) - 0.307 cos <41-<2.42)3 

< (2) (3) /cos (-2.42)) 

F * 5.27° 

6. Ds = 5.27° - 1° = 10 days (for a perfectly 

0.402 deg/day aligned band) 

A more realistic value for Ds in actual operation would be about 
1/2 of Ds or 5 days. 

9.4.3 Equation of Time and Time of Solar Noon 

Example: Compute the time of solar noon at longitude 81° 

38 ' West on September 22, 1980. 


1. Determine the day of year (n). 
September 22, 1980 is day number 266. 

2. Determine West longitude in hours (H) . 


H = 81° 38' West = 81.633° X 24 hours 
= 5.4422 hours = 5 hours 27 min 


360° 


3. Determine time elapsed (t) in days since January 0, 
0 hour UT 

t = n ♦ (UT ♦ H)/24 


where UT = universal time = 12 for solar noon 

t = 266 + (12 ♦ 5.4422) /24 = 266.72676 

4. Determine equation of time (EOT) 

EOT = -7.64 sin <0.9893t) ♦ 0.56 cos <0.9863t) 

-9.37 sin C2 (0 . 9863t ) 3 - 2.83 cos C2(0.9863t)3 min. 

EOT = -7.64 sin (0.9893(266.73)) ♦ 0.56 cos 
(0.9863(266.73) ) 

-9.37 sin C2<0. 9863(266. 73) )3 - 2.83 cos 
C863 <266 . 73) ) 3 min. 

= 8.0 min. 


The above example uses the EQT for the year 1980. This 
should provide adequate north-south alignment for any year. The 
current EQT can be obtained from the Almanac for Computers, 1980, 
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34th and Massachusetts 


Nautical Almanac for computers, 1980, 

Avenue, N. W. , Washington, DC 20390. 

5. Determine local mean time (LMT) 

LMT = 12 h OO min - 8.0 min * 11 h 52 min 

6. Determine universal time <UT) 

UT = LMT + H = 11 h 52 min + 5 h 27 min = 17 h 19 min 

7. Determine local time (LT> 

LT * UT - dT 

where dT is the difference in time zones between H and Greenwich, 
England. 

= 17 h 19 min n (Eastern Daylit Time) 

The time of solar noon on September 22, 1980, at a longitude of 
810 3Q' West is 1:19 pm (EDT) . 

9.4.4 Data Reduction 

Use of a shadow band necessitates applying a correction 
factor to the data to allow for that part of the total diffuse 
radiation which is obstructed by the band. In addition, a 
correction may be necessary if the spectral response of the 
sensor ia not ideal due to the variation of spectral irradiance 
between blue sky and various cloud conditions (as in the case 
only of the LI-200SB pyranometer ) . 

The problem of correcting the data should be approached both 
theoretically and experimentally, although neither approach is 
entirely satisf actorly in itself because the diffuse radiation 
varies over the dome of the sky (International Energy Agency, 
1980). 

Table 9.1 contains theoretically derived correction factors 
for the band obstruction for isotropic sky conditions on the 16th 
of each month. An additional 4X additive correction is included 
in the table values to account for the effects of non-isotropic 
distribution of the radiance over the sky. It should be realized 
that these corrections are approximations for general sky 
conditions and are not a substitute for corrections derived 
experimentally at a given location. The measured values of 
diffuse radiation should be multiplied by the appropriate 
correction factor. 

The correction factor for clear sky conditions can be 
determined experimentally by comparing the diffuse measurement 
(as measured when the shaddow band is in its normal position), to 
the diffuse measurement when a shadow disk is used to shadow the 
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sensor instead of the sadow band. The difference between these 
two measurements is the portion of diffuse radiation that is 
obstructed by the shadow band (International Energy Agency, 
1980). 

Table 9.1. Correction factors for solar radiation obstructed 
by the shadow band. 


Lat . °N Jan 
Lat . °S Jul 

Feb 

Auq 

Mar 

Sep 

Apr 

Oct 

May 

Nov 

Jun 

Dec 

Jul 

Jan 

Aug 

Feb 

Sep 

Mar 

Oct 

Apr 

Nov 

Mav 

Dec 

Jun 

0 1.12 

1.15 

1.17 

1.15 

1.14 

1.12 

1 . 12 

1.14 

1.16 

1.15 

1 . 14 

1 . 12 

10 1.11 

1.14 

1.16 

1.16 

1.14 

1.13 

1.14 

1.15 

1 . 16 

1.14 

1 .12 

1.11 

20 1.10 

1.12 

1.15 

1.16 

1.15 

1 . 14 

1.14 

1.15 

1.15 

1.13 

1 .11 

1.09 

30 1.09 

1.11 

1.14 

1.15 

1.15 

1.14 

1.15 

1.15 

1 . 14 

1.11 

1.09 

1.08 

40 1.07 

1.09 

1.12 

1.14 

1.15 

1.14 

1.15 

1.15 

1.13 

1.10 

1.08 

1.07 

50 1.06 

1.08 

1.11 

1.13 

1.14 

1.14 

1.14 

1.14 

1.11 

1.09 

1.07 

1.05 

60 1.05 

1.06 

1.09 

1.11 

1.14 

1.14 

1.14 

1.14 

1 . 11 

1.09 

1.07 

1.05 

70 

1.05 

1.07 

1.10 

1.13 

1.15 

1.14 

1.11 

1.09 

1.05 

1.04 

— 

80 

— 

1.05 

1.09 

1.14 

1.15 

1.14 

1.11 

1.07 

1.04 

— 

— 

90 

_ _ _ _ 

_ 

1.09 

1.14 

1.16 

1.15 

1.11 

1.06 








The LI-COR LI-200SB pyranometer sensor does not have an 
ideal spectral responsivity curve over the spectral irradiance 
range of blue sky and cloud coverrcf the following were known: 

1) Spectral irradiance of the sky ceor calibration; 

2) Spectral irradiance of the sky conditions at the time of 
data collection; 

3) LI-200SB pyranometer relative spectral responsivity curve. 

In reality, this is not practical because of the difficulty 

and expense involved in obtaining spectral correction factor 
experimentally. This can be derived from these measurements: 

G: Global solar radiation (total sun plus sky radiation on 

a horizontal surface) using the LI-200SB pyranometer. 

D1 : Diffuse solar radiation (sky radiation) using the LI- 

200SB and 2010S Miniature Shadow band (corrected for band 
obstruction) . 

De : Diffuse solar radiation using an Eppley Precision 

Spectral Pyranometer (PSP) and an Eppley shadow band 
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(corrected for band obstruction). 


On Figure 9.5, values of Dl/G (X) are plotted on the x-axis 
and the corresponding values of Dl/De (X) are plottedon the y- 
axis. The following equation is used to make the spectral 
correction of the LI-200SB pyranometer. All measurements were 
made in W m“2, although other units can be used since the 
correction factor is dimensionless. 

Dc * D1 

1.17 - 1 

1.2 «■ 11.8 <x) 

where Dc is the corrected diffuse radiation and x = Dl/G. The 
curve represented on Figure 9.4 is a plot of the denominator in 
the above equation. 

This equation applies only to solar radiation measured out 
doors and not greenhouse, growth chamber, artificial lighting 
conditions or under a plant or tree canopy. 

Example: Calculate the corrected diffuse solar radiation at 

a latitude of 30° N during March, where G = 800 W m-2 and the 
diffuse component measured by the L1-200SB and 20105S = 60 W m“2 
(uncorrected for band obstruction) . 

1. Correction for band obstruction (Table 9.1): 1.14(60) = 

68.4 W ra-2. 

2. Spectral correction: 

D c = 68.4 = 95.4 W m"2 

1.17 - 1 

1.2 + 11.8(68.4/800) 

Corrected diffuse radiation (D c > - 95 W m-2. 

A plot of the corrected diffuse radiation determined by the 
LI-C0R pyranometer versus diffuse radiation determined with the 
Eppley PSP pyranometer is given in Figure 9.5. 

IMPORTANT: When using the LI-200SB Quantum Sensor or LI- 

200SB Photometric Sensor, only the band obstruction correction is 
needed since these ensor3 have spectral responsivity curves that 
match very closely their respective ideal response curves. 


m 

: 



LI-200S VS EPPLEY DIFFUSE COMPARISON 



A31dd3/S002IH 


Figure 9.4. Spectral correction for the LI-200SB pyranometer. 
Values of Dl/G (JO are plotted on the x-axis and the 
corresponding values of Dl/De (JO are plotted on the y-axis. 
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Figure 9.7. Overlay of offset voltage modifications to ADC 1. 
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Figure 9.9. ADC-1 line carrier <BSR> control logic. 
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9.6 Surface Energy Balance Station Operating Instructions 
9.6.1 Sample Screen Display with Channel ID's 


1 

CHAN ENG 

CHAN 

ENG 

CHAN 

ENG 

1 

1 

1 

G 

2 

Q» 

3 

Kdn 

1 

1 

4 

Kup 

5 

D 

6 

THRt 

1 

1 

7 

Home 

6 

Udir 

9 

Tar 

1 

1 

10 

Twr 

11 

Tal 

12 

Twl 

1 

1 

13 

, Tsoil 

14 

Tthr 

15 

U 

1 

1 

16 

THRb 





1 



Screen 

Display 

* ' 

legend 



Chan Abbr 


Description 

• 



/ 

Units 


i 

G 

= 

soil heat flux 

W m-2) 

2 

Q* 

= 

net short and long wave radiation 

(W m-2) 

3 

Kdn 

= 

downward shortwave radiation 

(W m-2) 

4 

Kup 

= 

upward shortwave radiation 

(W m-2) 

5 

D 

= 

downward diffuse (sky) radiation 

(W m-2) 

6 

THRt 

— 

signal from top of total 
hemispherical radiometer 

(W m-2) 

7 

Home 

= 

AEM home signal 

(mv) 

8 

Udir 

= 

wind direction 

(degrees) 

9 

Tar 

= 

air temperature, right side 

(Celsius) 

io 

Twr 

= 

wet bulb temperature, right side 

(Celsius) 

11 

Tal 

= 

air temperature, left side 

( Celsius) 

12 

Twl 

= 

wet bulb temperature, left side 

(Celsius) 

13 

Tsoil 

= 

soil temperature, average top 5 cm 

(Celsius) 

14 

Tthr 

= 

temperature of THR 

(Celsius) 

15 

U 

= 

wind speed 

(m s-1) 

16 

THRb 


signal from bottom of total 
hemispherical radiometer 

(W m-2) 


9.6.2 Energy Balance Station Maintenance Checklist 

Daily maintenance and inspection is necessary in order to 
assure proper operation of the energy balance stations. The 
following checklist is designed to aid in systematic operation. 

1 . Observe system operation 

Check computer display 

Check time (set to MST) 

2. Water psychrometers 

Note water level in psycrometer reservoir (e.g. 

1/2 full, 2/3 full, etc.) 

RECEDING PAGE BLANK NOT FILMED 
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Left hand psychrometer 

Right hand psychrometer 

Refill bottles 

Before replacing bottle, squeeze it in 

inverted position until water no longer 
runs out after squeeze pressure is 
removed . 

Check feed tube position with gauge 

Insert feed bottle into psychrometer. If 

properly seated, bottom of bottle will be 
below lip of holder 

3. Battery check 

V Main storage battery. Replace if less 

than 11.5 V. 

< V Offset battery 

< V Recorder battery) 

4. Automatic exchange system (AEM) operation 

Observe AEM operation for at least two cycles 

Check belt tension 

NOTE: If for any reason the AEM is not operational, 

note the separation distance and height of the 
psychrometer pair for use in later analysis. 

5. Radiation sensors 

Check shadow band; is Licor sensor in shade? 

Check dessicant in net (Q> and total hemispheri- 
cal <THR) radiometers 

Clean radiometer domes if they are dirty 

6. Record aky/cloudineas conditions 

7. Program changes; should only be entered 

immediately after the system has written to tape 
to minimize data loss. 

9.6.3 Maintenance Equipment Checklist 


1. Distilled water 

2. Blank cassette tape 

3. Voltmeter 

4. Psychrometer gauge 

5. Paper towels/kleenex 

6. Black tape 
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10. Emergency storm shelter kit 

11. Rain gear 

12. Record book 



i 

C13. RAM cartridge) 

<14. NEC battery pack) 

15. _ Tools 



level 
compass 
scotch locks 

straight blade screwdriver 
Phillips screwdriver 
wire cutters 
pliers 

3/8 - 7/16 and 1/2 - 9/16 wrenches 
sandpaper strips 
rubber cement 



9.6.4 Operating the data acquisition program SAMP. BA 



9. 6. 4.1 Move the cursor to SAMP. BA on the menu and hit the 
"RETURN" key to begin program execution. Basic program operation 
from this point on is automatic and requires no operator inter- 
vention . 

9. 6. 4. 2 Six single keystroke commands have been implemented to 
control various system functions without interruption to data 
collection. These are activated by simply typing the single key 
as defined below. Five of the six commands are implemented in 
the current version of SAMP. BA. 

"R" : Shifts the display to raw data units. The display 

remains in this mode until changed by the operator. Data is from 
last instantaneous sample. 

**E" : Shifts the display to engineering units. The display 

remains in this mode until changed by the operator. 

Data is from last instantaneous sample. This is the 
default display. 


I 

fJ 

I 

I 

I 

I 


I 
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"C" : Shifts the display to show the most recently computed 

energy and radiation balance. The display reverts to 
the previously selected option after the next sample is 
completed . 

"G" : Shifts the display to show the air and wetbulb 

temperatures and gradients of the 30-second sample, 
present 6-minute average, past 6-minute average and the 
12-minute average. 

"P" : Directs program output to the parallel printer port as 

well as to the NEC display screen. WARNING: IF THIS 

OPTION IS SELECTED AND A PRINTER IS EITHER NOT IN- 
STALLED OR NOT ON LINE, THE PROGRAM WILL STOP AND DATA 
COLLECTION WILL BE TERMINATED. 

"0“ : Cancels the "P" command, directing program output to 

NEC display only. 

9.6.5 Changing cassette tape 

Data is recorded on standard audio tape cassettes. Tapes of 
longer playing time than C-90 are not recommended. 

a. Wait until the system is not writing to the tape, and 
ensure that adequate time is available before the system is 
scheduled to write data to the tape. This currently occurs on 
every hour and half hour. 

b. Depress the STOP/EJECT button to release the SAVE/LOAD 
function keys, and then a second time to eject the old tape 
cassette . 

c. Write the tape counter reading on the sticker on the 
upper left corner of the cassette, alonq with the date and time 
the cassette was removed. Ensure that the system number or site 
name <1 = PNL, 2 = SKIN, 3 = ESIDE, 4 = WSIDE and 5 = WPL) is 
written on the cassette. 

d. After writing the system number or site name, and cur- 
rent date and time on the new cassette, insert it into the re- 
corder and shut the cover. 

e. Rewind the new cassette (the full reel should be to the 
left as you face the recorder) . 

f. Reset the tape counter by pushing the black button to 
the right of the counter. 

g. Advance the tape past the leader using the FF key. A 
reading of 2 on the tape counter is sufficient. 

h. Simultaneously depress the SAVE and LOAD key, and the 
recorder is ready to accept data. 
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i. It is recommented that time be allowed for the operator 
to observe correct operation of the recorder. The signal from 
the computer can be heard by use of the MONITOR button in the 
lower left hand corner of the cassette. 

9.6.6 Main battery maintenance 

Battery life. Average current consumption of a station can 
be estimated by summing the current consumption for each compon- 
ent, taking into account the variable duty cycle for the psychro- 
meter blower and tape recorder. The resulting average current 
drain is about 233 ma (Table 9.2>. A 12 V battery with an 125 
Ampere-hour capacity could be expected to last a maximum of 536 
hours <125 Amp-hours/ .585 Amps), or about 23.4 days. It is 
recommended that the battery be replaced when about 25x of this 
capacity remains, or in the case of the above example, after 16.8 
days. x 

Table 9.2. Power consumption of the data acquisition system 
from a 12 Vdc power source. Efficiences of all regulators are 
included. -y>'‘ . 


Component 

Current 

<ma) 

Duty 

Cycle 

(X) 

Average 

Current 

<ma) 

Power 

<mW) 

ADC-1 

8 

100 

8 

96 

Auxiliary module 

7 

100 

7 

84 

Computer <32k RAM) 

75 

100 

75 

900 

Cassette recorder 

65 

5 

3 

36 

Psychrometer blowers 

100 

100 

lOO 

1200 

AEM 

500 

6 

40 

480 

Totals 

705 


233 

2796 


BATTERY REPLACEMENT PROCEDURE 

a. Remove the battery covers of both the old and replace- 
ment battery, and place the batteries side by aide. 

b. One at a time, remove the insulating cover from the 
spade lug on extra set of battery cables. Connect the black 
cable to the negative <-) terminal of the new battery, and then 
the red cable to the positive <■*■) terminal. 

c. One at a time, remove the battery cables from the old 
battery, covering the terminating spade lugs with the insulating 
covers . 

d. Remove the old battery, replace the covers. 

NOTE: Minimize the time that the two batteries are 
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connected in parallel, since connecting them in such a manner 
will tend to discharge the new battery. 

9.7 Reading Tape Cassettes 

Raw data is saved on cassette tape. These tapes can be read 
and transmitted via NEC computer's serial port to another compu- 
ter for processing using the program READT2.BA and the following 
instructions. These instructions assume that the receiving com- 
puter is a Sierra Data Sciences model SBC-100 computer <SDS) 
operation under the CP/M operating system. 

On NEC: 

1. Connect NEC computer in placie of the terminal to the SDS. 

2. Put the NEC into terminal mode. 

3. Move the cursor to "TELCOM" on the menu and hit "RETURN". 

4. Set the serial port protocol to "8N82XN" using the "STAT" 
command (f.4 key). 

5. Select TERM mode <f5 key). NEC is now acting as a terminal. 
On SDS (using NEC as terminal): 

6. Set the SDS computer to use XON/XOFF protocol using the 
GENMOD program (item 32 in the menu). 

7. Insert proper floppy to store data on. 

8. Type control C 

9. Type D: 

10. Type A: ED fname (invokes CP/M editor on file fname) 

11. Type I (insert mode) 

12. Type shift f5 twice to return control to NEC menu 

On NEC: 

13. Move cursor to READT2.BA; type RETURN. 

(data is now transmitted). 

(When the transmitted data fills the data buffer of the 
receiving computer, or the date changes, data transmission will 
cease and the program will stop and issue a buffer full message. 
If so, return the NEC to the terminal mode in TELCOM, and flush 
the data buffer or close the file and open a new one to store the 
data using the following procedure.) 

On NEC: 

14. Type CONTROL C (STOP key) when "Buffer full message" occurs. 

15. Return to menu (shift 5) 

16. Run TELCOM and TERM: by repeating steps 3 and 5 

On SDS: (To flush buffer and continue data transfer) 
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22. Note time and dat,e - of laa data record on display. 

23. Type SHIFT and STOP simultaneously. (Response: ?I0 ERROR 

i OK 

24. Return to menu (type shift £5). 

25. Run TELCOM and TERM: by repeating steps 3 and 5 / 

26. Type CONTROL Z then RETURN (exits insert mode) 

27. Type E, then RETURN (saves file on disk D:) 

28. Type A: PIP C : =D : f i lename (backs up file on the^C: disk.) 

29. Repeat process for another tape. 

To check for successful operation, enter "TYPE fname" and 
RETURN to list the data file. Editing may now be accomplished 
using ED.COM or Wordstar in the (N)ondocument mode. 

9.8 Program Listings 

The following programs are listed in this Appendix: 



SAMP. BAS (with cross reference listing) 

READT2.BA 
ADCTST . BA 
SAMPE.BAS 

SAMPP.BAS (with cross reference listing) 

SUMMARYE . BAS 
PL0T4.BAS 
PL0T5.BAS 

Cross reference listings, included for some programs, 
consist of two parts. The first gives line numbers referenced in 
program statements such as GOSUB's and GOTO's (first column), and 
the line number(s) in which those references occur in succeeding 
columns. The second gives an alcal list of all variables used in 
the program, followed by a list of line number (s) in which those 
references occur. 

The pages of each listing are numbered consecutively 
starting at page 1, with the name of each program or cross 
reference listing printed at the top of the page with the page 
number. Cross reference listings are distinguished by the 


9.8.1 

9.8.2 

9.8.3 

9.8.4 

9.8.5 

9.8.6 

9.8.7 

9.8.8 
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suffix "CRF" instead of the "BA" or "BAS" used for the program 
lisings . 
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Appendix 9.8.1 SAMPX.BA, A sampling program for the NEC computer. 

10 ' PROGRAM SAMPX.BA FOR NEC 8201 AND ADC-1 6/8/84 

2100 

IS ' LAST MODIFIED <SAMP2> 101211/19/84 

20 ' SYSTEM 1 LOWER VALLEY <PNL) 

25 ' USES CUSTOM INTP. SUBR. FOR COUNTERS 
40 ' 

45 MAXFILES=3 : CLEAR 100,61999! 

50 CLS: SCREEN 0,0 
55 DEFINT I-N 
60 GOSUB 115 : GOSUB 116 

98 ' 

99 GOTO 1000 * • ' * * * TO START OF PROGRAM » * * 


lOO 

P 





105 

P 




JUMP TABLE 

no 

P 





115 


GOTO 

9000 

P 

INITIALIZE 

116 


GOTO 

2100 

P 

MISC. CONSTANTS 

120 


GOTO 

500 

P 

CLOCK 

125 


GOTO 

300 

P 

SCREEN OUTPUT 

130 


GOTO 

400 

P 

E/RAD BALANCE OUTPUT 

135 


GOTO 

600 

P 

TC; MV TO C 

140 


GOTO 

700 

P 

SAMPLE A/D <3000) 

145 


GOTO 

800 

P 

GET DIGITAL INPUTS 

150 


GOTO 

900 

P 

A/D TO MV CONVERSION 

155 


GOTO 

1300 

P 

VECTOR WIND DIR. 

160 


GOTO 

1500 

P 

RTD, MV TO C 

165 


GOTO 

1600 

P 

DIGITAL OUTPUT 

170 


GOTO 

1700 

P 

BREB 

175 


GOTO 

1900 

p 

RADIATION BALANCE 

180 


GOTO 

9100 

P 

CONTROL PARAMETERS 

185 


GOTO 

9200 

0 

CALIBRATION FACTORS 

190 


GOTO 

2000 

P 

MISC. FUNCTIONS 

195 


GOTO 

950 

P 

A/D UNITS TO MV 

200 


GOTO 

12000 

P 

ERROR TRAPPING 

205 


GOTO 

1200 

P 

AUDIO CASSETTE DUMP 

210 


GOTO 

2500 

P 

PRINTER OUTPUT 

215 


GOTO 

1400 

P 

DIGITAL CASS. DUMP 

290 

P 

- - 

- - - - 

- 


300 

P 





301 

P 

DISPLAY SUMMARY 

302 

P 






305 CLS: LOCATE 0,1 

310 IF DSP THEN MCOL=4 ELSE MC0L=3 

315 FOR 1=1 TO MCOL 

316 IF DSP THEN PRINT "CHAN RAW ELSE PRINT"CHAN ENG 

317 NEXT: PRINT 

320 FOR 1=1 TO M7 STEP MCOL 

322 FOR K=I TO MCOL+I-1 

323 IF K>M7 THEN 335 

325 IF DSP THEN 326 ELSE 328 
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326 

328 

330 

335 

340 

350 

400 

401 

402 
405 
407 
410 
415 
420 
430 
435 
450 
455 
460 
465 
470 
475 
490 

499 

500 

501 

502 
505 
510 
515 
520 


PRINT USING "## ##### " ; C < K > ;A(K) ; : GOTO 330 

PRINT USING "### #####. ###";C<K);INT<1000*A1<K>>/1000 

next:print 

NEXT 

IF DSP THEN LOCATE 0,7: PRINT CC-OC> /N2; C-OC; C; 

RETURN 


DISPLAY ENERGY, RAD BALANCE 


CLS:GOSUB 120 
LOCATE 0,2 

PRINT " H E B Q G'* 

PRINT USING F1S;H;E;B;Q;G 

PRINT " KDN KUP LUP LDN THR“ 

PRINT USING F2S; KDN ; KUP ; LUP; LDN ; THR 

PRINT ** TT TWT E DT DE" 

PRINT USING F3S;T9;W9;E9;DT;DE; 

LOCATE 0,0:PRINT " U DIR" 

PRINT USING F1$;A1(17) ;A1<8> ; 

RETURN 

PRINT “ CP L S GAMMA RHO HOME" 

PRINT USING F4S;CP;XL/1E«-06;S0»1000;G1*1000;R0;P1; 
RETURN 

' CLOCK - HOURS/MINUTES/SECONDS (HR/MIN/SEC) 


SEC=VAL (MIDS (TIMES, 7 ,2) ) 

IF S8=59 AND SEC=0 THEN 570 
IF S80SEC THEN 515 ELSE RETURN 
LOCATE 23,0 
PRINT TIMES; " DATES; 


522 DSPS=INKEYS: IF DSPS="" THEN 545 
525 IF DSPS=CHRS<27> THEN E2=l ' ESC? 

530 IF DSPS="T" THEN TPE=1 

533 IF DSPS="R" THEN DSP=1:G0SUB 125 

534 IF DSP$="E" THEN DSP=0:G0SUB 125 

535 IF DSPS="C" THEN GOSUB 130 

540 IF DSPS="P" THEN PRT=1 

542 IF DSP$="0" THEN PRT=0 

545 S8=SEC 

550 HR=VAL (MIDS (TIMES, 1 , 2> ) 

555 YR=VAL<MIDS (DATES, 1,2)) : M0= VAL < MIDS < DATES, 4 , 2> ) : 

DA=VAL <MIDS <DATES ,7,2) ) 

570 MIN = VAL < MIDS (TIMES ,4,2) ) : RETURN 
600 ' 

605 ' THIS SUBROUTINE CONVERTS READINGS FROM A THERMOCOUPLE 
AND 

610 ' REFERENCE JUNCTION IN A/D UNITS TO DEG. C. 

615 ' 

625 V=A1 < IC) «■ AKREF) 

627 A1<IC)=B1 * V*B2*V~2*B3*V~3+B4* V~4 

640 RETURN 


*** SAMPLE A/D (ADC-15; CONVERT TO DECIMAL *** 


700 ' 

704 ' 

705 ' 

710 FOR K2=l TO NANLG :CN=C1 (K25 
715 FOR Jl = l TO N5 

725 POKE IOB, CN : POKE ITOB,l: ' SELECT CHANNEL; 

START A/D 

730 X=PEEK(IB> GARBAGE CHARACTER 

735 POKE IOB, 161 : POKE IT0B,1 ' GET ADC-1 HIGH 

BYTE/STATUS 

745 HBYTE=PEEK(IB5 SAVE HIGH BYTE 

750 IF (HBYTE AND 128) <> O THEN 735 ' CHECK STATUS FOR 

A/D FINISHED , , 

755 POKE I0B,145:P0KE IT0B,1: ' GET ADC-1 LOW BYTE 

765 LBYTE=PEEK(IB> : ' SAVE LOW BYTE '; 

770 HMASK=HBYTE AND 15 ' MASK 4 HIGH ORDER 

BITS FROM A/D , 

775 Y=LBYTE+256*HMASK ' ' COMBINE ALL 12 BITS 

FROM A/D 

780 IF (HBYTE AND 165=0 THEN Y=-Y ' FIX SIGN^tF 

NEGATIVE FLAG SET 

782 IF J1=N5 THEN 785 ELSE 787 

785 A(K2)=Y 

787 NEXT 

790 NEXT: RETURN 

800 ' 

801 ' *** SAMPLE & RESET COUNTERS *** 

802 ' 

805 OC=C 

810 Cl =PEEK (115 :C2=PEEK (125 :C3=PEEK(I35 
820 C=C1 ■*-256*02+65536 ! »C3 
822 IF C9= 1 THEN RETURN 
825 A (M75 =C-OC 

827 IF A (M75 <-lE+06 THEN A (M75 = A (M7 5 +1 .6843E+07 
830 A1 (M75 = A(M75 /N2 
880 RETURN 

900 ' 

901 ' A/D UNITS TO MV 
905 ' 

910 FOR 1=1 TO NANLG 
915 A1(I)=A(I5 *G2 (15 +B1 (15 

925 NEXT 
940 RETURN 

949 ' 

950 ' MV TO ENG. UNITS, LINEAR 

951 ' 

960 A1(IC5=A1(IC5*G(IC5+B(IC5 

970 RETURN 
1000 ' 

1005 ' MAIN SAMPLING LOOP 
1010 ' 

1015 LOCATE 0,0:PRINT "WAIT FOR SECONDS = O " 




1020 GOSUB 120: IF SEC>2 THEN 1020 
1025 GOSUB 120: IF SEC=0 THEN 1025 

1030 LOCATE 0,0:PRINT “SAMPLING INITIATED “ : ' ON ERROR GOTO 
12000 ' ???? 

1032 J9=0:P0KE I1,0:P0KE I2,0:P0KE 13,0 

1035 J9=J9*1:N6=0:H9=0: LOCATE 18,0:PRINT J9 

1040 FOR Kl=l TO M : D<K1)=0 : A2(K1)=0 : NEXT 

1045 A1=0:A2=0 ' ZERO VECTOR COMPONENTS OF WIND DIRECTION 

1050 GOSUB 120: IF TPE THEN GOSUB 215 

1052 IF INT(SEC/N2) OSEC/N2 THEN 1050 ' UPDATE CLOCK TILL 

TIME TO SAMPLE 

1055 GOSUB 145 ' SAMPLE COUNTERS 

1060 GOSUB 140 ' SAMPLE A/D' S 

1065 IF INT < <MIN+N1-1-N8)/N1) <>INT < <MIN + N1-1)/N1) THEN LOCATE 

0,0:PRINT "WAIT FOR EQUILI. " : GOTO 1050 ' SKIP 1ST N8 POINTS 

1068 GOTO 1085 ' SKIP HOME CHECK 

1070 IF ABSCA(HOME)) >400 THEN 1085 

1075 IF INT <MIN/N1 )=MIN/N1 THEN 1085 

1080 H9=H9+1 : IF H9^=2 THEN 1050 

1085 N6=N6+1 

1087 GOSUB 150 ' A/D UNITS TO MV 

1090 FOR IC = 1 TO M7 

1093 D<IC)=D<IC) ♦A(IC) ' SUM RAW DATA 

1095 ON N < IC) GOSUB 135,160,155,195,195 

HOO A2 < IC) =A2 <IC)+A1<IC) ' SUM ENG UNITS 

1110 NEXT 

1120 GOSUB 125 ' UPDATE DISPLAY 

1125 LOCATE 0,0: PRINT "SAMPLE BELOW SAVED" ; J9 ;: LOCATE 

35, 7: PRINT N6 ; 

1130 GOSUB 120 ' GET TIME 

1135 IF E2=l THEN 1145 ' EXIT IF "ESCAPE" LAST KEY PRESSED 

1140 IF INT <MIN/N1 ) =MIN/N1 AND SEC+N2>59 THEN 1142 ELSE 1050 

1142 C9=l :C5=C4: GOSUB 145:C9=0 ' COUNTER 

1143 C4=C : DC=C4-C5 : IF DC<-lE+06 THEN DC=DC+ 1 . 6843E+07 

1144 DSX<J9,M7)=DC/N1 :A2(M7) =DC*G < M7) / <N1 *60) ♦B(M7) 

1145 IF N6<10 THEN I=N6 ELSE 1=0 

1146 DSX<J9,M-1) =HR*1000! >MIN*10«-I ' HRS/MIN 

1150 DSX< J9,M) = MO *100 «• DA 

1155 GOSUB 165 ' REVERSE BOWEN RATIO DEVICE 

1160 FOR 1=1 TO NANLG 
1165 A2<I)=A2<I)/N6 

1170 IF N < I ) <> 3 THEN DSX < J9 , 1 ) =D < I > /N6 ELSE DSX ( J9, I ) =D < I ) 

1175 NEXT 

1180 J=NCRTD : K9= 1 : IF TPE <>1 THEN GOSUB 170 ' BREB 

1182 IF PRT=1 THEN GOSUB 210 ' PRINTER OUTPUT 

1185 IF E2=l THEN STOP 

1190 IF MIN MOD N4*N1=0 THEN GOSUB 205 

1191 IF J9>=N3 THEN J9=0 

1192 GOTO 1035 
1195 ' 

1200 ' SAVE RAW DATA ON CASSETTE TAPE 
1205 ' 
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1210^IS=J9-N4-*-l:IE = J9:TPE = 0:IF IS<1 THEN IS=1 

1220 OPEN "CAS: DATA" FOR OUTPUT AS #2 

1225 CLS: LOCATE 0,0 

1230 PRINT "WRITING TO TAPE 

1235 FOR J1 = IS TO IE 

1240 FOR 1=1 TO M 

1250 PRINT #2, DS*e< J1 , I ) ; 

1260 NEXT 

1265 LOCATE 23,0:PRINT TIMES;" REC= ";J1; 

1270 NEXT: CLOSE #2 : RETURN 
1300 ' 

1305 ' VECTOR AVG WIND DIRECTION 
1315 ' 

1330 A7= < A1 < IC > *G < IC) +B < IC) ) /DPR 
1340 A1=A1 + C0S< A7) : A2=A2-*SIN<A7) 

1345 IF AlOO THEN A3=ATN < A2/A1 ) ELSE A3=SGN < A2) *PI/2 
1350 IF SGNCA1X0 THEN A3 = A3«-PI 

1360 IF SGN < A1 ) >0 AND SGN<A2)<0 THEN A3=*A3 + 2«PI 
1380 D(IC)=A3*DPR:A1<IC)=D<IC> 

1390 RETURN 

1400 ' DIG CASSETTE DUMP 

1410 IF I0>0 THEN 1440 

1420 POKE IP, 201 ' TURN OFF COUNTER 

1425 OPEN "COM :6N82NN" FOR OUTPUT AS #3 

1430 PRINT"DISC. ADC-1; CONN. DIG REC." 

1432 INPUT"RECORDER ON; PRESS RETURN" 

1435 io=i:js=o 

1440 J8= J8+ 1 

1445 LOCATE 0,0:PRINT “DIG 0/P IN PROG. ";J8; 

1450 FOR 1=1 TO M 

1455 PRINT #3,DSSe( J8,I) ; 

1460 NEXT 

1475 IF J8>=N3 THEN 1480 ELSE 1440 

1480 10=0 :TPE=OJ CLOSE #3 

1481 PRINT"DISC. DIG REC.; CONN. ADC-1" 

1482 INPUT” RECORDER OFF; PRESS RETURN" 

1500 ' 

1505 ' CONVERT RTD READINGS TO DEG. C. 

1510 ' 

1545 T=(Al(IC))/RCa)/B2(IC) 

1550 A1 ( IC) =-245 . 665 +T* (235. 476*10 . 189 *T ) 

1565 RETURN 
1600 ' 

1601 ' PULSE BOWEN RATIO DEVICE 

1602 ' 

1605 ' CHANNEL: 123456 ALL OFF 

1607 ' 

1610 POKE 10B , 65 : POKE ITOB, 1 : X=PEEK < IB) 

1615 FOR 1=1 TO 500 : NEXT 

1620 POKE I0B,64:P0KE ITOB, 1 : X=PEEK < IB) 

1625 RETURN 
1700 ' 
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1705 ' ONLINE CALCULATIONS 

1715 ' SUBS, 6 = PRESENT VAL., SUB7,8 = PAST VAL. , SUB9,0 = 
RUNNING AVE. 

1720 ' 

1725 Q5=A2<2) :G5=A2<1) 

1730 TAV5=<A2(J)+A2<J+2) )/2:WAV5= <A2< J + l) +A2< J+3 ) )/2 
1735 P1=SGN (A (HOME) ) : IF Pl = -1 THEN IALT=0 ELSE IALT=2 
1740 T5=A2< J + IALT) : T6 = A2 < J-IALT + 2) : W5 = A2 < J + I ALT + 1 ) :W6=A2<J- 
IALT+3) 

1744 ' 

1745 ' *** FIND RUNNING AVERAGES **» 

1746 ' 

1750 Q= (Q7 <K9) +Q5) /2:G= <G7 <K9> +G5) /2 

1755 T= (TAV5+TAV7 <K9) ) /2 : TW= <WAV5+WAV7 <K9> ) /2 

1760 T9= CT5+T7 <K9) ) /2 : T0= (T6 + T8 <K9 ) > /2 : W9= < W5 + W7 <K9> ) /2: 

W0= <W6 + W8<K9> ) /2 

1764 ' 

1765 ' *** SAVE PRESENT VALUES *** 

1766 ' 

1770 G7 <K9) =G5 :Q7 <K9) =Q5 

1775 TAV7 < K9 > =TAV5 : WAV7 C K9) =WAV5 

1780 T7 ( K9 ) =T5 : T8 < K9 ) =T6 : W7 < K9 ) = W5 : W8 ( K9 ) = W6 

1784 ' 

1785 ' *** MISCELLANEOUS PARAMETERS *** 

1786 ' 

1790 TT=T:W=TW:GOSUB 2015:EA=EFN 

1792 CP= <239.9«-440.9* .622*EA/ <P-EA> >/ .2388 

1795 XL = 2 . 5013E + 06-2366**TW : GOSUB 2030:S0 = S 

1800 Gl=P»CP/< . 622*XL) : R0=3 . 4838* (P- . 378*EA ) / (T+273 . 16) 

1805 S1=9.81/CP:TT=TW:G0SUB 2030 

1810 S2=9 . 81 * < l/CP+3 . 4857E~03*EA/ < 273 . 16+T) /G1 ) / < 1 +S/G1 > 
1815 S7=9 .81*3. 4857E-03*EA/ < 273 . 16 + T) 

1819 ' 

1820 ' *«* GRADIENTS *»« 

1821 ' 

1825 TT=T9 : W=W9 : GOSUB 2015:E9=EFN 

1826 TT=T0:W=W0:G0SUB 2015:E0=EFN 
1830 DT=T9-TO+Sl*DELZ<K9> 

1835 DE=E9-E0+S7*DELZ(K9) 

1859 ' 

1860 ' *** BOWEN RATIO USING T, E *** 

1861 ' 

1865 B=G1*DT/DE 

1870 H = < -Q-G) /(l+l/B) :E = H/B 

1872 ' 

1873 ' «»» RADIATION BALANCE *** 

1874 ' 

1877 KUP=-A2<4) :KDN=A2<3) 

1880 IF KDN<=0 THEN A=0 ELSE A=-KUP/KDN ' ALBEDO 
1890 THR = SIGMA« < A2 ( 14 > -*-273 .16) ''4 + A2 <6) 

1895 LUP=-THR-KUP:LDN=Q5-KUP-KDN-LUP 
1910 RETURN 
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1950^' 

2000 ' * * * MISCELLANEOUS FUNCTIONS * * * 

2005 ' 

2015 

ESAT = <E < 1) + W«* (E <2) + W» <E<33 + W« <E <43 +W» <E <53 +W* <E<63+W*<E<7333333)) 
/10 

2020 EFN=ESAT-6.6E-04*<1+1.15E-03*W3*P*<TT-W3 

2025 RETURN 

2030 

S= <S < 1 3 +TT* <S <23 +TT* <S <33 +TT * (S <4 5 +TT» <S <53 + TT«* <S<63 + TT* <S<733333 

333/10 

2035 RETURN 

2050 ' 

2100 ' * * * MISCELLANEOUS CONSTANTS * * * 

2105 ' 

2115 E < 1 3 =6 . 1078 . 

2116 E<23= .44365185# 

2117 E <33 =.014289456# / 

2118 E<43 =2.65064850-04 

2120 E (5) =3 . 031240400000003D-06 

2121 E <6> =2 . 0340809D-08 “>>>•'. 

2125 E <73 =6 . 136820900000027D-1 1 

2126 ' 

2130 S <13 =.44381 

2131 S (2 3 =.028570026# 

2132 S < 3 3 =7 . 93805E-04 

2133 S <4 3 =1 . 2152151D-05 

2135 S <53 =1 . 0365614D-07 

2136 S <6 3 =3 . 532421800000003D-10 

2140 S < 7 3 = -7 . 090244800000048D-13 

2141 ' 

2145 Bl=25. 661297# 

2146 B2=- .619548690000003# 

2147 B3= .022181644# 

2148 B4= -3 . 5509E-04 
2150 RETURN 

2500 ' 

2505 ' PRINT SUMMARY 
2510 ' 

2515 LPRINT TIMES;” " ; DATES 

2520 MC0L=4 

2525 FOR 1=1 TO MCOL 

2530 LPRINT "CHAN RAW ENG 

2535 NEXT: LPRINT 

2540 FOR 1=1 TO M7 STEP MCOL 

2545 FOR K= I TO MCOL+I-1 

2550 IF K>M7 THEN 2575 

2565 LPRINT USING ”### ##### 

#####. ###";C<K3 ;DSX<J9,K3 ; I NT < 1000« A2 <K3 3 /1000; 

2570 NEXT: LPRINT 

2575 NEXT: LPRINT 
2580 ' 


ENG 
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DISPLAY ENERGY, RAD BALANCE 


■ 

■ 

i 

a 

a 

a 

a 

a 

a 


a 

a 

■ 


2600 ' 

2601 ' 

2605 LPRINT " H E B Q G"; 

2610 LPRINT ** KDN KUP LUP LDN THR" 

2615 LPRINT USING FIS ; H ; E ; B ; Q ; G ; 

2620 LPRINT USING F2S ; KDN ; KUP; LUP; LDN ; THR 

2625 LPRINT ” TT TWT E DT DE“ ; 

2630 LPRINT ” U DIR" 

2635 LPRINT USING F3S ; T9 ; W9 ; E9 ; DT; DE ; 

2640 LPRINT USING FIS ; A2 < 17 > ; A1 <8) 

2645 FOR 1 = 1 TO lOrLPRINT **- - - " ; : NEXT : LPRINT : RETURN 
2650 LPRINT " CP L S GAMMA RHO HOME” 

2655 LPRINT USING F4S ; CP ; XL/1E+06 ; SO* 1000 ; G1 *1000 ; RO ; PI 
2660 RETURN 
9000 ' 

9005 ' * INITIALIZE CONTROL PARAMS • 

9010 ' 

9020 OPEN "COM :8N82NN" FOR INPUT AS #1 
9030 DPR=57 . 2958 ' DEGREES/RADIAN 

90S0 DELZ < 1 ) = 1 1 ' PSYCHROM. SEP. 

9055 SIGMA =5 . 6697E-08' BOLTZMAN CONST 

9070 PN=192 ' SERIAL PORT DATA 

9080 PI = 3.14159 

9095 REF=0 ' TC REF CHANNEL 

9100 HOME=7 ' HOME CHANNEL 

9105 ELEV= 1804 ' ELEVATION, M 

9106 P=101 .3- .01055*ELEV ' ASSUME STD ATMOSPHERE 
9110 OPEN "INDAT1” FOR INPUT AS #2 

9112 INPUT #2 , XS ' SKIP LABEL 

9115 INPUT #2 ,M,N1,N2,N3, N4 , N5, N8 , GO, M7 

9120 N4=N4/N1 ' SET N4=# OF RECORDS/DISK UPDATE 

9125 DIM D(M) , A <M7) ,A1 <M7) , A2<M) 

9130 DIM C<M> ,C1 CM> ,G<M7) ,B<M7) ,G2<M7) 

9132 DIM B2 (M7) , N <M7 ) , B1 <M7 ) 

9135 INPUT #2 , XS ' SKIP LABEL 

9136 INPUT #2, LG, HG, HOME, REF , 01 , 02 , RC < 1 ) , NCRTD 

9137 INPUT #2,XS ' SKIP LABEL 

9140 FOR K=1 TO M7 

9145 INPUT #2 ,C<K) , Cl <K ) ,G<K) ,B<K) ,N<K) ,XS 

9150 IF Cl < K) =0 THEN G2(K)=1/LG ELSE G2<K)=1/HG 

9155 IF Cl <K) =1 THEN B1(K)=01 

9160 IF Cl ( K ) =3 THEN B1(K)=02 

9165 Cl <K) =C1 <K) *16+C<K) -1 

9168 IF N < K) =2 THEN NRTD = NRTD + 1 

9170 IF N<K) =3 THEN NWD =K 

9175 IF N < K ) =5 THEN NDIG = NDIG-1 

9180 NEXT 

9190 FOR K=NCRTD TO NCRTD +NRTD-1 

9192 INPUT#2,B2(K) 

9193 NEXT : CL0SE#2 
9195 NANLG=M7-NDIG 

9235 F1S= ”#####.# #####.# ###.### #####.# #####.#” 
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9240 F2S= "#####.# #####.# #####.# #####.# #####.#" 

9245 > F3S= "###.### ###.### ###.### ##.#### ##.####■• 

9250 F4$ = "#### . # ##.### ###.## ###.## ###.## ##.#" 

9260 ' CALC DATA BUFFER SIZE 
9265 N3=<FRE<0) -1600) /<2»<M+1>) 

9270 DIM DSX<N3,M> 

9275 LOCATE 0,7: PRINT N3*Nl/60;" HOURS OF DATA IN BUFFER**; 

9300 ' INIT UART INTERRUPT HANDLER 

9310 IF PEEK <62000 !) =51 AND PEEK <62115 !) =201 THEN 9320 ELSE 9315 

9315 PRINT "LOAD DIG IN ROUTINE" : STOP 

9320 Il = -3413 :I2=I1+1: 13=11+2 

9325 IB=-3420:I0B=IB+l:IT0B=IB+3 

9330 IP= -3188 

9335 POKE IP , 195 : POKE IP+1 , 48 :POKE IP+2,242 
9340 OUT PN,129:' START COUNTERS 

9350 POKE I0B,64:P0KE ITOB,l:'ALL DIG. O/P'S OFF 

9500 RETURN 

12000 ON ERROR GOTO O 

12002 IF INKEYsS=CHR$<27> THEN E2=l 

12005 PRINT "ERROR ** ; ERR ; " IN STATEMENT " ; ERL 

12020 RESUME 1015 ^ 
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Appendix 9.8.5 READT2.BAS, A program for the NEC which reads 
casette tape data into the editor of the SDS computer. 

10 * READT2 READS SAMPX DATA FILE AND 
15 * WRITES IT TO SERIAL PORT 9/18/84 
20 ' LAST MODIFIED 1530 9/18/84 

25 ' 

30 MAXFILES=2 : CLEAR 100,62335! 

35 POKE -3188,201 

50 M = 19 ' FIELDS/RECORD 

60 DIM DSsc<M> ISCREEN 0,0 

100 OPEN "COM : 8N82XN" FOR OUTPUT AS #2 

120 CLS: LOCATE 0,0 

130 N =0 

140 OPEN "CAS: DATA" FOR INPUT AS #1 
200 N=N+1 

205 DATE=DSX<M) : TIME= INT <DSsc <M-1 ) /lO) 

210 PRINT N; 

220 FOR 1=1 TO M 

225 IF EOF < 1 > THEN GOTO 400 

230 INPUT #1,DSX(I> 

250 PRINT USING ”#####"; DS* ( I ) ; 

260 PRINT #2, USING "######“; DS* < I ) ; 

270 ' IF EOF ( 1 ) THEN GOTO 400 
280 NEXT: PRINT"" 

285 PRINT #2,"" 

290 IF DATEODSX <M> AND N>1 THEN CLOSE #l:FOR 1=1 TO 20:PRINT 
CHR$ < 7 ) ; : NEXT : STOP 
295 GOTO 200 
400 PRINT"" 

405 PRINT #2,"" 

410 CLOSE #1 

415 IF N>285 THEN FOR 1=1 TO 20:PRINT CHRS ( 7) ;: NEXT : INPUT "ED 
BUFFER FULL; USE #W ";XS 

420 IF DATEODSX CM) AND N>1 THEN FOR 1 = 1 TO 20:PRINT 
CHR$ < 7 ) ; : NEXT : STOP 
425 GOTO 140 



m 
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Appendix 9.8.6 ADCTST.BAS, A Test program for the ADC-1 using 
the NEC computer. 

10 ' ADCTST: TEST FOR ADC-1 6/7/84 

15 ' LAST MODIFIED 3/ 6/85 

20 ' 

25 CLS: CN = 16.* POKE -3188,201 
30 OPEN "COM : 8N82NN" FOR INPUT AS 1 

40 DIM C < 16) ,M< 16) , N< 16) , A ( 16) , OFST (16) ,0(16) ,S< 16) 

90 PN=192 SERIAL PORT DATA ADDRESS 

105 X=INP(PN) CLEAR INPUT PORT OF OLD BYTES 

107 OS=0:NO=1:N1=10:N3=1:C(1)=1:GOTO 170 
110 INPUT "GAIN /OFFSET" ; OS 

120 PRINT "A/D STABILITY AND CALIBRATION TEST" 

130 INPUT “NO. OF CHANNELS TO TEST “;N0 . 

140 INPUT "NO. OF CHANNELS TO AVERAGE' ";N1 

145 INPUT "NO. OF SCANS/SAMPLE ";N3 

150 PRINT "SPECIFIY EACH CHANNEL TO TEST " 

160 FOR K= 1 TO NO:INPUT "?";C(K) :NEXT 

165 IF C<1)=0 THEN FOR 1=1 TO 16 : C ( I ) =1 : NEXT V*’. 

170 Yl=-lE+38 : Y2=lE+38 
180 PRINT 
185 N2= lO 

190 FOR L=1 TO NO: M <L) =-10000 : N <L) =10000 : NEXT 
200 FOR L=1 TO N1 
225 GOSUB 800 

228 FOR K=1 TO NO 

229 IF L= 1 THEN OFST(K)=A(X) 

230 S(K)=S<K)+A<K) -OFST <K) :Q < K) =Q <K ) + < A < K ) -OFST (K) ) ~2 

240 IF A <K ) >M (K) THEN M(K)=A(K) 

245 IF A(KXN(K) THEN N(K)=A(K) 

250 NEXT 
260 NEXT 

270 FOR L=1 TO NO 

280 Q ( L) =SQR < ABS ( < Q < L ) -S<L) ~2/Nl ) / <N1 -1 ) ) > 

290 S(L)=S(L) /Nl+OFSTCL) 

300 NEXT 

305 PRINT "CH NO. AVE STD DEV MAX MIN" 

310 FOR L=1 TO NO 

320 PRINT USING "###";C<L>, 

330 PRINT USING "#######.## ####### . ##" ;S(L),Q(L), 

335 PRINT USING "####### #######"; M <L) , N <L) 

340 Q < L> =0 : S (L) =0 

350 NEXT 
355 PRINT 
360 GOTO 190 
800 ' 

801 ' *** SAMPLE A/D (ADC-1 >; CONVERT TO DECIMAL *** 

802 ' 

805 FOR K2=l TO NO 

810 XS=INKEY$:.IF XSO""THEN Cl=ASC(XiS) 

811 IF XSO"" THEN IF Cl>57 THEN Cl=Cl-7 
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815 IF Cl >48 THEN C<l)=Cl-48 

816 CN=C<K2) + 0S-1 
818 FOR 11=1 TO N3 
820 OUT PN, CN 
A/D 

826 Y=INP<PN) 

827 FOR K=1 TO 200: NEXT 

830 OUT PN , 128+32 

BYTE/STATUS 

835 * OUT PN , 129 
840 HBYTE=INP(PN) 

A/D 

845 IF (HBYTE AND 128) <> O THEN 830 
FINISHED 

850 OUT PN, 129+16 
855 ' OUT PN , 129 
860 LBYTE=INP (PN ) 

865 HMASK=HBYTE AND 15 
FROM A/D 

870 Y=LBYTE+256*HMASK 
FROM A/D 

875 IF (HBYTE AND 16) =0 THEN Y=-Y 
FLAG SET 

877 IF I1=N3 THEN 880 ELSE 883 
880 A(K2)=Y:PRINT USING "#####"; Y; 
883 NEXT 

885 NEXT: PRINT:' HBYTE ;HMASK;LBYTE 

890 RETURN 

lOOO Cl =VAL< INKEYS) 

1010 IF CIO 0C1 THEN CN = C1 
1015 0C1=C1 

1020 PRINT CN : GOTO lOOO 
1050 GOTO lOOO 


SELECT CHANNEL; START 
: ' GARBAGE CHARACTER 
GET ADC-1 HIGH 

: ' SAVE HIGH BYTE FROM 

:' CHECK STATUS FOR A/D 

:' GET ADC-1 LOW BYTE 

: ' SAVE LOW BYTE FROM A/D 
:' MASK 4 HIGH ORDER BITS 

:' COMBINE ALL 12 BITS 

FIX SIGN IF NEGATIVE 
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Appendix 9.8.8 SAMPEE.BAS, a program for the AT computer which 
converts the raw data from the NEC computer into engineering 
units. 

4 'note changed lines 4470-4485 removed 

5 ' SAMPE.BAS modified for AT computer and Epson LQ-1000 
printer 

8 

06/25/86 0443 

10 ' SAMPB2.BAS RANDY data analysis program 

06/13/85 1130 

15 ' Based on PROGRAM SAMPC2.BAS 

12/17/84 1143 

20 ' For Hanford Site study, Washington 

100 ' 3981 for thermal conductivity and.. 3892 for correct soil 
heat flow. ^ 

120 ' Last modified 

5/7/86 . 

140 DEFINT I-N : M=18 

150 DIM T <50) , IFLG0<30) , IFLG07 <30) „ 

155 DIM N < 25) , D <25) , F < 17) , A$ < 50) , L <50> , T$ (13) , C < 50, 4) /B$<50> 
160 DIM A2 < 20) , CH < 20) , Cl <20> , G <20) , B < 20) , G2 (20) 

165 DIM NT < 20) , B1 < 20) , FLS < 120) , NS> ( 9) 

170 DIM 

DELZ<2) ,07 (2> ,G7<2) , TAV7 <2> , WAV7 <2> ,T7<2> ,T8<2> ,W7<2) ,W8(2) 

175 DIM RC<2> ,E<9> ,S<7> 

200 ' 

lOOO GOSUB 9000 j FS»= “ " ' Microstat init 

101O GOSUB 32000 ' init clock 

2000 ' 

2005 ' MAIN PROGRAM 
2010 ' 

2030 M3=o:M2=o:Mi=o:S4=o:S5=o:s7=o 

2040 ' 

2045 ' read input file name 
2050 ' 

2055 ICOUNT = ICOUNT + 1 : NS = DXIS + FL$ < ICOUNT) : N1$=NS + ,, R" 

2056 N9$=NS-*-" . DAT" 

2065 05=4 

2070 GOSUB 6200 ' set system 

specific info 
2100 

2540 L1=M :M1=L1 ' NOUT 

2560 

2565 C=1 : D=1000 

2605 N3M=D-C«-1 

2608 

2617 N3S=DX0S+FL$ < ICOUNT ) :ZS = N3$ + FSS 

2700 

2870 FOR L= 1 TO LI 
2880 L (L) =L 

2890 NEXT 
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read data system 


2900 GOSUB 8300 ' 
parameters 

2990 PRINT F$:PRINT*’FILE: " N3$ ** IS NOW BEING OUTPUT **:J1=0 

2995 ' 

2996 ' open input and output files 

2997 ' 

3000 OPEN "I" , #1 , N93 ' OPEN **R" , #1 ,Nls3,Q5 

3050 ' FIELD #1,05 AS TS 

3060 OPEN "R‘*,#2,ZS,Q5 

3080 FIELD #2,05 AS U$ 

3100 ' 

3105 ' main computation loop 
3110 ' 

3150 FOR J = C TO D 
3155 ' 

3160 ' read data into TO 
3165 ' 

3170 FOR K=1 TO LI 

3175 INPUT #1 , T<K) ! IF E0F(1) THEN 3540 

3210 NEXT 

3427 GOSUB 3600 

3440 Jl=Jl + l:PRINT GHRS < 13) J1 INT<T <M-1 > /lO) T< M) ** **; 

3450 ' 

3455 ' write out full TO array 

3456 ' 

3460 FOR L = 1 TO LI 

3480 IF 05=4 THEN LSET US=MKS$ <T <L <L > > > ELSE LSET 

US=MKD$<T<L<L> > ) 

3510 PUT #2, L+ <L1 * < Jl-1 ) > 

3520 NEXT 

3530 NEXT: PRINT FS 

3540 PRINT "END OF FILE OUTPUT** : NS=NlsS 

3570 CLOSE #1: CLOSE #2: PRINT 

3575 GOSUB 6400 'create output 

file directory 

3580 NS=N3$:G0SUB 8300 

3590 IF ICOUNT<IFILES GOTO 2000 

3595 CHAIN "SAMPP** : END 

3600 ' 

3605 ' MAIN SAMPLING LOOP 
3610 ' 

3670 FOR 1=1 TO Ml 

3675 IF NT < I ) =3 OR I>Ml-2 THEN A2 < I ) =T ( I ) : GOTO 3689 'No action 
Time or Udir 

3680 IF NT < I ) =3 AND TCMXVDATE AND INT <T<M-1>/10) <VTIME THEN 
A2< I > =T <1 ) ♦VANE : GOTO 3689 

3685 A2(I)=T(I)»G2(I)+B1<I) ' A/D UNITS TO MV 

3689 NEXT 

3690 FOR IC=1 TO M7 

3695 ON NT < IC ) GOSUB 6700,5000,6800,6700,6700 
3700 NEXT 
3714 ' 
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Home signal 


3775^ TIME * INT ( A2 < 18) /10) 

3785 GOSUB 4400 ' 
processing 
3800 ' 

3898 FOR 1=1 TO LI : T(I)=A2(I) : NEXT 

3982 TK= . 64 + 1 . 63*CS0IL- < . 64- . 135) *EXP C - < < 17*CS0IL) A 2) ) S ' TK is 
thermal conductivit 

3983 PRINT "TK= ** , TK 

3984 T(1)=T<1)* <1-1 .92* .138* (1- CTK/.48) ) ) / <1-1 .92* . 138* <1- 
< . 94/ . 48) ) ) 

3986 'Above is heat flow correction-see Fritschen and Gay 
4300 ' 

4305 ' Checks: Tw < O, dT or dTw < .005, Tw -> T and 4095 < 
signal < -4095 
4310 ' 

4330 IFLGW7=IFLGW: IF W5<0 OR W6<0 THEN TFLGW=1 ELSE IFLGW=0 
4335 IF IFLGWOIFLGW7 THEN IPFLG=1 : PRINT USING 
**#####••; TIME; :PRINT TAB < 13) PRINT USING "##.## / 

##.##";W5,W6; .‘PRINT CHRS<13); 

4339 ' 

4340 IFLGDT7 = IFLGDT : IF ABS CT9-T0) < . 005 OR ABS(W9-W0) < .065 THEN 
IFLGDT=1 ELSE IFLGDT=0 

4345 IF IFLGDTOIFLGDT7 AND J>C THEN IPFLG=1 :PRINT USING 
**#####" ;TIME; :PRINT TAB C39) PRINT USING **##.#### ##.####**; T9- 
T0,W9-W0; SPRINT CHRS<13>; 

4346 ' 

4347 IFLGD7=IFLGD: IF ABS (T5-T6) < .02 OR ABS<W5-W6) < .02 THEN 
IFLGD=1 ELSE IFLGD=0 

4348 IF IFLGDOIFLGD7 AND J>C THEN IPFLG = 1 .'PRINT USING 
••#####••; TIME; .-PRINT TAB (57) PRINT USING **##.#### ##.####"; T5- 
T6.W5-W6; SPRINT CHRS<13>; 

4349 ' 

4355 Il=-l:FOR 1=1 TO Ml-3 S IFLG07 < I ) =IFLGO ( I ) : IF ABS <T < I > ) > =4095 
THEN IFLGO < I ) =1 ELSE IFLG0(I)=O 

4360 IF IFLGO (I) OIFLG07 < I ) AND J>C THEN IPFLG = 1 S 11 = 11-1 SPRINT 
USING '*#####’*; TIME ; SPRINT TAB ( 100 + 11 *10) ; S PRINT USING **### 
#####•*; I , T < I ); SPRINT CHRS<13); 

4362 NEXT 

4365 TWD7=TWD S TWD=ABS ( (T5-W5) - (T6-W6) ) S FLG=1.5 
4370 IF TWD >FLG AND TWD7<=FLG THEN GOTO 4390 
4375 IF TWD<=FLG AND TWD7 >FLG THEN GOTO 4390 
4380 IF IPFLG=1 THEN PRINT 
4385 RETURN 

4390 PRINT USING "##### ##.## ##.## ##.## ##.## 

## . ###** ; TIME , T5 , W5 , T6 , W6 , TWD 
4395 RETURN 
4400 ' 

4405 ' Home signal processing.* HM2-> t-2, HMl-> t-1, HM-> t 
CHMO-> t+1] 

4410 ' 

4415 

HM2=HM1 S HM1 =HM S HMM2 = HMM1 S HMM1 = HMM : JFLGHO=IFLGHOS JFLGH1 = IFLGH1 S JFL 
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GH2= IFLGH2 

4420 HM=A2<H0ME> SHMM=HM 

4425 IF ABS < HM> <HMAX1 THEN HMM=0 

4430 IF ABS < HM1 ) >HMAX9 AND HMM = 0 THEN HMM = -SGN <HM1 ) *HMAX 
4432 ' 

4435 IF ABS<HMM)=HMAX AND ABS<HMM2> OHMAX THEN IFLGH1=1 
4437 IF ABS(HMM) OHMAX AND ABS CHMM2) = HMAX THEN IFLGH1=0 
4440 IF HMM=0 THEN IFLGH0=1 ELSE IFLGH0=0 

4445 IF ABS CHM-HM1 ) <HMAX1 AND HM>HMAX9 THEN IFLGH2 = 1 ELSE 
IFLGH2=0 

4450 P1=SGN CHMM) 

4455 IF J=C THEN PRINT: PRINT "SYSTEM ";ISYS A2<19>;" 1984 
", DATES, TIMES;" ";N3M;" RECORDS" 

4460 IF J=C THEN PRINT " Tw dryout / Tw < O 

dT or dTw < 0.02 HOME" 

4465 IF J=C THEN PRINT " TIME T5 Tw5 T6 Tw6 dT-dTw 
dTavg dTwavg dT dTw PI raw" 

4470 ' IPFLG=OSll=0 

4475 * IF IFLGHOO JFLGHO THEN IPFLG=1 :PRINT USING 
"#####"; TIME; :PRINT TAB<74*I1*9) ; SPRINT USING "### 

#####";P1,HM; SPRINT " lnter*ed";CHRS<13> ; :I1=I1*1 
4480 'IF IFLGHIO JFLGH1 THEN IPFLG= 1 S PRINT USING 
••#####••; TIME; SPRINT TAB <74*11 *9) ; S PRINT USING "### 

#####" ; PI , HM ; S PRINT " bad aw " ;CHRS( 13) ; S 11=11*1 
4485 'IF IFLGH20 JFLGH2 THEN IPFLG=1 SPRINT USING 
••#####•• -TIME; SPRINT TAB<74 + I1«9) ; SPRINT USING "### 

#####»;P1,HM; SPRINT ” Homed ";CHRS(13> ; :I1=I1+1 
4490 RETURN 
5000 ' 

5005 ' CONVERT RTD READINGS TO DEG. C. 

5010 ' 

5045 T= <A2 (1C ) ) /RC < 1 ) /G <CH < IC) ) 

5050 A2<IC) =-245. 665+T- <235.476*10. 189*T) 

5065 RETURN 
6200 ' 

6205 ' read system specific data 
6207 ' 

6210 ISYS=VAL<MIDS<N1S,4, 1) ) S IF ISYS=9 THEN HMAX=300O ELSE 
HMAX=1000 

6215 HMAX9=.9*HMAX S HMAX1 = . 1 *HMAX S HMAX=5000 

6217 JFLGH0=0 : IFLGH0=0 S JFLGH1=0 S IFLGH1 = 0 S JFLGH2 = 0 S IFLGH2=0 

6218 

IFLGW7=0 S IFLGW = 0 S IFLGDT7=0 S IFLGDT = 0 S IFLGD7=0 S IFLGD = 0 S TWD7 = 0 S TWD = 0 

6219 FOR 1=1 TO Ml S IFLGO < I ) =0 S IFLG07 < I > =0 S NEXT 

6220 ON ISYS GOSUB 6230,6280,6280,6280,6280,6280,6245,6260,6275 
6222 PRES =101 .3- . 01055»ELEV ' ASSUME STD ATMOSPHERE, ELEV = 
ELEVATION <M> 


6223 

6230 

RETURN 
VANE = 

8.3 S 

VDATE=919 

S VTIME= 600 

6235 

6245 

RETURN 
VANE = 

5.2 S 

VDATE=915 

S VTIME= 920 

6250 

RETURN 
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6260 VANE = 11.3 : VDATE=915 : VTIME= 1800 
6265 RETURN 

6275 VANE = 11.5 : VDATE=915 : VTIME= 1250 
6280 RETURN 
6300 ' 

6305 INFLSs="INDAT"-*-RIGHTS (STRsS < ISYS) , 1 > ♦ " . DO" 

6310 OPEN "I”, #1,INFLS:NDIG=0:NRTD=0 

6312 INPUT #1,X9S ' SKIP LABEL 

6315 INPUT #1,M,N1,N2,N3,N4,N5,N8,G0,M7 

6320 N4=N4/N1 ' SET N4=# OF RECORDS /DISK UPDATE 

6335 INPUT #1,X9S3 ' SKIP LABEL 

6336 INPUT #1 ,LG! , HG , HOME , REF , 01 ,02,RC<1) ,NCRTD 

6337 INPUT #1,X9S ' SKIP LABEL 

6338 INPUT #1,DELZ<1),ELEV,CS0IL,DZ, REF, HOME 

6339 INPUT #1,X9S 

6340 FOR K8 = l TO M7 .1, 

6345 INPUT #1 , CH <K8 ),C1<K8),G <K8) , B <K8 > , NT (K8) , X9S 
6350 IF C1(K8)=0 THEN G2<K8)=1/LG! ELSE G2<K8)=1/HG 
6355 IF Cl <K8) =1 THEN 81<K8)=01 

6360 IF Cl C K8) =3 THEN B1(K8)=02 ^ 

6365 Cl <K8> =C1 <K8) *16-*-CH <K8) -1 
6368 IF NT <K8) =2 THEN NRTD=NRTD+1 
6370 IF NT <K8 ) -3 THEN NWD ~K& 

6375 IF NT <K8) =5 THEN NDIG=NDIG+1 

6380 NEXT 

6395 NANLG=M7-NDIG 

6399 CLOSE #1: RETURN 

6400 ' 

6405 ' create output directory file 
6408 ' 

6410 Q5=4 : N= J-l : M = L1 

6425 GOSUB 32100 : DSs = "Manhattan" -TIMEsS-*- ” " ♦DATES 
6430 OPEN "0",#1,N3$ 

6440 PRINT#l,Q5;",";N;”,";M;",";DiS;","; :F0R L=1 TO LI: PRINT 

#1 , AS ( L <L) ) :NEXT 

6450 PRINT #1,ZS : CLOSE #1 

6460 PRINT FS5:NS = N3S: GOSUB 7100:PRINT 

6495 RETURN 

6700 ' 

6705 ' MV TO ENG. UNITS, LINEAR 

6710 ' 

6715 A2 < IC) =A2 < IC) **G < IC) + B < IC) 

6800 RETURN 
7000 ' 

7060 IF LEN <G$) =0 THEN G = 1:PRINT GlRETURN 
7070 G=VAL<GS) : PRINT : RETURN 
7090 ' 

7100 PRINT"HEADER DATA FOR: ";NS TABOO) "LABEL: " D$ 

7110 PRINT" NUMBER OF CASES: " N TABOO) "NUMBER OF VARIABLES: " 
M: RETURN 
7120 ' 

7200 ON ERROR GOTO 7250 
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7210 OPEN "I",#1,NS:INPUT #1,Q5,N,M,DS 
7220 FOR J=1+M1 TO M+M1:INPUT #1 , AS< J> :RSET 
SPS=AS(J) :AS<J)=SPS:NEXT J:INPUT #1,ZS 
7230 CLOSE #l:ON ERROR GOTO O : RETURN 

7250 PRINT: IF ERR=53 THEN PRINT "FILE NOT FOUND" : PRINT JS 

7255 IF ERR053 THEN PRINT "ERROR # " ; ERR ; " IN LINE " ; ERL 

7260 INPUT "NEW FILE NAME :", NS : NS=HS+N$ : CLOSE #1 

7270 GOSUB 8300 

7280 GOTO 7210 

7300 ' 

7400 PRINT: INPUT; "ENTER BEGINNING CASE NUMBER: " ,C 
7410 INPUT", ENDING CASE NUMBER: ",D 
7420 G=C:H1 = 1 :H2=D:G0SUB 8200: IF WOl THEN 7440 
7430 PRINT JS:GOTO 7400 

7440 G=D:H1=C:H2=N:G0SUB 8200:IF WOl THEN RETURN ELSE 7430 
8000 ' 

8010 ' *S-R* 

8020 ' 

8030 PRINT 

8035 PRINT"ENTER OPTION: " , : GS=INPUTS < 1 ) 

8040 IF ASC(GS) =13 THEN GS=MIDS (T1S, 1 , 1 ) 

8050 G=ASC(GS) -64: PRINT GS; 

8060 HI * ASC < LEFTS (T1S , 1 ) ) -64 : H2= ASC < MIDS (T1S , 2,1)) -64 : GOSUB 8200 
8080 IF WOl THEN RETURN ELSE 8035 
8090 ' 

8200 IF G> = H1 AND G<=H2 THEN W=0:RETURN 
8210 PRINT JS;CHRS< 13) ; :W=1 : RETURN 
8215 ' 

8300 OPEN "R" , #1 , "PARMD" , 38 

8310 FIELD #1,19 AS XS,9 AS NNS : GET #1,1:LSET XS=X$:LSET 
NNS=NS : PUT #1,1 
8320 CLOSE #1: RETURN 
8325 ' 

8400 IF LEFTS <NS,6)="< NONE ) " THEN 8430 
8410 PRINT:PRINT"OPEN FILE: " CHRS<34) NS CHRS<34);" 

8420 PRINT" (PRESS " CHRS<34) "RETURN" CHRS<34) " TO USE 
FILE)" 

8430 PRINT"ENTER FILE NAME: " ; : N9S="" : FOR J=1 TO 10 
8432 XXS=INPUTS<1) : IF XXS=CHRS(13) THEN 8438 ELSE PRINT 
8434 N9S=N9S+XXS 
8436 NEXT 

8438 IF LEN < N9S) =0 THEN PRINT" " NS : PRINT : RETURN 

8439 IF MIDS (N9S ,2,1)=":" THEN NS=N9S:G0T0 8450 

8440 NS=HS+N9S 

8450 GOSUB 8300: PRINT" NS : PRINT : RETURN 

8455 ' 

8600 PRINT :PRINT TAB (10) " VARIABLE NUMBERS AND NAMES 

": PRINT 

8620 A=A+6:B=B+6:IF B>M THEN B=M 
8630 FOR K= A TO B 

8640 PRINT USING "###"; K; :PRINT" . " AS(K) " ";:NEXT K:PRINT:IF 
B<M THEN 8620 


OPEN 

XXS; 
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8670- RETURN 
8675 ' 

8900 D9$=”N0 YES“ : IF D9 = l THEN D9S="YESN0 ": PRINT 
8910 PRINT Q1S; : QS=INPUTS < 1 > 

8920 IF QS=MIDS<D9S,4,1> THEN PRINT RIGHTS (D9S , 3 ): 0=1 : RETURN 
8930 PRINT LEFT$(D9S,3) :Q=0:RETURN 
8935 ' 

8950 PRINT: PRINT " " • 

8960 FOR L=1 TO LI : PRINT " " AS<L<L> > ; :NEXT 

L : PRINT : PRINT : RETURN 
8965 ' 

8970 PRINT :PRINT"PRESS 'RETURN' TO CONTINUE: 

: QS= INPUTS < 1 ) : L3 = l : PRINT FS 
8980 RETURN . - 

8985 ' 

8990 IF P=1 THEN RETURN 

8992 PRINT :PRINT"PRESS ANY KEY TO CONTINUE: " ; : QS = INPUTS < 1 > 

8994 PRINT CHRS < 13) ;: RETURN . y 

9000 ' 

9010 ' » INIT* 

9020 ' ^ 

9120 RS=CHRS<13)+" 

9300 ' 

9305 ' * INITIALIZE CONTROL PARAMS * 

9310 ' 

9330 DPR=57 .2958 ' DEGREES/RADIAN 

9355 SIGMA=5 . 6697E-08' BOLTZMAN CONST 

9365 PI=3. 14159 

9370 OPEN "I" ,#1 , “PDS . FIL" 

9380 INPUT #1,PGS:IF PGSO”SAMPE"THEN 9380 
9517 INPUT #1 , ICFLG, IS, IE, DXIS, DXOS, FSS, FTS , MSGS 
9520 IFILES=0 

9525 I FILES= IFILES +1 : INPUT #1 , FLS < IFILES) : IF EOF (1) THEN 9540 
ELSE 9530 

9530 IF FLS <IFILES> =“END" THEN IFILES=IFILES-1 :GOTO 9540 
9535 PRINT IFILES; FLS < IFILES) GOTO 9525 
9540 CLOSE #1:PRINT IFILES;FL$<IFILES) 

9799 ' 

9800 ' Field (variable) names 
9805 ' 

9810 FOR 1=1 TO M : READ AS<I>:NEXT 

9820 DATA "G ","Q ", "KDN ",“KUP " , **D ","QDN ","HOME 

** , "UDIR ** 

9830 DATA "TAR ","TWR ", "TAL ",“TWL " , "TSOIL" , "THAT ","U 
","OUP ” 

9850 DATA "TIMER" , "DATE " 

9890 RETURN 
32000 ' 

32005 ' Time and date routine for SDS MV5.0c 
10/15/84 11:40 

32007 ' last modified 11/ 

1/84 7:50 
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32010 ' inlt 
32015 ' 

32020 DIM M ( 12) 

32025 FOR 1=1 TO 12: READ M(I):NEXT:IF INT (DATE/4) *4 = DATE THEN 
M ( 2) =29 

32030 DATA 31,28,31,30,31,30,31,31,30,31,30,31 
32035 ' 

32040 ' main program 
32045 ' 

32100 RETURN: " ' 'DATE=PEEK (S.H40) : IF DATE=0 THEN GOSUB 32170 
32105 DATES= “19“ ♦RIGHTS ( STRS < DATE) ,2) 

32120 

TIMES = RIGHTS (STRS (PEEK (S.H43) ) ,2) ♦ “ : “+RIGHTS (STRS (PEEK (6.H44) ) , 2> + “ 
: “+RIGHTS (STRS < PEEK (&H45) ) ,2) 

32125 DA Y= PEEK (&H41 ) +256* (PEEK (&H42) ) : M0=0 

32135 FOR 1 = 1 TO 12 : MO = MO + M ( I ) :IF M0> =DAY THEN MO = MO-M ( I ) : GOTO 
32145 

32140 NEXT 
32145 DAY =DAY- 

MO:MO=I :DATES = R.IGHTS(STRS(MO) ,2) ♦“/“+RIGHTS(STRS(DAY> ,2) ♦“/“♦DATE 

s 

32150 PRINT DATES, TIMES 
32160 RETURN 
32165 ' 

32170 RETURN:' set time 
32175 ' 

32180 INPUT “MONTH, DAY, YEAR ”, MO,DAY, YR : M0=M0-1 

32185 INPUT “HOUR, MINUTE, SECOND “ , HR , MIN, SEC 

32190 JDAY=0:F0R 1=1 TO MO : JDAY= JDAY+M ( I > : NEXT 

32195 JDAY= JDAY+DAY : JDAYH=INT ( JDAY/256) : JDAYL= JDAY- JDAYH-256 

32200 POKE S.H41 , JDAYL : POKE &H42,JDAYH 

32205 POKE &H43,HR:P0KE &H44 , MIN : POKE &H45,SEC 

32210 PRINT “depress <CR> to start clock :XCS=INPUTS<1 ) 

32215 IF XCSOCHRS (13) THEN GOTO 32185 ELSE POKE &H40,YR 
32220 RETURN 


Appendix 9.8.9 SAMPP.BAS, a program for the AT computer which 
converts the output of SAMPE.BAS into 6-minute energy and 
radiation balances. 

6 

6/25/86 0445 

10 ' SAMPB3.BAS RANDY data analysis program 

6/13/85 1330 

20 ' Based on PROGRAM SAMPC3 1/ 

7/85 2230 

30 ' Last modified 

5/5/86 

32 ' For Hanford Site study, Washington 

35 ' Reduced soil heat capacity ala DeVries (1963) 

40 'CONVERT INDAT CSOIL IN XH20(G/G) TO CSOI=HEAT CAPACITY, 
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BATTELLE 4/86 
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55 
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100 

140 
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155 

160 


ICFLG = 


include IS point running mean of G in top 10 


' . 1 -> exclude G calculation in top 10 cm 

'IS = no. of points in soil heat storage running mean 

'FSS = output file name extension ( .MF or R> 

'FTS = not used 
S4=0:S5=0 :S7=0 

DEFINT I,J,L-N : NS=19:NST=52:N0UT=34 
DIM T <54 ) , IFLG0 (30) , IFLG07 <30) 

DIM N <25> , D < 25) ,F<17) ,AS<53) , L <53) , TS <13) , C <50 , 4) , BS <51 ) 
DIM A2 < 20) , CH < 20) , Cl < 20) , G < 20) , B <20) , G2 < 20) 

DIM NT <20 ) ,B1<20) ,FL3<120> ,NS<9) 


165 
170 DIM 

DELZ <2) ,07<2> ,G7<2> , TAV7 <2) , WAV7 <2) , T7 <2) , T8 <2) ,W7<2) ,W8<2) 

175 DIM RC <2> , E < 9) , S < 7) , GP < 10) , RCC< 250) 100) 

180 RC2=0 :RCC=0: 'record counter for loop and the first three 

records , 

315 GOSUB 6100 ’ ' set constants 

1000 GOSUB 9000 :FS="" ' Microstat init 

1010 'GOSUB 32000 ' init ^jlbck 

2000 ' 

2005 ' MAIN PROGRAM 
2010 ' 

2030 M3=0:M2=0:M1=0:S4=0:S5=0:S7=0.*FOR 1 = 1 TO IS: GP < I > =0 :NEXT 
2040 ' 

2045 ' read input file directory 

2055 ICOUNT = ICOUNT + 1 : NSs=DXIS + FLS < ICOUNT) : 

nis=ns+".mf": print "Nss=",nss 

2060 GOSUB 7200:N1M=N:M1=M:PRINT FSIGOSUB 7100:D1$=D$ 


2070 'GOSUB 6200 ' set system 

specific info 
2100 ' 

2540 LI =NOUT 
2565 ' 

2570 C=1 : D=N ' IF ICFLG=1 THEN D=N : PRINT "OUTPUT ALL CASES <" 
N1M •■>•* 

2575 ' IF ICFLG=0 THEN C=IS : D=IE : PRINT "OUTPUT 

SUBSET OF CASES <" IS "TO" IE ">" 

2605 N3M=D-Cf 1 
2608 ' 

2617 N3S=DX0S+FL$< ICOUNT) +"P" :Z$=N3$+FS$ 

2700 ' 

2870 FOR L= 1 TO LI 
2880 L < L) = L«-M1 
2890 NEXT 

2895 GOSUB 6300.1 'READ INDAT? FILES 

2900 GOSUB 6400 ' create output 

file directory 

2990 PRINT FsS:PRINT“FILE: " N3S “ IS NOW BEING OUTPUT ":J1 = 0 

2995 ' 

2996 ' open input and output files 
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2997 ' 

3000 OPEN "R",#1,N1S,05 
3030 OPEN "R" , #2,Z$,Q5 
3050 FIELD #1,Q5 AS TS 
3080 FIELD #2,Q5 AS US 
3100 ' 

3105 ' main computation loop 
3110 ' 

3150 FOR J = C TO D:ON ERROR GOTO 3205 
3155 ' 

3160 ' read data into TO 
3165 ' 

3170 FOR 1 = 1 TO Ml : GET #1 , 1 + (Ml * ( J-l ) > 

3190 A2 ( I ) =CVS ( TS) 

3200 GOTO 3210 

3205 PRINT "ERROR #" ERR " OCCURED IN LINE" ERL " 

3208 ON ERROR GOTO O 

3210 NEXT 

3230 TIME = I NT ( A2(NS-2) 710) 

3427 GOSUB 3600 

3440 Jl=Jl-*-l :PRINT CHRS(13) J1 TIME A2(NS-1> " 

3450 ' 

3455 ' write out full TO array 

3456 ' 

3460 FOR L= 1 TO LI 

3480 LSET US=MKSS(T(L(L> > > 

3510 PUT #2 ,L+ <L1»(J1-1>> 

3520 NEXT 

3530 NEXT:PRINT FS 

3540 PRINT "END OF FILE OUTPUT" : N$=N1S 
3570 CLOSE #1:CL0SE #2:PRINT 
3580 NS=N3S: ' GOSUB 8300 
3590 IF ICOUNT <IFILES GOTO 2000 

3595 ' IF ISTP=0 THEN CHAIN "SUMMARYE" ELSE CHAIN "SC" 

3596 CHAIN "SUMMARYE" : END 
3600 ' 

3605 ' MAIN SAMPLING LOOP 
3610 ' 

3785 PI =SGN < A2 (HOME) ) : REM GOSUB 4400 ' Home signal 

processing 

3890 ' 

3895 J8=NCRTD : 19 = 1 : GOSUB 4000 ' Energy and 

radiation balance 
3898 7 

3900 T ( NS) =TIME : T(NS+1) =QSTAR : T(NS+2> =H : 

T (NS+3) = E 

3910 T ( NS+4 ) =GP : T(NS+5> =KDN : T(NS+6) =KUP : 

T ( NS+7) = A2<5> 

3920 T ( NS + 8 ) =LDN : T(NS+9> =LUP : T(NS-*-10) =A2(15> : 

T (NS+11 ) =A2(8) 

3930 T (NS+12) =T9 : T (NS+13) =W9 

T (NS+15) = W0 
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T (NS+14 > =T0 


: T < NS+ 17 > =E9 
T <NS+21 > = QDN 
: T <NS+25) = P1 


3940 T<NS+16) = A2 (13) 

T <NS+19) =DT 
3950 T <NS+20) =DE 
T <NS+23> =Q9-QN 
3960 T <NS+24 ) =RHB 
T < NS + 27 ) a A2 < HOME ) 

3970 T (NS + 28 ) =A2 <NS-1 ) :T<NS+29)=BR 
T<NS+31 ) =EALT 

3972 T <NS+32) =CV# : T<NS+33)=RB 

3980 RETURN 
4000 ' 

4005 ' Bowen ratio energy balance - 2 point running 


T < NS+ 18 ) =E0 
T <NS+22) =QUP 
T <NS+26 ) -GS 
T <NS+30) = HALT 


mean 

4010 * SUB5,6 = PRESENT VAL. , SUB7,8 = PAST VAL., SUB9,0 = 

RUNNING AVE. ^ 

4015 ' 

4020 Q9=A2<2) :G5=A2<1) - , 

4025 S4=S7 :S7=S5 : S5= A2 < 13) ' S7=Taoil at TIME-6 mins; *S4 at TIME- 
12 mins 

4030 TAV5=<A2<J8)+A2<J8 + 2) )/2:WAV5= <A2<J8 + 1) +A2< J8 + 3)''3#2 
4035 IF Pl=l THEN IALT=0 ELSE IALT=2 

4040 T5 = A2 < J8 + IALT) : T6=A2 < J8-IALT+2) : W5 = A2 < J8+IALT+1 ) :W6=A2<J8- 
IALT +3) 

4051 ' 

4052 ' *** FIND RUNNING AVERAGES *** 

4053 ' 

4054 QSTAR= <Q7 < 19) +09) /2 : GP= <G7 < 19) +G5) 72 : QN9= <QN7 + QN5)72 

4055 T= ( TAV5 + TAV7 < 19) ) 72 : TW= < WAV5+WAV7 < 19 ) )72 

4059 

T9= < T5+T7 <I9))72:T0=(T6+T8<19) ) 72 1 W9= < W5+W7 <I9))/2: W0= < W6+W8 ( 19) ) 
72 

4060 IF 1 > RCC THEN GOSUB 5670 :' STARTUP AVERAGES 

4061 IF QSTAR < O THEN QSTAR=1 . 0621 *QSTAR 

4062 IF ISYS=1 THEN QSTAR=QN9 : ' THR NET 

4063 'GOSUB 4300 ' wet bulb processing 

4064 ' 

4065 ' *** SAVE PRESENT VALUES *** 

4066 ' 

4070 G7<19) =G5;Q7(I9) =09:QN7=QN5 : 'THR NET 

4075 TAV7 < 19) =TAV5 : WAV7 < 19) =WAV5 

4080 T7 < 19) =T5 :T8 < 19) =T6 : W7 < 19) =W5 : W8 < 19 ) =W6 

4084 ' 

4085 ' *** MISCELLANEOUS PARAMETERS **« 

4086 ' 

4090 TT-T : W1 =TW : GOSUB 6015:EA=EFN 

4092 CP= (239 .9+440 . 9* . 622 *EA/ < PRES-EA ) ) / . 2388 

4095 XL=2501300t-2366*TW: GOSUB 6030:SO=S 

4100 G4=PRES*CP7 < . 622* XL) : R0=3 . 4838* < PRES- . 378*EA) /< T+273 . 16) 

4105 Sl=9. 810001/CP : TT=TW : GOSUB 6030 

4110 S2 = 9 .810001 * ( 1/CP+ .0034857*EA/ <273. 16 + T) /G4) 7 < 1 +S/G4) 

4115 S3 s 9 . 810001 « .O034857*EA7 <273. 16 + T) 

4119 * 
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4120 ' *** GRADIENTS *** 

4121 ' 

4125 TT=T9:W1=W9: GOSUB 6015 : E9=EFN : W1 =T9 : GOSUB 
6015 : RHT=100»E9/ESAT 

4126 TT=T0:W1=W0: GOSUB 6015 : EO=EFN : W1 =T0 : GOSUB 
6015 : RHB=10O»E0/ESAT 

4130 DT=T9-T0«-S1 *DELZ ( 19 ) 

4135 DE=E9-E0+S3*DELZ(I9) 

4159 ' 

4160 ' *** BOWEN RATIO USING T, E *** 

4161 RCC=RCC+1 

4162 'Convert #H20<G/G) to volumetric and calc heat capacity. 

4164 GS = -CS0I*DZ« <S5-S4) / <2*N1*60) : IF RCC <3 THEN GS=0:'heat 
storage. 

4166 BR = G4 « DT / DE : Q A V = QST AR + GP + GS 
4168 H = (-QAV)/ (li-l/BR) :E = H/BR 

4170 GOSUB 5005 

4171 IF SGN <E) <> SGN <DE) THEN E=EALT : H=HALT 

4172 IF (-.6 > BR) AND < BR > -1.25) THEN E = EALT : H = HALT 

4174 ' *** RADIATION BALANCE *** 

4175 ' 

4177 KUP=-A2<4) :KDN=A2<3) 

4180 IF KDN<=0 THEN A=0 ELSE A=-KUP/KDN ' ALBEDO 
4200 f 

4205 ' Diffuse correction, per LI-COR 2010S shadow band manual 
4210 ' NOTE: Eppley and not LI-COR used for total solar radiation 
4215 ' 

4220 IF KDN< =0 THEN 4245 ELSE A2 (5) =A2 <5 ) * 1 . 13 ' Table I 

4225 A2<5)=A2<5)/<1.17-<1/<1.2+11.8*CA2<5)/KDN) ))) ' Spectral 

correction 
4235 ' 

4245 IF KDN<0 THEN KDN=0 

4250 IF KUP>0 THEN KUP=0 

4255 IF A2 <5) <0 THEN A2<5)=0 

4257 QUP=-SIGMA* (A2<14) +273 •16) / '4-A2<16) 

4260 QDN = SIGMA« < A2 < 14 ) +273 . 16) ~4-*-A2 < 6 > 

4261 IF QDN > 3000 THEN QDN=3000:IF QUP < -3000 THEN QUP=-3000 

4262 IF KDN = 0 THEN QUP = 1.062-QUP 

4263 IF KDN = O THEN QDN = 1.062*QDN 

4265 LUP=QUP-KUP:LDN=QDN-KDN:QN=QDN+QUP:QN5=QN 
4270 IF ISYS=7 THEN GOSUB 5650 
4280 RETURN 
4300 ' 

4305 ' Checks: Tw < 0, dT or dTw < .005, Tw -> T and 4095 < 
signal < -4095 
4310 ' 

4330 IFLGW7 = IFLGW : IF W5<0 OR W6<0 THEN IFLGW = 1 ELSE IFLGW=0 
4335 IF IFLGWOIFLGW7 THEN I PFLG = 1 : LPRINT USING 

; TIME ; :LPRINT TAB < 13) LPRINT USING "##.## 

W5, W6 ; :LPRINT CHRS(13); 

4339 ' 

4340 IFLGDT7=IFLGDT: IF ABS(T9-T0) < .005 OR ABS<W9-W0) < .005 THEN 
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IFLGDT=1 ELSE IFLGDT=0 

4345 IF IFLGDTOIFLGDT7 AND J>C THEN IPFLG=1 rLPRINT USING 
•*#####••; TIME ; rLPRINT TAB < 39 >; r LPRINT USING "##.#### ##.####"; T9- 
T0,W9-W0; rLPRINT CHRS<13>; 

4346 ' 

4347 IFLGD7= IFLGD : IF ABS (T5-T6) < .02 OR ABS < W5-W6) < .02 THEN 
IFLGD= 1 ELSE IFLGD=0 

4348 IF IFLGD OIFLGD7 AND J>C THEN IPFLG= 1 : LPRINT USING 
••#####•« ;TIME; rLPRINT TAB (57) ; .‘LPRINT USING "##.#### ##.####";T5- 
T6,W5-W6; rLPRINT CHRS<13>; 

4349 ' 

4355 Il = -l:FOR 1 = 1 TO Ml-3 : IFLG07 C I > =IFLGO < I > .* IF ABS ( A2 ( I > ) >=4095 
THEN IFLGO < I ) =1 ELSE .IFLGO < I ) =0 

4360 IF IFLGO < I ) OIFLG07 < I > AND J>C THEN IPFLG = 1 : 11 = 11 + 1 rLPRINT 
USING "#####"; TIME; rLPRINT TAB<100+I1»103 VrLPRINT USING "### 
#####••; I, a 2( I) ; rLPRINT CHRS<13> ; 

4362 NEXT . . ' -* 

4365 TWD7=TWD r TWD=ABS< (T5-W5) - <T6-W6> > r FLG=1.5 
4370 IF TWD >FLG AND TWD7<=FLG THEN GOTO 4390 
4375 IF TWD<=FLG AND TWD7 >FLG THEN GOTO 4390 
4380 IF IPFLG=1 THEN LPRINT 
4385 RETURN 

4390 LPRINT USING "##### ##.## ##.## ##.## ##.## 

##.###"; TIME, T5,W5,T6,W6, TWD 
4395 RETURN 
4400 ' 

4410 ' 

4415 

HM2=HM1 r HM1 =HM r HMM2 = HMM1 r HMM1 = HMM r JFLGHO= IFLGHO r JFLGH1 = IFLGH1 r JFL 
GH2=IFLGH2 

4420 HM=A2(HOME) ;HMM=HM l GET #1,7+M1»J : HMO = CVS<Ts3> 

4425 IF ABS<HM> <HMAX1 THEN HMM=0 

4430 IF ABS < HM1 > >HMAX9 AND ABS <HMO) >HMAX9 AND HMM = 0 THEN HMM = - 
SGN < HM1 ) *HMAX 
4432 ' 

4435 IF ABS<HMM) =HMAX AND ABS <HMM2) OHMAX THEN IFLGH1=1 
4437 IF ABS(HMM) OHMAX AND ABS <HMM2) =HMAX THEN IFLGH1=0 
4440 IF HMM=0 THEN IFLGHO=l ELSE IFLGH0=0 

4445 IF ABS < HM-HM1 ) <HMAX1 AND HM>HMAX9 THEN IFLGH2=1 ELSE 

IFLGH2=0 

4450 P1 = SGN <HMM) 

4455 IF J=C THEN PRINT rPRINT "SYSTEM ";ISYS A2(19);" 1984 
’* , DATES , TIMES; " ";N3M;” RECORDS" 

4460 IF J=C THEN PRINT " Tw dryout / Tw < O 

dT or dTw < 0.02 HOME" 

4465 IF J=C THEN PRINT " TIME T5 Tw5 T6 Tw6 dT-dTw 
dTavg dTwavg dT dTw PI raw" 

4470 IPFLG=0: 11=0 

4475 IF IFLGHOO JFLGHO THEN IPFLG=1 rPRINT USING 
••#####••; TIME ; rPRINT TAB <74 + 1 1 *9) ; r PRINT USING "### 

#####";P1,HM; rPRINT " interred" ; CHRS < 13) 11=11 +1 
4480 IF IFLGH1 <> JFLGH1 THEN IPFLG=1 r PRINT USING 


dT-dTw 
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•*#####••; TIME ; :PRINT TAB (74 + 1 1 *9) PRINT USING **### 

#####"; PI, HM; :PRINT •• bad aw ; CHRS < 13) 11 = 11 + 1 

4485 IF IFLGH2< > JFLGH2 THEN IPFLG= 1 : PRINT USING 
"#####"; TIME; :PRINT TAB(74+I1*9) ; :PRINT USING "### 

#####"; PI , HM ; :PRINT Homed CHRS< 13) 11=11 +1 

4490 RETURN 

5005 'ALTERNATE CALCULATIONS OF H AND E 

5006 WS = A2 < 15) : 'PRINT “WS=",WS 

5008 CV# = - <QAV) / ( <WS*CP*DT) + < WS*XL* . 622*D£/PRES ) ) 

5501 RB =9 . 810001 *DT* 3 . 24/ < (TT+273) *WS~2) ; '3.24= (2-ZO) ~2 

5506 PRINT ••CV="*,CV#,"RB= ,, ,RB 

5508 IF RB > .006 THEN GOTO 5515 

5510 CVA#=-2.567*RB + .0246: GOTO 5540 

5515 CVA#=-.0123*RB + .0246 

5540 HALT=CVA#*WS*CP*DT 

5550 EALT=CVA#*XL*WS* .622*DE/PRES 

5615 PRINT "HALT=" , HALT, **EALT= ” , EALT 

5634 RETURN 

5650 QUP=A2<16) :QDN=A2(6) :LUP=0:LDN=0: 'QUP AND QDN ARE 

PRYANOMETERS 

5660 RETURN 

5670 'STARTUP AVERAGES 

5671 QSTAR=09:GP=G5:QN9=0N5:T=TAV5:TW=WAV5 

5672 T9=T5 : T0=T6 : W9 = W5 : WO=W6 : RETURN 
6000 ' 

6005 ' • * • MISCELLANEOUS FUNCTIONS * * « 

6010 ' 

6015 

ESAT =<E(1)+W1»(E(2)+W1»<E(3)+W1« <E<4) +W1 * < E (5) + W1 * (E <6 ) +W1 * < E < 7 ) ) 
))))))/10 

6020 EFN = ESAT- . 00066* < 1+ . OOH5*Wl ) *PRES* (TT-W1 ) 

6025 RETURN 
6030 

S= <S(1)+TT* <S<2) +TT * <S < 3 ) +TT * ( S < 4 ) +TT « < S < 5 ) +TT *» (S ( 6) +TT* <S ( 7 ) ) ) ) ) 
) ) )/10 
6035 RETURN 
6050 ' 

6100 ' * * * MISCELLANEOUS CONSTANTS * • * 

6105 ' 

6115 E < 1 ) =6 . 1078 

6116 E<2)=. 44365185# 

6117 E<3)=. 014289458# 

6118 E<4) = .00026506485# 

6120 E <5) =3 .031240400000003D-06 

6121 E< 6) =.000000020340809# 

6125 E<7) =6 . 136820900000059D-11 

6126 ' 

6130 S(l)=. 44381 

6131 S < 2) = .028570026# 

6132 S (3) =7 . 93805E-04 

6133 S(4) =.000012152151# 

6135 S(5)=. 00000010365614# 
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6136 S<6> =3 . 532421800000003D-10 
6140 S(7)=-7. 090244800000164D- 13 
6150 RETURN 
6200 ' 

6205 ' read system specific data 
6207 ' 

6210 ' ISYS=VAL<MID3(N1S,4,1) > : 'IF ISYS=9 THEN HMAX=15 ELSE 

HMAX=5 

6215 ' HMAX9= .9*HMAX : HMAX1= . 1 *HMAX : HMAX=30 

6217 ' JFLGH0=0 : IFLGH0=0 : JFLGH1 =0 : IFLGH1 =0 : JFLGH2=0 : IFLGH2=0 

6218 ' 

IFLGW7=0 : IFLGW=0 : IFLGDT7=0 : IFLGDT=0 : IFLGD7=0 : IFLGD=0 : TWD7=0 : TWD=0 
6219 ' FOR 1=1 TO Ml : IFLGO < I ) =0 : IFLG07 ( I > =0 : NEXT 
6300 ISYS=VAL<MIDS<Ni$,4,l>) 

6305 INFLS="INDAT"fRIGHT$(STR3<ISYS) , 1 ),+?•. DO" 

6310 OPEN "I", #1 , INFLS : NDIG=0 : NRTD=0 

6312 INPUT #1,X9S ' SKIP LABEL ' , 

6315 INPUT #1 ,MA,N1 ,N2,N3,N4,N5,N8,G0,rf7 

6320 N4=N4/N1 ' SET N4 = # OF RECORDS/DISK UPDATE 

6335 INPUT #1,X9$ ' SKIP LABEL ‘V??- 

6336 INPUT #1 , LG , HG, HOME , REF, 01 , 02 , RC < 1 > , NCRTD 

6337 INPUT #1 , X93 ' SKIP LABEL 

6338 INPUT #1 , DEL2 ( 1 ) , ELEV, CSOIL, DZ , REF , HOME 

6339 INPUT #1,X9$ 

6340 CLOSE #1 

6342 CSOI= < .402*2+4. 23-CSOIL) *10-6: 'CONVERT XH20 TO HEAT CAPACITY 
6344 PRES=101. 3-. 01055-ELEV: 'ASSUME STANDARD ATMOSPHERE 
6395 RETURN 
6400 ' 

6405 ' create output directory file 
6408 ' 

6410 Q5=4 : N=N3M : M=L1 

6425 GOSUB 32100 : D$= "Manhattan "+TIME5+" "+DATES 
6430 OPEN "0",#1,N3$ 

6440 PRINT#1 ,05 ;",**• N ; ", **;M ; D$ ; " , " ; : FOR L = 1 TO LI : PRINT 

#1,AS<L<L> > : NEXT 

6450 PRINT #1,2$ : CLOSE #1 

6460 PRINT F$:N$=N3$: GOSUB 7100:PRINT 

6495 RETURN 

7000 ' 

7060 IF LEN <G$) =0 THEN G=l:PRINT GrRETURN 
7070 G=VAL<G$> : PRINT; RETURN 
7090 ' 

7100 PRINT"HEADER DATA FOR: ";NS TAB (30) "LABEL: " D$ 

7110 PRINT"NUMBER OF CASES: " N TABOO) "NUMBER OF VARIABLES: " 

M: RETURN 
7120 ' 

7200 ON ERROR GOTO 7250 

7210 OPEN "I",#1,N$:INPUT #1 , 05 , N , M , DS 
7220 FOR J = 1 ♦Ml TO M+M1:INPUT #1 , A$ < J) : RSET 

sp$=a$(J) :a$<j)=sps:next j:input # 1 , 2 $ 

7230 CLOSE #i:ON ERROR GOTO O : RETURN 
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1111 





7250 PRINT : IF ERR=53 THEN PRINT "FILE NOT FOUND" : PRINT JS 

7255 IF ERR<>53 THEN PRINT "ERROR # ";ERR; " IN LINE " ; ERL 

7260 INPUT "NEW FILE NAME :", NS : NS=HS+NS : CLOSE #1 

7270 ' GOSUB 8300 

7280 GOTO 7210 

7300 ' 

7400 PRINT: INPUT; "ENTER BEGINNING CASE NUMBER: ",C 
7410 INPUT" , ENDING CASE NUMBER: ",D 
7420 G=C : HI =l:H2=D: GOSUB 8200:IF W<>1 THEN 7440 
7430 PRINT JS:GOTO 7400 

7440 G=D:H1=C:H2=N:G0SUB 8200.* IF WOl THEN RETURN ELSE 7430 
8000 ' 

8010 ' *S-R* 

8020 ' 

8030 PRINT 

8035 PRINT"ENTER OPTION: " ; :GS=INPUTS < 1 > 

8040 IF ASC<GS>=13 THEN GS=MIDS<T1S, 1 , 1 > 

8050 G=ASC<GS> -64: PRINT GS; 

8060 HI = ASC < LEFTS < T1S , 1 > > -64 : H2= ASC < MIDS < T1S , 2 , 1 > ) -64 : GOSUB 8200 
8080 IF WOl THEN RETURN ELSE 8035 
8090 ' 

8200 IF G>*H1 AND G<=H2 THEN W=0:RETURN 
8210 PRINT JS;CHRS<13> ; :W=1 :RETURN 
8215 ' 

8300 OPEN "R" , #1 , "PARMD" , 38 

8310 FIELD #1,19 AS XS,9 AS NN$:GET #l,l:LSET XS=XS:LSET 
NNS=NS : PUT #1,1 
8320 CLOSE #1: RETURN 
8325 ' 

8400 IF LEFTS <NS , 6) =" (NONE) " THEN 8430 

8410 PRINT :PRINT"OPEN FILE: " CHRS(34) NS CHRS<34);" 

8420 PRINT" (PRESS ” CHRS<34> "RETURN" CHRS<34) “ TO USE OPEN 
FILE)" 

8430 PRINT"ENTER FILE NAME: ” ; : N9S="" : FOR J=1 TO lO 

8432 XXS=INPUTS<1) :IF XXS=CHRS<13> THEN 8438 ELSE PRINT XXS; 

8434 N9S=N9$+XXS 
8436 NEXT 

8438 IF LEN <N9S) =0 THEN PRINT" " NS : PRINT : RETURN 

8439 IF MIDS <N9S, 2, 1 ) =" : " THEN NS = N9S:G0T0 8450 

8440 NS=HS+N9S 

8450 GOSUB 8300: PRINT" NS: PRINT: RETURN 

8455 ' 

8600 PRINT :PRINT TAB<10> " VARIABLE NUMBERS AND NAMES 

PRINT 

8620 A= A+6 : B=B-*-6 : IF B>M THEN B=M 
8630 FOR I=A TO B 

8640 PRINT USING “###"; I ;: PRINT" . " AS<I) “ ";:NEXT I: PRINT: IF 
B<M THEN 8620 
8670 RETURN 
8675 ' 

8900 D9S="N0 YES": IF D9=l THEN D9S="YESN0 ": PRINT 
8910 PRINT Q1S; : QS=INPUTS < 1 ) 
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8920 IF QS=MIDS(D9S,4,1) THEN PRINT RIGHTS(D9S, 3) :Q=1 : RETURN 
8930 PRINT LEFTS (D9S, 3) :Q=0: RETURN 
8935 ' 

8950 PRINT: PRINT “ 

8960 FOR L=1 TO LI : PRINT " " AS ( L <L) ) ; : NEXT 

L : PRINT : PRINT : RETURN 
8965 ' 

8970 PRINT : PR I NT" PRESS 'RETURN' TO CONTINUE: 

**; : QS= INPUTS < 1 ) :L3=1:PRINT FS 
8980 RETURN 
8985 ' 

8990 IF P=1 THEN RETURN 

8992 PRINT :PRINT"PRESS ANY KEY TO CONTINUE: ";: QS = INPUTS < 1 ) 

8994 PRINT CHRS ( 13) RETURN 
9000 ' 

9010 ' *INIT* 

9020 ' / 

9120 RS=CHRS (13) 

9300 ' 

9305 ' « INITIALIZE CONTROL PARAMS * 

9310 ' 

9315 ' NCRTD=9 ' Channel number of 1st RTD 

9320 ' Nl=6 ' Basic data rate (min) 

9330 DPR = 57 . 2958 ' DEGREES/RADIAN 

9355 SIGMA=5. 6697E-08 ' BOLTZMAN CONST 

9365 PI=3. 14159 

9370 ' DZ=.l ' depth of Ts avg <m) 

9380 ' H0ME=7 ' HOME CHANNEL 

9500 ' 

9505 ' Initialize 
9510 ' 

9515 OPEN "I” , #1 , '*PDS . FIL” 

9516 INPUT #1,PGS:IF PGSO"SAMPP" THEN 9516 

9517 INPUT #1,ICFLG,IS,IE,DXIS,DX0S,FSS,FTS,MSGS 
9520 IFILES=0 

9525 IFILES=IFILES+1 : INPUT #1 , FLS < IFILES) : IF EOF (1) THEN 9540 
ELSE 9530 

9530 IF FLS<IFIL£S) =“END" THEN IFILES=IFILES-1 :GOTO 9540 
9535 PRINT IFILES ; FLS ( IFILES ),: GOTO 9525 
9540 CLOSE #l:PRINT IFILES; FLS ( IFILES) 

9799 ' 

9800 ' Additional Field (variable) names 
9805 ' 

9810 FOR I=NS TO NST : READ AS(I):NEXT 

9820 DATA "TIME ",**0 ","H ","E ","GP ** , ”KDN 

•• , ”D 

9830 DATA "LDN ","LUP ","U “,"UDIR 
■* , "TATOP" , **TWTOP“ , "TABOT" , "TWBOT" 

9840 DATA "TSOIL" , "EATOP" , "EABOT" , "DT ","DE ","QDN 

" , **QERR " 

9850 DATA "RHBOT","M ",”GS " , "HMREC" , "DATE " BR 

" , "EALT " 


", "KUP 

" , "QUP 
", "HALT 


75 


9860 DATA "CV ","RB 
9890 RETURN 
32000 ' 

32005 * Time and date routine for SDS MV5.0c 
10/15/84 11:40 

32007 ' laat modified 11/ 

1/84 7:50 

32010 ' init 
32015 ' 

32020 DIM M < 12) 

32025 FOR 1=1 TO 12: READ M(I>:NEXT:IF INT < DATE/4) *4 = DATE THEN 
M ( 2) =29 

32030 DATA 31,28,31,30,31,30,31,31,30,31,30,31 
32035 ' 

32040 ' main program 
32045 RETURN 

32100 DATE=PEEK < &H40) : IF DATE=0 THEN GOSUB 32170 
32105 RETURN: ' DATES =" 19" ♦RIGHTS (STRS( DATE) ,2) 

32120 

TIMES=RIGHTS (STRS (PEEK (&H43) ) , 2) + " : "♦RIGHTS (STRS C PEEK < &H44) ) ,2) ♦ " 
: "♦RIGHTS (STRS (PEEK (S.H45) ) ,2) 

32125 DAY=PEEK(&H41) +256* (PEEK (&H42) ) : M0=0 

32135 FOR 1=1 TO 12 : MO=MO+M ( I ) : IF MO>=DAY THEN MO=MO-M ( I ) : GOTO 
32145 

32140 NEXT 
32145 DAY=DAY- 

MO:MO=I:DATES=RIGHTS(STRS(MO) ,2) ♦"/"♦RIGHTS (STRS (DAY) ,2) ♦“/"♦DATE 
S 

32150 PRINT DATES, TIMES 
32160 RETURN 
32165 ' 

32170 ' set time 
32175 RETURN 

32180 INPUT "MONTH, DAY, YEAR ", MO, DAY , YR : MO=MO-l 

32185 INPUT "HOUR, MINUTE, SECOND ", HR, MIN, SEC 

32190 JDAY=0 : FOR 1=1 TO MO : JDAY= JDAY+M ( I ) : NEXT 

32195 JDAY=JDAY^DAY: JDAYH= INT ( JDAY/256 ) : JDAYL= JDAY-JDAYH*256 

32200 POKE &H41, JDAYL.'POKE &H42,JDAYH 

32205 POKE &H43, HR : POKE &H44 , MIN : POKE &H45,SEC 

32210 PRINT "depress <CR> to start clock ” ; :XCS=INPUT$ ( 1 ) 

32215 IF XCSOCHRS (13) THEN GOTO 32185 ELSE dokE &H40,YR 
32220 RETURN 
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Appendix 9.8.10. SUMMARYE . BAS , a program for the AT computer 
which summarizes the the output of SAMPP.BAS <6-minute data) 
into 30-minute averages. 


5 'SUMMARYE. BAS modified for AT computer and Epson LQ-1000 
printer 

6 

06/25/86 

10 'Program SUMMARYB . BAS 
06/14/85 

11 'Based on SUMMARY5.BAS of 1517 

01/08/85 

-> 'output 30 min averages 
-> input 6 min average/ 
list data as is 


FTS = 


•‘S** 

= ".MF' 
ICFLG = O 


-> 


-> recompute EB w/o soil heat storage, in top 10 
-> recompute EB w/ modified Gs term 




last modified 0919 


12 * 

13 

14 

15 =1 

cm 

20 =2 

1/ 8/85 

25 'NASA Konza Prairie study 
30 
6/16/85 

80 ON ERROR GOTO 30000 

90 DEFINT 1 : WIDTH "LPT1 : “ , 255 : XS="" : GS="" ; NS="" 

100 DIM B*<64) , AS<100) # AIS(IOO) ,D<53) ,FLS<53) ,B5X<53> 

105 DIM 

AVG <34) ,AVG1 <34) , AVG2 <34) ,SUM<34) :TS=CHRS<27) + "S"-*-CHRS <0> :VS=CHRS 

<27)+"T" ' Superscript on/off 

110 DIM IQ<34) ,NQ<34) 

GOSUB 31000 :G0SUB 32000 
IC0UNT=0 


120 

150 

300 

500 

510 

540 


I COUNT = I COUNT >1 : GS=FLS < ICOUNT ) -FTS 
NS=DXIS+GS: GOSUB 7200 : B1X=N : B5X = M : GOSUB 7100 
GOSUB 32100 

601 PI=3 . 1314159# : DPR =57 . 2958 

602 LPRINT CHRS <27) ; "M" ; CHRS <27 ) ; ,, 2" ; ' 12 cpi, 6 lpi 

603 LPRINT TAB<5) DATES " M TIMES " file " NS *• label: " DS 

604 LPRINT: LPRINT MSGS:LPRINT :LPRINT : LPRINT 

605 LPRINT TAB<24) "UNIVERSITY OF WASHINGTON FOREST METEOROLOGY" 

606 LPRINT 

607 LPRINT TAB<24> " ENERGY/RADIATION BALANCE SUMMARY " 

608 LPRINT TAB < 24) " Konza Prairie - 1986" 

609 LPRINT 

610 SYS = VAL <MIDS (GS ,2,1)) :0N SYS GOTO 
611,618,618,618,618,613,615,616,617 

611 LPRINT TAB (24) " System 1: North Facing Slope": GOTO 

618 

613 LPRINT TAB <24) " System 6: Tempe Az demo":G0T0 618 

615 LPRINT TAB < 24) " System 7: East Facing Slope" :G0T0 618 

616 LPRINT TAB (24) " System 8: West Facing Slope" :G0T0 618 
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617 

LPRINT 

TAB (24) 

618 

LPRINT 

TAB (43) 

619 

LPRINT 

USING " 

620 

LPRINT 


621 

LPRINT 

• « 

LDN 

LUP 

Ta 

622 

LPRINT 


623 

LPRINT 

• « 

Wm" 

;TS;"-2 

" ; VS ; " 


System 9: South Facing Slope" :GOTO 618 


TIME Q 
Udir" 


H 


KDN 


KUP 


Wm" ; TsS ; "-2“ ; VS ; " Wm" ; TS ; "-2" ; VS ; " 
;TS;"-2";VS;" Wm" ; TS ; "-2" ; VS; ” Wm";TS;"- 

2" « . •• Wm" • T$ • "— 2** ; V$ • , 

624 LPRINT " Wm" ; TS; “-2" ; VS; " Wm";TS; "-2"; VS;TS; " o";VS;"C“;" 
ms" ; TS ; " - 1 " ; VS ; " deg . " : LPRINT 

625 LPRINT CHRS<27);"0" ' 8 lpi 

626 F1S = ••############ ###### ###### ###### ##### ###### ###### 
##### ###### ###### ###.## ###.# ######" 

627 IF FTS="S" THEN 720 ' create Mstat directory 

628 N1S=DX0S+MIDS<GS, 1 ,LEN<GS) -1 ) +"S" : Q5 = 4 : N = 49 : M = B5X : 


Z0S=N1S+"R" 

629 DS="K0NZA “ + TIMES-*-" "+DATES 

630 OPEN "0",#1,N1S ' Create directory file 

631 PRINT #l,a5;»,";N;",";M;",";DS;","; :F0R J=1 TO M.'PRINT 
#1,AS(J) .‘NEXT J 

632 PRINT #1,Z0S 

633 CLOSE #1 

634 OPEN "R",#1,N1S+"R",Q5:FIELD #1,05 AS T9S 
720 OPEN "R" , #2 , ZS , Q5 : FIELD #2,Q5 AS N9S 
800 FOR RX=1 TO BIX 

802 FOR KX = 1 TO M 

805 GET #2,Kx+M* <RX-1) :D<KX) =CVS(N9S) 

2000 NEXT : GOSUB 10000 

2003 NEXT; CLOSE: FOR 1=1 TO 1 1 : PRINT : NEXT 

2004 GOSUB 33000 

2005 IF ICOUNTCIFILES GOTO 500 

2006 CHAIN "PLOTRE" : END 
2020 GOTO 30020 

7090 ' 

7100 PRINT“HEADER DATA FOR: ";NS TABOO) "LABEL: " DS 

7110 PRINT" NUMBER OF CASES: ” N TABOO) "NUMBER OF VARIABLES: " 
M: RETURN 
7120 ' 

7200 ' 

7210 OPEN "I",#1,NS:INPUT #1,Q5,N,M,DS 

7220 FOR J = 1 TO M:INPUT #1 , AS < J ) : NEXT J:INPUT #1,ZS 

7230 CLOSE #1: RETURN 

7300 ' 

10000 ' 30 min summary routine SUMMARY1.CMF 

10001 ' 9/24/84 Last modified 0914 

10/05/84 

10002 ' 

10003 IF D<1)=0 AND RX=1 THEN RETURN 

10004 FOR 1=1 TO B5X 

10005 AVG<I)=D<I)+AVG(I) 
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10006' NEXT : IF FT$='*S" THEN 10017 

10007 ' 

10008 ' VECTOR AVG WIND DIRECTION 

10009 ' 

10010 A7 = D <125 /DPR 

10011 A1 = A1 + C0S<A7> : A2 = A2-*-SIN<A7) 

10012 IF AlOO THEN A3 = ATN < A2/A1 ) ELSE A3=SGN < A2) »PI/2 

10013 IF SGN(AlXO THEN A3=A3+PI 

10014 IF SGN ( A1 ) >0 AND SGN(A2)<0 THEN A3=A3+2*PI 

10015 AVG < 12) = A3* DPR 

10016 IF D < 1 > =0 THEN D<1>=2400 

10017 HR = INT <D ( 1)/100) 

10018 MIN=D< 1 ) -HR*100 

10019 NR=NR+1 

10020 IF MIN MOD 30 = 0 OR R%=Blx THEN GOTO 10023 

10021 RETURN 

10022 ' , 

10023 XI =X9: X9=HR+MIN/60 : XD= <X9-X1 ) *2 ' XD = # missing records + 
1 

10024 FOR 1 = 1 TO BSx 

10025 IF I<> 12 THEN AVG < I ) =AVG < I > /NR 

10026 NEXT 

10027 IF AVG <6) <0 THEN AVG<6)=0 

10028 IF AVG(8XO THEN AVG<8>=0 

10029 IF AVG <7) >0 THEN AVG<7)=0 

10030 FOR 1=1 TO XD : XT=X1+ .5*1 

10031 AVG <1 ) =40* I NT (XT ) +60* XT 

10032 AVG2 ( 1 > =AVG < 1 ) 

10033 FOR J=2 TO B5X 

10034 AVG2 ( J) =AVG1 < J ) + < AVG ( J ) - AVG1 <J) >/XD*I 

10035 ' IF AVG <1) >1830 AND AVG(l) <2030 THEN PRINT 
AVG(J),AVG1(J) ,AVG2(J) ,XD 

10036 NEXT 

10037 IF IOXD THEN EF$=" *" ELSE EF$ = ,, “ 

10038 IF ICFLG > 1 THEN GOSUB 10100 

10039 GOSUB 20000 : ISUM=ISUM+1 

10040 FOR K= 1 TO B5* : IF FT$="S M THEN 10043 

10041 IF 05=4 THEN LSET T9$=MKS$ < AVG2 (K) ) ELSE LSET 
T9S=MKD$<AVG2<K> ) 

10042 PUT #1 , K+ <B5x* < ISUM-1 ) ) 

10043 IF K<11 OR K=27 THEN SUM (K > =SUM < K ) ♦ AVG2 CK ) * . 0018 ELSE 
SUM <K> =SUM <K) +AVG2 < K) 

10044 NEXT 

10045 NEXT 

10046 FOR K=1 TO B5x 

10047 AVG1 <K) =AVG <K) ' previous averages 

10048 AVG<K) =0 

10049 NEXT 

10050 NR=0:A1=0: A2=0 

10051 IF Rst<Bl JK THEN RETURN 

10052 FOR K = ll TO B5*:IF K<>27 THEN SUM (K > =SUM <K) / I SUM 

10053 NEXT 
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10054 GOSUB 20015 

10055 GOSUB 20017 

10056 GOSUB 20029 

10057 GOSUB 20031 

10058 GOSUB 20035 

10059 GOSUB 20037 

10060 ISUM=ISUM+1 : IF FTsS="S" THEN 10065 

10061 FOR K= 1 TO B5* 

10062 LSET T9$ = MKSsS<SUM<K> ) 

10063 PUT #1,K«-<B5X*<ISUM-1> ) 

10064 NEXT 

10065 FOR K=1 TO B5X :SUM (K) =0 : NEXT : ISUM=0 : X9=0 

10066 CLOSE #l: RETURN 
10100 ' 

10104 IF ICFLG= 1 THEN FACT0R=0 : GOTO 10115 

10105 ON SYS GOTO 10106,10107,10108,10109,10110 

10106 FACTOR= .65: GOTO 10115 

10107 FACTOR= .76: GOTO 10115 

10108 FACTOR= i :GOTO 10115 

10109 FACTORS i ;GOTO 10115 

10110 FACTORS. 62 

10115 AVG2<27> =FACT0R*AVG2 (27) 

10120 BETA = AVG2 <3) /AVG2 (4) 

10125 A VG2 < 4 ) = - < A VG2 < 2 ) ♦AVG2 < 5 ) +AVG2 < 27 ) ) / < 1 ♦ BET A ) 
10130 AVG2 (3) =BETA*AVG2 <4 ) 

10180 RETURN 
10190 ' 

20000 LPRINT USING 

20001 LPRINT USING 

20002 LPRINT USING 

20003 LPRINT USING 

20004 LPRINT USING 

20005 LPRINT USING 

20006 LPRINT USING 

20007 LPRINT USING 

20008 LPRINT USING 

20009 LPRINT USING 

20010 LPRINT USING 

20011 LPRINT USING 

20012 LPRINT USING 

20013 LPRINT EFS 

20014 RETURN 

20015 FL*= 14 : LPRINT 


••############•• ; AVG < 1 ) , 
"###### ,, ;AVG2(2> , 

••######••; AVG2<3) , 

"######" ;AVG2(4> , 

"#####•*; AVG2<5) +AVG2<27> , 

"######" ; AVG2(6) , 

"######*•; AVG2(7) , 

"#####"; AVG2<8) , 

"######"; AVG2<9) , 

"######“; AVG2<10) , 

••###.##••; AVG2<13) , 

"###. #**;AVG2<11) , 

, •######• , ; AVG2 < 12) , 

' removed 


• removed ; RW 


20016 

20017 

RETURN 
FLX= 14 

: LPRINT 

20018 

FLX= 6 

: LPRINT 

USING 

20019 

FLX= 6 

: LPRINT 

USING 

20020 

FLX= 6 

: LPRINT 

USING 

20021 

FLX= 5 

: LPRINT 

USING 

20022 

FLX= 6 

: LPRINT 

USING 


Totals:"; ; 

"###. ##";SUMC2) , 

"### . ; SUM (3) , 

M ###. ##";SUM<4) , 
"##.##••; SUM < 5) +SUM < 27 ) , 
##";SUM<6) , 


RW 
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20023 FLX= 6 ILPRINT USING "###.##”; SUM < 7) , 

20024' FLX= 5 ILPRINT USING "##.##*'; SUM (8) , 

20025 FLX= 6 : LPRINT USING *’### .##**; SUM (9) , 

20026 FLX= 6 ILPRINT USING '*### . ##*' ; SUM ( 10) , 

20027 FLX= 14 I LPRINT ” <MJ m**;TsS; ••-2" ; VS; ") " 
removed ; -RW 

20028 RETURN 

20029 FLX= 14 ILPRINT ' removed ; -RW 

20030 RETURN 

20031 FLX= 14 ILPRINT " Averages (units as in column 

headings) I •*; ; 

20032 FLX= 6 ILPRINT USING "### . ##" ; SUM ( 13) , 

20033 FLX= 5 ILPRINT USING "###.#"; SUM < 11 ) , 

20034 RETURN 

20035 FLX= 14 ILPRINT i'LPRINT ' removed ; -RW 

20036 RETURN ^ ; 

20037 FLX= 14 ILPRINT " <•*; 

20038 FLX= 14 ILPRINT CHRS<34); j 

20039 FLX = 14 ILPRINT 

20040 FLX= 14 ILPRINT CHRS<34); 

20041 FLX= 14 ILPRINT " indicates interpolated values inserted 
for missing data)”; 

20042 

LPRINT I LPRINT I LPRINT I LPRINT I LPRINT I LPRINT I LPRINT I LPRINT I LPRINT I LP 
RINT 

20043 RETURN 
30000 ' 

30010 IF ERROS1 THEN 30030 
30020 CLOSE 

30025 PRINT I PRINT I CHAIN "PL0T4" 

30030 ON ERROR GOTO 0 
30050 END 
31000 ' 

31005 ' Initialize 
31010 ' 

31015 OPEN "I",#l, "PDS.FIL" 

31016 INPUT #1,PGS;IF PGSO"SUMMARYE" THEN 31016 

31017 INPUT #1,ICFLG,IS,IE,DXI$,DX0S,FS$,FT$,MSG$ 

31020 IFILES=0 

31025 IFILES= IFILES+1 I INPUT #1 ,FLS(IFILES) I IF EOF <1) THEN 
31040 ELSE 31030 

31030 IF FLS < IFILES) = ‘’END’* THEN IFILES = IFILES-1 1 GOTO 31040 
31035 PRINT IFILES; FLS< IFILES) , I GOTO 31025 
31040 CLOSE #1 iPRINT IFILES;FLS(IFILES) 

31050 RETURN 
32000 ' 

32005 ' Time and date routine for SDS MV5.0c 
10/15/84 11140 

32007 ' last modified 

11/30/84 18102 
32010 ' init 
32015 ' 
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32020 DIM M < 12) 

32025 FOR 1=1 TO 12: READ M < I ) : NEXT : IF INT (DATE/4 ) *4 = DATE THEN 
M < 2) =29 

32030 DATA 31,28,31,30,31,30,31,31,30,31,30,31 
32035 ' 

32040 ' main program 
32045 ' 

32100 RETURN: 'DATE=PEEK (&H40) :IF DATE=0 THEN GOSUB 32170 
32105 DATES="19" +RIGHTS C STRS (DATE) ,2) 

32120 

TIMES= RIGHTS (STRS (PEEK <&H43) ) ,2) + " : "+RIGHTS (STRS ( PEEK (&H44) ) , 2) *" 
: "+ RIGHTS (STRS (PEEK < &H45) ) ,2) 

32125 DAY=PEEK (&H41> -*-256* (PEEK (S.H42) ) : M0=0 

32135 FOR 1=1 TO 12 : MO=MO+M ( I ) : IF M0> =DAY THEN MO = MO-M ( I ) : GOTO 
32145 

32140 NEXT 
32145 DAY = DAY- 

MO:MO=I :DATES=RIGHTS(STRS(MO) ,2) -*-"/"-*-RIGHTS (STRS (DAY) , 2) +*V + DATE 
S 

32150 PRINT DATES, TIMES 
32160 RETURN 
32165 ' 

32170 RETURN : ' set time 
32175 ' 

32180 INPUT "MONTH, DAY, YEAR ”, MO , DAY , YR : MO=MO-l 

32185 INPUT "HOUR, MINUTE, SECOND ", HR, MIN, SEC 

32190 JDAY = 0 : FOR 1 = 1 TO MO : JDAY = JDAY+M < I ) : NEXT 

32195 JDAY=JDAY+DAY: JDAYH=INT C JDAY/256) : JDAYL= JDAY- JDAYH *256 

32200 POKE &H41,JDAYL:P0KE &H42, JDAYH 

32205 POKE &H43, HR : POKE &H44 , MIN : POKE S.H45.SEC 

32210 PRINT "depress <CR> to start clock " ; : XCS= INPUTS< 1 ) 

32215 IF XCSOCHRS < 13) THEN GOTO 32185 ELSE POKE &H40,YR 
32220 RETURN 

33000 PRINT "STARTING TIME = ", TIMES 
33010 OPEN "R" , #1 , DXIS+GS+" . MF” , 4 
33020 OPEN "O" , #2 , DXIS+GS+ " . TXT" 

33030 FIELD #1,4 AS NS 
33050 FOR L = 1 TO 241 
33060 FOR IQ = 1 TO 34 
33070 GET #1 
33080 NQ < IQ) = CVS (NS) 

33130 NEXT IQ 
33135 WRITE #2, 

NQ < 1 ) ,NQ(2),NQ(3) ,NQ<4> ,NQ<5) ,NQ<6) , NQ (7 ) , NQ <8) ,NQ<9) , NQ ( 10) , NQ ( 1 
1) ,N0(12> ,NQ(13) ,NQ(14) ,NQ(15) ,N0(16) ,N0(17) ,NQ(18) ,NO(19) ,NO(20) 

, NQ (21 ) ,NQ (22) , NQ < 23) , NQ <24 ) , NQ (25) , NQ <26) , NQ ( 27) , NQ < 28) , NQ < 29) , N 
Q (30) , NQ (31 ) ,NQ <32) ,NQ<33) ,NQ<34) 

33140 NEXT L 

33150 CLOSE #l:CLOSE #2: 

33160 OPEN "R" , #1 , N1S+"R" , 4 
33170 OPEN "0",#2,N1S+".TXT" 

33180 FIELD #1,4 AS NS 
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33190 FOR L = 1 TO 49 
33200 FOR IQ = 1 TO 34 
33210 GET #1 
33220 NQ(IQ)= CVS (NS) 
33230 NEXT IQ 
33235 WRITE #2, 


NQ(l),NQ(2),NQ(3),NQ(4),NQ(5),Na(6),NQ(7),Na(8), 

1 ) , NO ( 12) , NQ (13) , NQ < 14) , NQ (15) , NQ < 16) , NQ < 17) , NQ < 

,NQ(21) rN Q(22),NQ<23),NQ(24), N Q(25),NQ<26),NQ(27;;NaV28;; i iQ(29ri 

Q (30) , NQ (31 ) , NQ ( 32) , NQ < 33) , NQ (34) ' 


NQ (9) , NQ ( 10) , NQ ( 1 
18) , NQ ( 19) , NQ ( 20) 


Q 

33240 

33250 

33255 

33260 


NEXT L 
CLOSE #1: CLOSE #2: 

PRINT "ENDING TIME=", TIMES 
RETURN 


j 
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Appendix 9.8.11. PLOTRE.BAS, a program for the AT computer which 
converts the output o£ SUMMARYE . BAS into line printer plots of 
radiation balances. 

3 ' PLOTRE.BAS modified for AT computer and Epson LQ-1000 
printer 

4 ' 

06/25/86 

5 ' PL0T4K Konza Prairie study 
6/14/85 

6 ' PL0T4C removed refs to MOVE.COM for compilation 
10/26/84 

7 ' PL0T4 combines EBPLOT.CMF (ASCOT) for multiple files 
10/22/84 

9 ' PL0T1 1/12/84 minimum fixes to run on SDS 

10 ' from Program PLOT 1/13/82 (Northstar) 

11 ' Added SUBR 9100 to program Okidata 84 printer 1/12/83 

12 ' Added plot vert, height as variable (HP) 11/82 

20 ' last edited 0946 

6/16/85 

80 ' 

90 DEFINT IrWIDTH "LPT1 : " , 255 : XS= "" : GS = "" : N$="" 

100 DIM AS(IOO) ,D<50) , AVG(35) ,FLS(35) ,PS(145) 

120 G0SUB 31000: 'G0SUB 32000 ' Initialization, Time of Day 

routines 

150 IC0UNT=0 

250 

275 VRES=5 * Vert res. in 144ths of an inch 

276 MARG=17 ' margin offset 
280 WP=5 . 65 : WP=INT ( WP*17 . 1 ) +1 
290 HP = 5:HP = I NT (HP* 144/ VRES ) -•- 1 
305 ' 

500 I COUNT = I COUNT + 1 : GS=FLS ( IC0UNT) + “P" : AS = GS 

505 NS=DXIS+GS:G0SUB 8200:B1J< = N:B5J« = M:G0SUB 8100 
510 OPEN ”I",#1,NS:INPUT #1 , 05 , N , M , DS : CLOSE #1 
515 ' G0SUB 32100 

520 XMIN=1 .ie+38:ymin=xmin:xmax=xmin:ymax=xmin 

525 FOR 1 = 1 TO HP : PS ( I ) =SPACES ( WP) : MIDS (PS (I), 1,1)=*’ I" 

526 MIDS (PS ( I ) , WP , 1 ) = " I : NEXT 

530 FOR 1=1 TO WP:MIDS(PS(1) , I , 1 )="-": MIDS (PS (HP) , 1 , 1 )="-": NEXT 
535 OPEN "R" , #2 , DXIS+GS+ " . MF" , Q5 : FIELD #2,05 AS N9S 

623 YMAX=1000:YMIN=- 

1360 :XMAX=24: XMIN=0 : TS=CHRS (27 ) + CHRS ( 83 ) +CHRS (0 ) : VS=CHRS (27 ) + CHRS 
(84) 

624 LPRINT CHRS (27) CHRS (50) ; CHRS (27) CHRS ( 80 ) ; CHRS ( 15 ) ; 

625 LPRINT TAB (5) DATES " " TIMES " file " NS " label: " DS 

626 LPRINT: LPRINT: LPRINT: LPRINT: LPRINT: LPRINT 

627 LPRINT TAB (MARG+27 ) "UNIVERSITY OF WASHINGTON FOREST 
METEOROLOGY" 

628 LPRINT TAB (MARG+27 ) " RADIATION BALANCE PLOT " 

629 LPRINT " 
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i 
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630 LPRINT TAB (MARG+12) "Konza Prairie" TAB ( MARG+62) "o = Net 

radiation CW m" ; TS; "-2" ; VS ; "3 " 

631 LPRINT TAB (MARG+62) "« = Shortwave rad. down CW m";TS;"- 
2 ";VS;"3" 

632 LPRINT TAB (MARG+17) ; :LPRINT USING 

"##./##/ 85" ; VAL (MIDS C AS ,3,2) ) ; VAL (MIDS < AS , 5 , 2) ) ; : LPRINT 
TAB ( MARG + 62) "x = Shortwave rad. up CW m" ; TS; "-2" ; VS; ”3 " 

633 LPRINT TAB(MARG+62) "= = Long wave rad. down CW m";TS;"~ 

2"; VS; "3" 

634 SYS = VAL (MIDS (AS, 2 , 1 ) ) : LPRINT TAB (MARG) ; :ON SYS GOTO 
635 , 642 , 642 , 642 , 642 , 636 , 639 , 640 , 641 

635 LPRINT ** System l: North Facing Slope": GOTO 642 

637 LPRINT " * System 6: Tempe Az demo" ;:GQTO 642 

639 LPRINT " System 7: East Facing Slope": ;G0T0 642 

640 LPRINT " System 8: West Facing Slope": IGOTO 642 

641 LPRINT " System 9: South Facing Slope": GOTO 642 

642 LPRINT TAB (MARG + 62) “ + = Long wave rad. up CW V*;TS; 

2"; VS; "3 ** 

643 LPRINT TAB (MARG + 62) = Diffuse rad. CW^;TS;"~ 

2"; VS; "3" 

644 LPRINT TAB (MARG+62) "t = Air temperature 
C";TS;"o";VS;"C3 " 

645 LPRINT TAB (MARG+62) "w = Wet bulb temperature 
C" ; TS ; "o" ; VS ; ”C3 " 

646 LPRINT 

783 SX= (WP-1) / (XMAX-XMIN) :SY= (HP-1 ) / (YMAX-YMIN) 

795 LPRINT TAB (MARG) "XMIN=" ; XMIN ; " XMAX=" ; XMAX ; " YMIN=" ; YMIN ; " 
YMAX=" ; YMAX 

796 GOSUB 9100 : LPRINT 
798 COUNT = 1 

800 FOR R%=1 TO Bl« 

802 FOR K* = l TO M 

804 GET #2,K*+M»(RX-1) : D ( Kx> =CVS ( N9S > 

806 NEXT: GOSUB lOOOO 

810 NEXT 

820 FOR I =HP TO 1 STEP -1 
825 LPRINT TAB (MARG) PS(I) 

830 NEXT 

835 FOR 1 = 1 TO 3 .'LPRINT : NEXT .'GOSUB 9110 
840 CLOSE :F0R 1=1 TO 10 : LPRINT : NEXT 
850 IF ICOUNT<IFILES GOTO 500 
855 CHAIN "PLOTEE": END 
900 ' 

5000 GOTO 30020 

7200 IF SW=1 THEN SW=0: RETURN 

7201 IF SX=0 THEN 7250 

7202 IX= (X-XMIN) *SX+1 : IY= ( Y-YMIN) »SY+1 
7205 IF (IX<=0 OR IY<=0) THEN RETURN 
7210 IF (IX>WP OR IY>HP) THEN RETURN 
7215 MIDS (PS ( IY ) , IX, 1 ) =PCS : RETURN 
7250 IF X>EXMX THEN EXMX=X 

7255 IF X<EXMN THEN EXMN=X 


I 
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7260 IF Y>EYMX THEN EYMX= Y 
7265 IF Y <EYMN THEN EYMN = Y 
7270 RETURN 
8090 ' 

8100 PRINT "HEADER DATA FOR: ";NS TAB<30> "LABEL: " D© 

8110 PRINT" NUMBER OF CASES: " N TABOO) "NUMBER OF VARIABLES: " 
M: RETURN 
8120 ' 

8200 ' 

8210 OPEN "I",#l ,N©:INPUT #1,Q5,N,M,D© 

8220 FOR J=1+M1 TO M+M1:INPUT #1 , A© < J) : RSET 
SP©=A©(J) :A©(J)=SP©:N£XT J : INPUT #1,2© 

8230 CLOSE #1: RETURN 
8300 ' 

9100 ' Epson LQ-1000 printer set up 

9105 LPRINT CHR©<15) ; JLPRINT CHR© < 27) CHR© < 51 ) CHR© < 7 ) ; 

9107 RETURN 

9110 LPRINT CHR© <27 ) CHR© (80) ; CHR© <27 ) CHR© < 50 ) : RETURN 

10000 ' Plotting/averaging routine RADPLOT.CMF 

10001 ' 9/24/84 Last modified 1629 

10/034/84 

10002 ' 

10003 IF D < 1 ) =0 AND R*=l THEN RETURN 

10005 FOR 1=1 TO B5% 

10006 AVG<I)=D<I)+AVG<I> 

10007 NEXT 

10008 HR=INT <D< 1)/100) 

10009 IF D C 1 ) =0 THEN D<1)=2400 

10010 MIN=D< 1 ) -HR* 100 

10011 X=HR+MIN/60 

10012 IF MIN=0 OR MIN=12 OR MIN=30 OR MIN=42 THEN Y=0:PCS="- 
":GOSUB 7200: Y= -660: GOSUB 7200 : Y= - 1 160 : GOSUB 7200 

10013 NR=NR+1 

10014 IF MIN MOD 30 = O THEN GOTO 10017 

10015 RETURN 

10016 ' 

10017 FOR 1=1 TO B5* 

10018 AVG < I ) = AVG < I ) /NR 

10019 NEXT 

10020 IF X=0 AND REC0RD*>2 THEN X=24 

10021 Y=AVG<6) : PC©= "*": GOSUB 7200 

10022 Y = AVG<7) :PC©="x".*GOSUB 7200 

10023 Y=AVG<9) : PC©=" =": GOSUB 7200 

10024 Y=AVG<10) :PC©="+" :GOSUB 7200 

10025 Y= AVG (8) : PC©= "#" : GOSUB 7200 

10026 Y = AVG<2) :PC© = "o" .'GOSUB 7200 

10027 Y=AVG<13) *20-1 160 : PC©= "t" : GOSUB 7200 

10028 Y=AVG<14) *20-1160 : PC©= "w" : GOSUB 7200 

10029 FOR 1=1 TO B5% 

10030 AVG < I .) =0 

10031 NEXT 

10032 NR=0 
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10033 RETURN 
30000 ' 

30020 CLOSE 

30030 IF ISTP=0 THEN STOP ELSE CHAIN "SC" 

30040 STOP 
31000 ' 

31005 ' Initialize 
31010 ' 

31015 OPEN "I",#l , "PDS.FIL" 

31016 INPUT #1,PGS:IF PGSO"PLOTRE" THEN 31016 

31017 INPUT #1,ICFLG, IS , IE,DXIS, DXOS , FSS , FTS , MSGS 
31020 IFILES=0 

31025 IFILES=IFILES + 1, : INPUT #1 , FLS ( IFILES) : IF EOF <1> THEN 
31040 ELSE 31030 

31030 IF FLS<IFILES)="END" THEN IFILES=IEILES-1 :GOTO 31040 

31035 PRINT IFILES; FLS ( IFILES) , :GOTO 31025 

31040 CLOSE #1 .'PRINT IFILES ;FLS < IFILES ) ' , 

31050 RETURN 
32000 ' 

32005 ' Time and date routine for SDS MV5.0c 
10/15/84 11:40 

32007 ' last modified 11/ 

1/84 7:50 

32010 ' init 
32015 ' 

32020 DIM M < 12 ) 

32025 FOR 1 = 1 TO 12: READ M<I).'NEXT:IF INT < DATE/4) *4 = DATE THEN 
M (2) =29 

32030 DATA 31,28,31,30,31,30,31,31,30,31,30,31 
32035 ' 

32040 ' main program 
32045 ' 

32100 DATE=PEEK<&H40> : IF DATE=0 THEN GOSUB 32170 
32105 DATES=" 19"+ RIGHTS (STRS (DATE) ,2) 

32120 

TIMES=RIGHTS (STRS (PEEK (&H43) ) , 2) +" : "♦RIGHTS (STRS (PEEK (S.H44) ) , 2) ♦" 
: "♦RIGHTS ( STRS ( PEEK (&H45) > ,2) 

32125 DAY = PEEK (&H41 ) +256 « (PEEK (&H42) ) : M0 = 0 

32135 FOR 1=1 TO 12 : MO=MO+M ( I ) : IF MO>DAY THEN MO=MO-M ( I ) : GOTO 
32145 

32140 NEXT 
32145 DAY = DAY - 

MO:MO=I :DATES=RIGHTS(STRS(MO> ,2) ♦"/"♦RIGHTS (STRS (DAY) ,2) ♦"/"♦DATE 
S 

32150 PRINT DATES, TIMES 
32160 RETURN 
32165 ' 

32170 ' set time 
32175 ' 

32180 INPUT "MONTH, DAY, YEAR ", MO , DAY , YR : M0=M0-1 
32185 INPUT "HOUR, MINUTE, SECOND ", HR, MIN, SEC 
32190 JDAY = 0 : FOR 1 = 1 TO MO : JDAY = JDAY+M < I ) : NEXT 
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32195 JDA Y= JDAY+DAY : JDAYH=INT< JDAY/256) : JDAYL= JDAY- JDAYH*256 

32200 POKE &H41 , JDAYL : POKE &H42,JDAYH 

32205 POKE &H43, HR : POKE 6.H44 , MIN : POKE &H45,SEC 

32210 PRINT "depress <CR> to start clock " ; : XCS=INPUT$ < 1 ) 

32215 IF XC$OCHR$< 13) THEN GOTO 32135 ELSE POKE &H40,YR 

32220 RETURN 
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Appendix 9.8.12. PLOTEE.BAS, a program for the AT computer 
which converts the output of SUMMARYE.BAS into line printer plots 
of energy balances. 

3 ' PLOTRE.BAS modified for AT computer and Epson LQ-IOOO 
printer 

4 

06/25/86 

5 ' PL0T5K Konza Prairie study 
6/14/85 

6 ' PL0T5SC removes MOVE.COM references for compilation 
12/18/84 

7 ' PL0T5S combines EBPLOT.CMF (ASCOT) for multiple files 
10/22/84 

9 ' PL0T1 1/12/84 minimum fixes to runSon SDS 

10 ' from Program PLOT 1/13/82 (Northstar) 

11 y Added SUBR 9100 to program Okidata 84 printer 1/12/83 

12 ' Added plot vert, height as variable (HP) 11/82 

13 ' ASCOT: uses already averaged data (" *S“ files) , y, 

14 ' last edited 1012 

06/16/85 

15 ' 

20 'FTS = "S** -> input 30 min averages 

22 ' = " . MF" -> input 6 min averages 

24 'ICFLG =0 -> plot data as is 

26 ' 1 -> recompute EB w/o soil heat storage in top 10 

cm 
28 

80 ON ERROR GOTO 30000 

90 DEFINT I .'WIDTH "LPT1 : " , 255 : XS= ■••• : GS= •••• : NS= 

lOO DIM B*e < 64) , AS ( 100) , BS ( lOO) , D ( 50) , AVG < 35) , FLS C 35 ) , PS < 145 ) 
lio * OPEN “R", #3, “MOVE.COM" : GET #3 : GT = VARPTR ( #3) : PT = GT +11: CLOSE 
#3 

120 GOSUB 31000: GOSUB 32000 ' Initialization, Time of Day 

routines 

150 IC0UNT=0 

272 ' 

275 VRES=5 ' Vert res. in 144ths of an inch 

276 MARG=17 ' margin offset 
280 WP=5 . 65 : WP=INT (WP*17 . 1 ) +1 
290 HP=5 : HP = I NT ( HP»144/VRES) +1 
305 ' 

500 IC0UNT=IC0UNT+1 : GS=FLS ( ICOUNT) *"P“ : A$ = GS 
505 N$=DXIS+G$ : GOSUB 8200 : B1* = N : B5* = M : GOSUB 8100 
510 OPEN “I",#1,NS:INPUT #1 , Q5 , N , M , DS : CLOSE #1 
515 ' GOSUB 32100 

520 XMIN=1 .1E+3S: YMIN=XMIN:XMAX=XMIN: YMAX=XMIN 
525 FOR 1=1 TO 

HP : P£ ( I ) = SPACES ( WP ) :MID$(PS(I) ,1,1) = “ I “ : MIDS (PS (I) , WP, 1 ) = ’* I ** I NEXT 
530 FOR 1 = 1 TO WP:MIDS(PS(1) , I , 1 ) =" -” : MIDS ( PS (HP) , 1 , 1 )="-": NEXT 
610 OPEN “R",#2,DXIS+GS+“.MF“,Q5:FIELD #2,05 AS N9S 
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624 YMAX=1000: YMIN=- 1360: XMAX=24 :XMIN=0: PI =3. 14159 :DPR=57. 2958 

625 LPRINT 

CHRS < 27 ) CHRS < 50 ) ; CHRS (27) CHRS (80) ; CHRS (15) ; :TS=CHRS(27) +CHRS(83) + 
CHRS(O) : VS=CHRS(27) +CHRS(84) 

626 LPRINT TAB (MARG) DATES " " TIMES " file " NS ” label: 

DS 

627 LPRINT:LPRINT:LPRINT:LPRINT 

628 LPRINT TAB(MARG+27) ; “UNIVERSITY OF WASHINGTON FOREST 
METEOROLOGY" 

629 LPRINT TAB (MARG+27 ) ; " ENERGY BALANCE PLOT ” : 

LPRINT 

630 LPRINT TAB < MARG + 12) "Konza Prairie" TAB < MARG + 62) "v = Net 

radiation CW m";TS; "-2" ; VS; "3 " 

631 LPRINT TAB (MARG +62) = Soil Heat Flux CW 

m" ; T$ ; " -2" ; VS; "3 " 

632 LPRINT TAB (MARG+17) ; :LPRINT USING 

"##/##/84" ; VAL < MIDS < AS ,3,2) ) ; VAL ( MIDS < AS, 5,2) ) ; : LPRINT 
TAB < MARG + 62) "o = Sensible Heat Flux CW jn" ; TS ; "-2" ; VS ; "3 " 

633 LPRINT TAB (MARG +62) "+ = Latent Heat Flux CW 

m";T$;"-2"; VS;"3 " 

634 SYS = VAL (MIDS (AS ,2,1)) : LPRINT TAB (MARG) ;: ON SYS GOTO 
635 , 642 , 642 , 642 , 642 , 636 , 639,640 , 641 

635 LPRINT " System l: North Facing Slope": GOTO 642 

636 LPRINT " System 6: Tempe Az demo"; :GOTO 642 

639 LPRINT " System 7: East Facing Slope": GOTO 642 

640 LPRINT " System 8: West Facing Slope" :GOTO 642 

641 LPRINT " System 9: South Facing Slope" : GOTO 642 

642 LPRINT TAB (MARG+62) "x = Wind direction [Degrees] " 

643 LPRINT TAB (MARG+62) "# = Wind speed Cm s“;TS;"- 

1"; VS; “3 " 

644 LPRINT 

783 SX= (WP-1 ) / (XMAX-XMIN) :SY= (HP-1 ) / ( YMAX- YMIN ) 

795 LPRINT TAB (MARG) "XMIN=" ; XMIN; " XMAX=" ; XMAX ; " YMIN= " ; YMIN ; " 
YMAX="; YMAX 

796 GOSUB 9100 : LPRINT 
798 C0UNT=1 

800 FOR R* = l TO BIX 
802 FOR K*=l TO M 

805 GET #2,Kac+M» (RX-1) ;D(K?«)=CVS(N9S) 

810 NEXT: GOSUB lOOOO 

815 NEXT 

820 FOR I =HP TO 1 STEP -1 
825 LPRINT TAB(MARG) PS ( I ) 

830 NEXT 

835 FOR 1=1 TO 3 : LPRINT :NEXT : GOSUB 9110 
840 CLOSE : FOR 1=1 TO 14 : LPRINT : NEXT 
850 IF ICOUNT<IFILES GOTO 500 
900 ' 

5000 GOTO 30020 

7200 IF SW=1 THEN SW=0: RETURN 

7201 IF SX=0 THEN 7250 

7202 IX= ( X-XMIN ) *SX+1 : IY= ( Y-YMIN) *SY+1 



i 
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7205 IF <IX<=0 OR IY<=0) THEN RETURN 

7210 IF <IX>WP OR IY>HP> THEN RETURN 

7211 Q1=Q1+1:S1=S1+X: S2 = S2-*- Y : S3=S3+X*X : S4=S4+X* Y : S5=S5+ Y« Y 
7215 MIDS C PS < IY) , IX, 1 ) =PCS : RETURN 

7250 IF X>EXMX THEN EXMX=X 
7255 IF X<EXMN THEN EXMN=X 
7260 IF Y>EYMX THEN EYMX=Y 
7265 IF Y<EYMN THEN EYMN=Y 
7270 RETURN 
8090 ' 

8100 PRINT “HEADER DATA FOR: ";NS TABOO) “LABEL: " DS 
8110 PRINT"NUMBER OF CASES: ” N TABOO) "NUMBER OF VARIABLES: " 
M: RETURN 
8120 ' 

8200 ' . 

8210 OPEN "1“ , #1 ,NS : INPUT #1,Q5,N,M,DS 

8220 FOR J=1+M1 TO M+M1:INPUT #1 , AS < J) : RSET , , 

SPS=AS<J) : AS < J ) =SPS : NEXT J : INPUT #1,Z$- 
8230 CLOSE #1: RETURN . 1 

8300 ' 'S?l- 

9100 ' Okidata 84A printer aet up 

9105 LPRINT CHRS<15) ; :LPRINT CHRS <27) CHRS <51 ) CHRS C7) 

9107 RETURN 

9110 LPRINT CHRS < 27) CHRS < 80) ; CHRS <27 ) CHRS (50) : RETURN 

10000 * Plotting/averaging routine EBPLOT.CMF 

10001 ' 9/23/84 1536 Last modified 1525 

10/22/84 

10002 ' 

10003 IF D < 1 ) =0 AND R%=1 THEN RETURN 

10004 IF FTS = “S" THEN GOTO 10018 

10005 FOR 1=1 TO B5X 

10006 AVG<I) =D(I) +AVG<I) 

10007 NEXT 

10008 ' 

10009 ' VECTOR AVG WIND DIRECTION 

10010 ' 

10011 A7 = D < 12) /DPR 

10012 A1 = A1 +COS < A7 ) : A2 = A2 + SIN<A7) 

10013 IF A1O0 THEN A3 = ATN < A2/A1 ) ELSE A3 = SGN < A2) »PI/2 

10014 IF SGN < A1 ) <0 THEN A3 = A3 + PI 

10015 IF SGN < A1 ) >0 AND SGN<A2)<0 THEN A3=A3+2*PI 

10016 AVG<12)=A3*DPR 

10017 IF D < 1 ) =0 THEN D<1)=24 

10018 HR=INT<D< 1)/100) 

10019 MIN=D < 1 ) -HR» 100 

10020 X=HR+MIN/60 : Y=0 

10021 IF MIN MOD 15 = 0 THEN Y=0 :PCS="-“ :GOSUB 7200 

10022 IF MIN MOD 15 = O THEN Y= -500 : PCS="-“ : GOSUB 7200 

10023 IF MIN MOD 15 = O THEN Y=-860:PCS="-“ :GOSUB 7200 

10024 IF FTS = “S" THEN GOTO 10031 

10025 NR=NR + 1 : IF MIN MOD 30 = O THEN GOTO 10027 

10026 RETURN 
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10027 ' 

10028 FOR 1=1 TO B5x 

10029 IF I<> 12 THEN D < I ) =AVG < I ) /NR 

10030 NEXT 

10031 IF ICFLG AND 1 THEN GOSUB 10100 

10032 Y=D<5)+D<27) : PCS«= GOSUB 7200 

10033 Y=D<3> :PCa="o": GOSUB 7200 

10034 Y = D < 4 ) : PCS="-*-" : GOSUB 7200 

10035 Y = D < 12) -860 : PCS= "x" : GOSUB 7200 

10036 Y = D<11) »50-1360:PCS = * , #":G0SUB 7200 

10037 Y=D<2> : PCS= "v" : GOSUB 7200 

10038 IF FTS = "S** THEN RETURN 

10039 FOR 1=1 TO B5X 

10040 AVG < I ) =0 

10041 NEXT 

10042 NR=0: A1=0: A2=0 

10043 RETURN 
10100 ' 

10110 BETA=D<3) /D<4) 

10120 D(4)=-<D(2)+D<5) ) / < 1 +BETA) 

10130 D (3) =BETA*D (4) 

10140 D (27) =0 



10180 RETURN 
10190 ' 

30000 ' 

30010 IF ERR051 THEN 30030 
30020 CLOSE 

30025 PRINT: PRINT: STOP 
30030 ON ERROR GOTO O 
30050 STOP 
31000 ' 

31005 ' Initialize 
31010 ' 

31015 OPEN ** I " , #1 , ”PDS . FIL” 

31016 INPUT #1 , PGS : IF PG$O M PL0TEE" THEN 31016 

31017 INPUT #1, ICFLG, IS, IE, DXIS, DXOS, FSS , FTS , MSGS 
31020 IFILES=0 

31025 IFILES=IFILES + 1 : INPUT #1 , FLS ( IFILES) .* IF EOF (1) THEN 
31040 ELSE 31030 

31030 IF FLS(IFILES) ^‘END** THEN IFILES=IFILES-1 :GOTO 31040 
31035 PRINT IFILES ; FLS (IFILES ) , : GOTO 31025 
31040 CLOSE #1 .‘PRINT IFILES;FLS<IFILES) 

31050 RETURN 
32000 ' 

32005 ' Time and date routine for SDS MV5.0c 
10/15/84 11:40 

32007 ' last modified 11/ 

1/84 7:50 

32010 ' init 
32015 ' 

32020 DIM M < 12) 

32025 FOR 1=1 TO 12 : READ M(I):NEXT:IF I NT < DATE/4) «4 = DATE THEN 


e 


Q- 
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M < 2) =29 

32030 DATA 31,28,31,30,31,30,31,31,30,31,30,31 
32035 RETURN ' 

32040 ' main program 
32045 ' 

32100 DATE=PEEK(&H40) :IF DATE=0 THEN GOSUB 32170 
32105 DATES="1 9” + RIGHTS <STR$<DATE> ,2) 
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ORIGINAL PAGE IS 
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HUE 0 H E Gp Rdn Hup D LPn Up U UDI8 litop lulcp libel l»bot Isoil Eatop Eabut PI PE Odn Gup SHbot H Gs 
•J/a2 «J/«2 «J/t2 «J/i2 ij/i2 «3/i2 a.l/<2 sJ/»2 iJ/*2 i/i Pej C C C C C kPi kPi C kPi ij/i2 «J/«2 X i3/«2 

I SUH 15.1 -6.2 -8.9 -0.6 21.9 -1.2 13.6 31.1 -10.9 59.0 -15.1 0.0 0.2 

Ison 16.6 -6.1 -8.5 -1.3 21.3 -1.1 13.1 18.6 -21.1 12.9 -28.2 0.0 -0.3 

«UK -0.9 0.3 -0.3 0.1 0.6 -0.1 0.1 15.6 -16.1 16.2 -16.9 0.0 0.5 


SYS'En 8, JULY 21, lies, OAT A r tin THE EOtliA P8H8IE. RAnSAS 



ORIGINAL PAGE IS 
OF. POOR QUALITY 



jtsuit 6. Ml ?J. 1986, DATA ES0I1 IKE KC.'IEA PSAJAIE, KANSAS 



ORIGINAL' PAGE IS 
OF POOR QUALITY 
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ORIGINAL PAGE IS 
OF, POOR QUALITY 



Tint 0 H £ Gp Kdn Kup D Ldn tup U UDlft Tatop Twtop Tatot Inset Tsoi 1 Eatop Eabot dT dE Qdft Oup RHbot h G$ Hflrec 
»J/a? Ij/i2 Ij/i2 tJ/i2 iJ/e2 cJ/i: «;/i2 tJ/4? c L 2 a/s dej C C C C C kPa kPa C kPa iJ/i2 iJ/a2 X tJ/>2 


SVSIIH 9, JULY 20, 1986, DATA fROrt IHfc KONZ A PRAIRIE, KANSAS 


ORIGINAL PAGE IS 
OF POOR QUALITY. 
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OF POOR QUALITY 



Tint 0 H E Gp Kdn Kup D Idn L.p U UDIft T a top Mop Tabot T wtot Tsoi 1 6a bop Cabot dT dE Odn Qup RHbot n Gs Hflrec 
iJ/a? *J/e? iJ/i2 iJ/i2 iJ/i? *J/i2 *J/*: *]/• 2 »; a2 i/s C C C C C kP* kPa C kPa tJ/«2 iJ/a2 % iJ/*2 


, 1986, DATA FROM THE MIA FRAME, KANSAS 


ORIGINAL PAGE IS 
OP POOR QUALITY 



ORIGINAL PAGE IS 
OF POOR QUALITY 




— . «m o oo irt >o 10 

i — — r».ojuiffl«3<NM'0'0-Or~r>-or«-f^^3a30»0(^aji 
I CO CM — ■ — I—OOOOOO --0000000<~'OK10' 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 1 


NNCMMMMN<MNNWriN«ONOiNM(NMNN(NNNN«NNNNWN 

oooooooocooaaoooooaooooaooooooooo 




I O * - ' 


K> -O (N O 
n M M 


r<03<N<Nr' — -<iniM'<jcooiin‘^'i — • **" -o 



o w» o- 





o o o o o o c? < 
o o o o ci o o ■ 






I o o <=» < 


■ N ~ <n i 


oo'uic^oN-iso^nNn^^M^'* 
-T V o-» in >d 4 -O 4> in lO Ul l/> Ml Ml Ml » ^ 


nnnnnnnnnnmnnpjmnnm 


a> O M3 • 

I KJ fM fO CM « 

I N N N N I 


Cm 


i cm oo ^ • 

I M - O ^ » I 

i cm cm cm cm cm < 








i m> c-~ r- i 


I CM CM CM CM CM « 







-o 03 r- <m o cm eo M <> 

Ol Ml >M o -M M CM CN CM — »o lO CO IO *o Ml -o o o 



TIME 0 H E Gp Kiln Hup D Ldn Lup U tl»!5 tatop Tutcp Tabot Tubot Tsoi l Eatop EaPot (It HE Odd Oup RKbot « 8s HBrpc 
«J/«2 ij/a2 «J/*2 »J/i2 «J/a2 tJ/>2 iJ/«2 aJ/a2 «J/«2 i/s deg C C C C C kPa kPa C kPa ij/d2 til .2 k |J/|2 

TSbll 22.0 -0.5 -11.5 -0.4 32.9 -5.8 11.1 35.7 -42.2 68.4 -47.9 0.0 -0.4 0.01 


StSTEII 9, JULY 24, 1985, BAJA FROli 1HE KOU2S POAIRIE. KM5SS 
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FLUX DENSITY (W/m2) FLUX DENSITY (W/m2) 






FLUX DENSITY (W/m2) 


ENERGY BALANCE at Ashland, KS., Sys. 8 
Matted wheat stubble, July 14, 1986 



ENERGY BALANCE at Ashland, KS., Sys. 8 
Matted wheat stubble, July 15, 1986 
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ENERGY BALANCE at Ashland, KS., Sys. 8 
Matted wheat stubble, July 18, 1986 



TIME (cdst) 



ENERGY BALANCE at Ashland, KS., Sys. 7 
Matted wheat stubble, July 18, 1986 
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FLUX DENSITY (W/m2) 


ENERGY BALANCE at Ashland, KS., Sys. 9 
Matted wheat stubble, July 14, 1986 
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FLUX DENSITY (W/m2) FLUX DENSITY (W/m2) 


ENERGY BALANCE on the Konza Prairie, KS. 
Grass, Sys. 1, July 20, 1986 




ENERGY BALANCE on the Konza Prairie, KS. 
Grass, Sys. 1, July 21, 1986 
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FLUX DENSITY (W/m2) 



ENERGY BALANCE on the Konza Prairie, KS. 
Grass, Sys. 1, July 24, 1986 
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FLUX DENSITY (W/m2) FLUX DENSITY (W/m2) 


ENERGY BALANCE on the Konza Prairie, KS. 
Grass, Sys. 8, July 20, 1 986 



ENERGY BALANCE on the Konza Prairie, KS. 
Grass, Sys. 8, July 21, 1986 
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FLUX DENSITY (W/m2) 
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ENERGY BALANCE on the Konza Prairie, KS. 
Grass, Sys. 9, July 24, 1986 
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